Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,581 A→C 20.8% intergenic (+47/‑149) hrpB → / → mrcB putative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913164,5810AC20.8% 445.3 / 57.7 300intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (97/135);  new base C (5/63);  total (102/198)
Fisher's exact test for biased strand distribution p-value = 1.48e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 310 total aligned reads displayed.

CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAA  >  NC_000913/164432‑164727
                                                                                                                                                         |                                                                                                                                                  
ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCa                                                                                                                                                    >  1:9839137/1‑150 (MQ=255)
ccGTGAGGTGCAAAAAGAGATGAAAGGGCGGTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCa                                                                                                                                                    >  2:7982351/1‑150 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCag                                                                                                                                                   >  1:1030053/1‑149 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCAt                                                                                                                                                   >  2:696140/1‑150 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCAt                                                                                                                                                   >  1:11499356/1‑150 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGATATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCag                                                                                                                                                   >  1:8178627/1‑149 (MQ=255)
   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCATCa                                                                                                                                                 >  2:5272754/1‑150 (MQ=255)
     tggTGCAAAAAGTTATGAAAGGGCTTTATCCCAAACTGGTCGGGCGGGACGACCCGGCAAATACTGCACCGACGCGACGGAGTAAAAATTATTGGTATGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGTG‑‑TGGGCGGGGGTGTGGGCATCAgc                                                                                                                                               <  2:1883118/149‑1 (MQ=255)
     aGGTGCAAAAAGAGATGAAAGGGCGTTATCCAAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCATCAgc                                                                                                                                               >  1:10650482/1‑150 (MQ=255)
       gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCATCAgcgc                                                                                                                                             >  2:8823141/1‑150 (MQ=255)
           aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                         >  2:11129548/1‑150 (MQ=255)
              aaGAGTTGAATGGGCTTTATCCCAAACATGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACGAAAATGTATTCGTACGTTTGTTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑TGGGCGGGGGTGTGGGCATCAGCGCGCACGtt                                                                                                                                      <  2:3962364/150‑1 (MQ=37)
                 agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAATTTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGTG‑‑AGGGCCGGGGTGAGGGCCTCAGCGCGCACGTtcac                                                                                                                                   <  2:10220506/150‑1 (MQ=255)
                 agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcac                                                                                                                                   >  2:3966806/1‑150 (MQ=255)
                   aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCATCAGCGCGCACGTggaccc                                                                                                                                 >  2:10897784/1‑150 (MQ=255)
                    tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCATCAGCGCGCACGtgggcccc                                                                                                                                >  1:4172979/1‑143 (MQ=255)
                    tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCATCAGCGCGCACGTggacccc                                                                                                                                >  1:11632140/1‑149 (MQ=255)
                     gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCATCAGCGCGCACGTtgacccac                                                                                                                               >  1:8908398/1‑148 (MQ=255)
                      aatgggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGTGGGCATCAGCGCGCACGTtcaccctca                                                                                                                              <  2:12077726/147‑1 (MQ=255)
                       acgggCGTTATCCCAAACCTCTCGGGCCGGACGCCCCGCCAACTCCTGCACCGCCGCGCCGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCACCCCCTGGGGGAG‑‑TGGGCCGGGGTTAGGGCATCAGCGCGCACGTtcaccctcac                                                                                                                             <  2:10989020/148‑1 (MQ=255)
                              ttATCCCACACTTGTGGGGCCGGCCGCCCCGGAAAATCCTGCCCCGAGGCCCCGGACGAAAAATTTTGCGT‑AGGTTGGTTTT‑CGTCCCTCTCCGGGGGGGTG‑‑TGGGCGGGGGGGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                      <  2:5598800/150‑1 (MQ=255)
                                 ccccAACCTTGTCGGGCGGGACGACCGGGAAAATACTGCCCCGACGCGCGGGACTAAAATTTTTTGTTATGGTTGTTTTT‑C‑CCCCTCGCCGTGGGGGTG‑‑TGGGCGGGGGTGTGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                                                                                   <  2:143747/149‑1 (MQ=255)
                                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCTCCCCCCCCCGAAcccccac                                                                                                                 >  2:3827328/1‑146 (MQ=38)
                                    cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCATCAGCGCGCACGGGCCCCCACACCCAAACCCTCAcc                                                                                                                >  2:9036754/1‑150 (MQ=38)
                                     aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGCCCGCACCCAAAaccaaaccc                                                                                                               >  1:7677908/1‑144 (MQ=37)
                                       aCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCAAAACCACAcccgc                                                                                                             >  1:11090466/1‑148 (MQ=38)
                                        cttgtgtgGCGGGCCGCCCCGGAAAATCGTGCACGGACGCGCGGTACTAAAATGTTTTGGTTTGGTTGGGTGC‑C‑GCCCTCTCCGTGGGGGTG‑‑TGGGCGGGGGTGGGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCa                                                                                                            <  1:2975893/148‑1 (MQ=255)
                                          tGTCTGGCCGGACGCCCGGGCAAATACTGCACGGAGGCGCCGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCCCCCCCTGTGGTAG‑‑TGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                          <  2:6062614/150‑1 (MQ=255)
                                             cTGGCCGGACGACCCGGCAAATACTGCACCGAAGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCATCAGCGCGCACGTGAACCCGCACCCTAACCCACACCCTCAAGGga                                                                                                       >  2:1927931/1‑150 (MQ=38)
                                              tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCCCACCCCAACCCACACCCCCACggggg                                                                                                      >  1:3568363/1‑148 (MQ=37)
                                              gggCCGGACCACCCGCAAAATCCTGCCCCCCCGCGACGGAAAAAAAATTATTCTTAA‑GTTGGTTTTAC‑CCCCACCCCCTGGGGGCG‑‑AGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                                                                                      <  2:2224819/149‑1 (MQ=255)
                                               ggccggccGACCCGCCAAACACGGCCCCGACGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTC‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGCCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                     <  2:322897/150‑1 (MQ=38)
                                                  cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAggg                                                                                                  <  2:9092828/150‑1 (MQ=255)
                                                  cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACCCCCCAAACCCCACCCACAAGGGGGAggg                                                                                                  >  1:5140842/1‑150 (MQ=37)
                                                  cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCAAACCCACACCCGCCCGGGAGAggg                                                                                                  >  1:11236973/1‑150 (MQ=37)
                                                          ccGGCAACTACTGCCCCGCCGCGCCGGACGAAAAAGTTTTCGTCAGGTGGGTTTT‑C‑TCCCGCGCCCGGGGGGAG‑‑TGGGCGGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                          <  2:11914849/150‑1 (MQ=255)
                                                          ccGGCAAATACGGCCCCGACGCGACGGACTAAAAAGTATTCGTCTGGTTGGTTTT‑C‑TCCCTCACCCGGTGGGAG‑‑AGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                          <  2:753558/150‑1 (MQ=38)
                                                          ccGGAAAATACTGCACCGACGCGACGGACGAAAAATTTGTCGTATGGTTGGTTTT‑C‑TCCCCCCCCCTGGGGGTG‑‑AGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                          <  1:2777862/150‑1 (MQ=255)
                                                           cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCCCCCCCCAAACCACGCCCCCCAGGGGGGGGGGGcccgggg                                                                                         >  1:5916338/1‑145 (MQ=21)
                                                             gCAAATACTGCCCCGCGGCGACGCCCGAAAAATTTTGGTTATGTTTGGTTTT‑C‑GCCCCCCCCGTGGGGTTG‑‑TGGGCGGGGGGGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                       <  2:11430917/150‑1 (MQ=255)
                                                             gCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGGCCCCACCCCCACACCCCCGCCCCCCAGGGGGGGGGGGCCGCGCGAg                                                                                       >  1:1894500/1‑150 (MQ=21)
                                                             gCAAATACTCCACCGACGCGACGTACGAAAAAGTATTCGTAAGGTTGGTTTA‑C‑ACCCACTCCCTGTGGGAG‑‑AGGCCCGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                       <  2:7159095/150‑1 (MQ=255)
                                                               aaaTACTGCACCGACGCGACGGACTAAAACTTATTCGTAAGGTTGGTTTT‑C‑TCCCTCGCCCTGTGGGAG‑‑AGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                                                     <  2:1652701/150‑1 (MQ=39)
                                                                         ccgccgCTACGGACTAATATGTATTCGTATGTTTTTTTTC‑C‑GCCCCCCCCGTGTGGTTT‑‑TGGGCCGGGGTGAGGGACTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  2:9286046/150‑1 (MQ=255)
                                                                         ccGAGGCGACGTACAAAAATGTTTGCGTATGGTTGTTTTT‑C‑CCCCTGTCCGTGTGGTAG‑‑TGGGCGGGGGTGTGGGACTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  2:10715191/150‑1 (MQ=255)
                                                                         ccGACGCGCCGGACGAAAATGTATTCGTACGGTTGGTTTT‑C‑TCCCTCGCCGTGGGGGTG‑‑AGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  2:633507/150‑1 (MQ=255)
                                                                         ccGACGCGACGGACGAATAAGTATTCGTACGGTGGGTTTT‑C‑ACCCTCTCCCTGTGGGTG‑‑AGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGATCGATCGAGCACAAATTTTga                                                                           <  1:7885842/150‑1 (MQ=255)
                                                                         ccGACGCGACGGACGAAAAATTATTCTTACGGTTGGTTTT‑C‑TCCCTCTCCCTGGGGTCG‑‑TGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  1:4988455/150‑1 (MQ=255)
                                                                         ccGACGCGACGGACGAAAAAGTAGTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGGG‑‑AGGGCCGGGGTGAGGGCATCAGCGCGCACGTACACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAACTTTga                                                                           <  2:5481216/150‑1 (MQ=255)
                                                                            ccgcgaggGACTAAAATTTATTCGTTCGGTTGGTTTC‑C‑GCCCTCCCCGGGGGGTTG‑‑TGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagag                                                                        <  1:653712/149‑1 (MQ=255)
                                                                            acgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCGAACCCTCGCCCGCCAGGGCGGGGGGGCCGATCGAGCACAAATTTTgcgcg                                                                        >  2:7862340/1‑146 (MQ=37)
                                                                             cgcgacgGACTAAAATGTATTCGTAAGGTTGGTTTT‑C‑TCCCGCTCCCTGTGGGTG‑‑AGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagaga                                                                       <  2:9880893/150‑1 (MQ=255)
                                                                                gacggacgAAAAAGTATTCGTAAGGTTGGTTTA‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCGCACCCAAAACCCCACCCACAACGGAGGGGGGACCGATCGAGCGCAAAGTTTGTGAGatat                                                                    >  1:10218494/1‑150 (MQ=25)
                                                                                   ggacgAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCATCCGCGCGCACGGGCCCCCACCCCCCAACCCACAACCACCAAGGGGAGGGGGCCGGACGAGCCCCAAATTTGGGAGgtaaatt                                                                 >  1:5804550/1‑146 (MQ=12)
                                                                                   ggacgAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCATCAGCGCGCAGGTTCACCCACACCCTAACCCACACCCACCAGGGAGAGGGGACCGCTCGAGCACCAATTTTGAGAGATAtatt                                                                 >  1:5619630/1‑150 (MQ=37)
                                                                                    gacgAAAACTTTTTGGTATGGTTGGTTTT‑C‑GCCCGCGCCCGGTGGGTG‑‑CGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGGTCGAGCACAAATTTTGAGAGATAtcttc                                                                <  2:6473051/150‑1 (MQ=255)
                                                                                    gacgAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCAAACCCAAAACCACTACCACAACGGAGGGGGGGCCGATCGAGCACCAAGTTTGTGAGataagatc                                                                >  2:11947225/1‑145 (MQ=21)
                                                                                      ctaaaaaGTATTCGTACGGTGGGTTTT‑C‑TCCCACTCCCTGTGGGTG‑‑CGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttctt                                                              <  2:6540797/148‑1 (MQ=39)
                                                                                      cgAAACTGTATTCGTACGGTTGGTTTT‑C‑GCCCTCTCCCAGTGGGAG‑‑AGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttctt                                                              <  2:11199299/150‑1 (MQ=38)
                                                                                       gAAAAAGTATTCTTAGGGTTGGTTTT‑C‑ACCCTCTCCCTGGGGTTG‑‑TGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttc                                                             <  2:3643269/150‑1 (MQ=255)
                                                                                        aaaaaGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACCCCCAAACCCACACCCACAAGGGAGAGGGGGCCCAACGAGCCCAAAATTTGATAGATAtattaatat                                                            >  1:11959528/1‑145 (MQ=21)
                                                                                           aaGTATCCGTATGGTTGGTTTT‑C‑TCCCTCCCCCTGTGGGCG‑‑TGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                         <  1:6038116/150‑1 (MQ=39)
                                                                                             ttATTCGTACGGTGGGTTTT‑C‑ACCCTCTCCCGGTGGTTT‑‑TGGGCGGGGGTGTGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                       <  2:4513483/149‑1 (MQ=255)
                                                                                             gTATTGGTACGGTTGGTTTG‑C‑GCCCACGCCGGGTGGGCG‑‑AGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                       <  1:482975/150‑1 (MQ=255)
                                                                                             gTATTCGTATGGTTGGTTTT‑C‑GCCCTAGCCGTGTGGGAG‑‑AGGGCCGGGTTTAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                       <  2:8710535/150‑1 (MQ=255)
                                                                                             gTATTCGTATGGTGGGTTGT‑C‑TCCCTCTCCCGGTGGGAG‑‑TGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                       <  2:7354655/150‑1 (MQ=255)
                                                                                             gTATTCGTACGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑TGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                       <  1:9050560/150‑1 (MQ=255)
                                                                                             gTATTCGTAAGGGTGGTTCT‑C‑TCCCTCTCCCTGTGGTAG‑‑AGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                       <  2:6806458/150‑1 (MQ=255)
                                                                                                 tCGTATGGTTGTTTTT‑C‑TCCCTCTCCCTGTGGTTG‑‑TGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                   <  1:2060768/150‑1 (MQ=255)
                                                                                                 tCGTACGGTTGTTTTT‑C‑TCCCACCCCGTGGGGGTG‑‑AGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                   <  2:11720039/150‑1 (MQ=255)
                                                                                                  gggtatgttggttTT‑‑‑TCCCCCTCCCCGTGGGTGTTGGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                   <  2:10307092/144‑1 (MQ=255)
                                                                                                  cgtctggttggttTT‑C‑TCCCGCTCCCGGGGGGAT‑‑TGGGCCGGTGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                  <  2:484961/145‑1 (MQ=255)
                                                                                                   gTATGGTTGGTTTT‑C‑TCCCTCTCCCGGGGGGAG‑‑TGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCTCCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                 <  2:10147176/150‑1 (MQ=39)
                                                                                                   gTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑TGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                 <  1:3962446/150‑1 (MQ=255)
                                                                                                     ttggtgggttTA‑C‑TCCCTCTCCGGGTGGGTG‑‑TGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                               <  1:10897784/148‑1 (MQ=255)
                                                                                                      cggtgggttTG‑C‑TCCCCCTCCCGGGGGGTG‑‑CGGGCCGGGGGGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                              <  1:1320239/149‑1 (MQ=255)
                                                                                                          tggttTT‑C‑TCCCTCTCCCTGGGGGAG‑‑AGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                          <  2:12234704/150‑1 (MQ=255)
                                                                                                               tt‑C‑TCCCTCTCCCGGTGGGAG‑‑AGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAG‑ACAGaaga                                    <  2:2971646/150‑3 (MQ=255)
                                                                                                                    acccgcccccgGTGGGTG‑‑CGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                                  <  1:6264703/143‑1 (MQ=255)
                                                                                                                     ccctgtccgtGTGGGTG‑‑CGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGTGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                 <  2:11567377/150‑1 (MQ=255)
                                                                                                                     ccctctccgtGGGGGTG‑‑TGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                 <  1:696140/150‑1 (MQ=255)
                                                                                                                     ccctctccctGTGGGAG‑‑AGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                 <  2:9648148/150‑1 (MQ=255)
                                                                                                                     ccctccccgtGGGGGAG‑‑CGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                 <  2:10611566/150‑1 (MQ=255)
                                                                                                                       ctcgccgtGTGGGAG‑‑TGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                               <  1:7957368/150‑1 (MQ=255)
                                                                                                                       cgctccgtGTGGGTG‑‑TGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                               <  2:4373661/148‑1 (MQ=255)
                                                                                                                              tgggggAG‑‑TGGGCGGGGGGTAGGGCCTCCGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc                        <  1:3266447/147‑1 (MQ=255)
                                                                                                                                tgGGAG‑‑TGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                      <  1:7280763/150‑1 (MQ=255)
                                                                                                                                  ggtg‑‑agGGCCGGGGTGCGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                    <  1:10006014/147‑1 (MQ=255)
                                                                                                                                    ggtgcgggCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGAGCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                    <  2:9509025/145‑1 (MQ=255)
                                                                                                                                           gCCGGGGTGAGGGCATCAGCGCGCACGTTCAACCGCACCCGAACCCACAACCACAAGGGAGATGGGGCCCAAAGAGCACAAAATTTGATAGATATATTAATATGTATTTTTAAACAATAAAAGAAAAAAGGGCATTTGAGCCAGAATAtt             >  2:7595808/1‑150 (MQ=255)
                                                                                                                                             gggggTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc           <  1:8871181/149‑1 (MQ=255)
                                                                                                                                             cGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc           <  2:8178783/150‑1 (MQ=255)
                                                                                                                                              ggggTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg          <  2:10005636/150‑1 (MQ=255)
                                                                                                                                              ggggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg          <  1:8478421/150‑1 (MQ=255)
                                                                                                                                              ggggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg          <  2:9799941/150‑1 (MQ=255)
                                                                                                                                                  tGAGGGCATCAGCGCGCACGGGCACCCACACCCAAAACCAAAACCACAAAGGGGATGGGGCCCCTACGGCACAAAATTTGTGAGATAAATTATTCTGTGTTGTAAAAAAACAAAAGAAAAAACGGCGTTTTCGCCGGAATATTTGGgaga      >  2:960105/1‑150 (MQ=255)
                                                                                                                                                    aGGGCATCAGCGCGCACGTTCACCCACACCCAAACCCTAACACTAAAGGGGGAGGCGACCGCTCCAGCACAAAAATTGAGAGATATAATCTTATGTATTTTAACAGAACAACAGAAAAACAGGCGTTTGAGCCTGAATAATTCggcggaa    >  2:5885596/1‑145 (MQ=255)
                                                                                                                                                     gggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa   <  1:2488566/150‑1 (MQ=255)
                                                                                                                                                      ggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa  <  1:4497416/150‑1 (MQ=255)
                                                                                                                                                      cccATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa  <  2:3394527/148‑1 (MQ=255)
                                                                                                                                                         |                                                                                                                                                  
CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑C‑TCCCTCTCCCTGTGGGAG‑‑AGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAA  >  NC_000913/164432‑164727

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