Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,638,058 C→T L455L (CTG→TTG)  creC → sensory histidine kinase CreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,638,0580CT98.2% 169.6 / ‑4.0 55L455L (CTG→TTG) creCsensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (23/30);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.69e-01

AAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGT  >  NC_000913/4637814‑4638299
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
aaagtttcTGGATAACGCCATCGAGTTTACCCCCTAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACATGTCACGCTTAAGGGGCTGTATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGTTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCa                                                                                                                                                                                                                                              <  2:389411/244‑1 (MQ=255)
      aCTGGATAACGTCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTCTCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGc                                                                                                                                                                                                                                        <  1:483093/251‑1 (MQ=255)
       cTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAACTGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGCGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGcg                                                                                                                                                                                                                                       >  2:20087/1‑250 (MQ=255)
             aaCGCCAGAGATGTTACCCCCGCGAGCTGTTGCATCCCGCTAAGCGCCGAAGTGGATCAGGAAAACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTGTGAACTCTTTTTCTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGGTCGTCAGTGAGGTCGCCCGTTTGTTTAACGCCGACGTCACGTTGCGCAACGTTCAGGAAgg                                                                                                                                                                                                                                 <  2:377704/251‑1 (MQ=255)
                         tttACCCCCGATAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGTTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTgg                                                                                                                                                                                                                      <  1:555300/250‑1 (MQ=255)
                                ccGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCtt                                                                                                                                                                                                              <  1:429889/251‑1 (MQ=255)
                                   agagCGGTTGCATCACGCTAAGCGCCGAAGTGGATCCGGAACACGTCACGCTTAAGGTGCTGGATAACGGTAGTGGTATTCCTGACTACGCGCTTTCACGTCTTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGACGCCCGTTTGTTTAACGGCGAAGTCCCGTGGCGCAACTTGCCGGAAGGGGGCGTTCTGGCCACGCTTCg                                                                                                                                                                                                            >  2:156353/1‑250 (MQ=255)
                                   agagCGGTGGCATCACGATAAGGTCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTTGTATGCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTGCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGGCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTc                                                                                                                                                                                                             <  2:517639/249‑1 (MQ=255)
                                        ggTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTc                                                                                                                                                                                                       >  1:440778/1‑250 (MQ=255)
                                                   gCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCCGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTc                                                                                                                                                                                            >  2:113350/1‑250 (MQ=255)
                                                      aaGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACAt                                                                                                                                                                                        <  1:492492/251‑1 (MQ=255)
                                                      aaGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACAt                                                                                                                                                                                        >  1:1360/1‑251 (MQ=255)
                                                       aGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATa                                                                                                                                                                                       <  1:208580/251‑1 (MQ=255)
                                                                             acGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATGTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAg                                                                                                                                                                 <  1:193298/251‑1 (MQ=255)
                                                                              cGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGt                                                                                                                                                                <  1:296114/251‑1 (MQ=255)
                                                                                 cACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCtt                                                                                                                                                             <  1:211705/251‑1 (MQ=255)
                                                                                  aCGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTc                                                                                                                                                            <  1:72025/251‑1 (MQ=255)
                                                                                  aCGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTc                                                                                                                                                            <  1:562833/251‑1 (MQ=255)
                                                                                      ttAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTAt                                                                                                                                                        >  1:198419/1‑251 (MQ=255)
                                                                                       tAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTAt                                                                                                                                                        <  1:561535/250‑1 (MQ=255)
                                                                                          ggTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCctg                                                                                                                                                    >  1:371757/1‑251 (MQ=255)
                                                                                                          aGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCGTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAgg                                                                                                                                     <  2:52260/250‑1 (MQ=255)
                                                                                                             ggTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGgaga                                                                                                                                  >  1:176518/1‑250 (MQ=255)
                                                                                                              gaaTTCATGACTACGCGCTTTAACCGATTTTTCAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAACATCAGCGGACTGGGGTTTGCGTTCGCCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGACACGTTGAGCAACGTGAAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCCCCGTCACTTAACATCGCTTCAAATTCTTCCCACATAGTCTTCGTCTCCTGCTGCCATTGCCAAGgagatg                                                                                                                                >  2:86941/3‑249 (MQ=255)
                                                                                                                   ccTGACTACGCGCTTTCACGTATGTTTGAACGCTTTTACTCTTGGCCTCGTTGAAATGGGCAAAAAAGCAGCGGTCTGGGGTGGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGGCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTAt                                                                                                                          <  2:476616/251‑1 (MQ=255)
                                                                                                                    cTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGTAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATCGCTTCAAATTCTTCCCACATAGTCTTCTTATCCTGCTGCCATTGCAAAGGAGAAGACTATg                                                                                                                         >  2:111416/1‑251 (MQ=255)
                                                                                                                      gACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGACGACTATGtt                                                                                                                       >  1:337295/1‑251 (MQ=255)
                                                                                                                      gACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGt                                                                                                                        <  1:27968/250‑1 (MQ=255)
                                                                                                                              gcTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATccc                                                                                                               <  1:423878/251‑1 (MQ=255)
                                                                                                                                  tCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCTCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTTCAAAGGAGAAGACTATGTTGAAATCCCCCCt                                                                                                           >  2:169273/1‑251 (MQ=255)
                                                                                                                                   cACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCt                                                                                                           <  1:308498/250‑1 (MQ=255)
                                                                                                                                        aTTTTTGAACGCTTTTACTCTTGGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCt                                                                                                      <  2:440778/250‑1 (MQ=255)
                                                                                                                                           tttGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGCCTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTCTCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGa                                                                                                   >  2:111585/1‑250 (MQ=255)
                                                                                                                                             tGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGaaa                                                                                                 >  1:23375/1‑250 (MQ=255)
                                                                                                                                                  gCTTTTACTCTTTGCCTCGTGCAAATGGGCCAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAATGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGCAAATgcct                                                                                           >  2:216418/1‑248 (MQ=255)
                                                                                                                                                        aCTCTTTGCCTCGTGCAAATGGGCAAAAAAGCCGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGCCTAGCCTg                                                                                     >  2:493950/1‑251 (MQ=255)
                                                                                                                                                          tctTTACCTCGAGCGAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTGGTCAGTGCGGCCGCCCGTTTGTTTAACGGCGAAGTCACGATGCGCAACGTGCAGGAAGGAGGCGTGACGGACTCGCTTCGACTTCACGGTACCTTCACAGAGCTTCAAGTTCTTCCCACATAGTATTCGTATCCTGCTGCCATTGCAACGGAGAAAACTCACATTAAATCCCCCATGTTCGGGaaaaagcatagcctagt                                                                                   >  2:531004/1‑238 (MQ=255)
                                                                                                                                                               gCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTg                                                                              <  2:371757/251‑1 (MQ=255)
                                                                                                                                                                          aTGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGTAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                                   >  1:375084/1‑251 (MQ=255)
                                                                                                                                                                             ggCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttg                                                                <  1:20087/251‑1 (MQ=255)
                                                                                                                                                                              gCAAAAAAGAAGCGGGCTGGGTTTGGCGTTCGGCTGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGTCGTGCTGGCCTCGCTTCGACTTCACCGTCCCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGTAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCGGTTTTGTGCAGTAGTGCtgttgt                                                               <  2:518177/251‑1 (MQ=255)
                                                                                                                                                                                    aaGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTc                                                         >  1:353087/1‑251 (MQ=255)
                                                                                                                                                                                     agcagcGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGACCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTcc                                                        >  1:439384/1‑251 (MQ=255)
                                                                                                                                                                                      gccgtccgcTGGGGTTGGTGTTCGGCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAAtt                                                          <  2:486656/240‑1 (MQ=255)
                                                                                                                                                                                              cTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAAtt                                                          >  1:486656/1‑240 (MQ=255)
                                                                                                                                                                                                ggggTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTGAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCGCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTGGGTGCATTATTGCTGTTGTTAATTCCGATAATGc                                                <  2:66631/248‑1 (MQ=255)
                                                                                                                                                                                                     tGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGa                                                                                                                     >  1:215698/1‑174 (MQ=255)
                                                                                                                                                                                                        cGTTCGTCATTGAGGTCGCCCGTTTGTTTAACGTCGAAGTCACTTTGCGCTACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAg                                     <  2:472808/251‑1 (MQ=255)
                                                                                                                                                                                                         gTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGTAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAgg                                    >  1:192060/1‑251 (MQ=255)
                                                                                                                                                                                                                      gTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg                       >  1:254995/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gtttAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTGCGACTTCACCGTCACGTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTAc             <  1:309472/250‑1 (MQ=255)
                                                                                                                                                                                                                                  tttAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAg                                                                           >  1:48425/1‑187 (MQ=255)
                                                                                                                                                                                                                                  tttAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCGGGCCTCGCTTCGACTGCACCGTCACTTCACATAGCTTCAAAGTCTTCCCACAGAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCGGGAAAATGACTAGCCTGTTTGGTGCAg                                                                           <  2:48425/187‑1 (MQ=255)
                                                                                                                                                                                                                                     aaCGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAg        <  1:471578/251‑1 (MQ=255)
                                                                                                                                                                                                                                           gAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGAtgt  <  1:113350/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
AAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGT  >  NC_000913/4637814‑4638299

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: