Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,581 A→C 18.6% 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,5810AC18.6% 401.4 / 32.2 217intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (59/113);  new base C (2/43);  total (61/156)
Fisher's exact test for biased strand distribution p-value = 2.38e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 225 total aligned reads displayed.

CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCA  >  NC_000913/164432‑164730
                                                                                                                                                     |                                                                                                                                                        
ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                          >  2:2373635/1‑150 (MQ=255)
ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                          >  2:7101457/1‑150 (MQ=255)
ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                          >  1:2471742/1‑150 (MQ=255)
 cGTGCGGTGCAAAACGAGATGAACGGGCGTTATCCCACCCCGGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                         <  2:8270984/150‑1 (MQ=38)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                         >  1:1140257/1‑150 (MQ=255)
   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                       <  2:6165692/150‑1 (MQ=255)
   agAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACAGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGTGCCGGGGTGAGGGCATCa                                                                                                                                                       <  2:7733643/149‑1 (MQ=255)
          cAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCa                                                                                                                                                >  1:206882/1‑150 (MQ=255)
           aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                               >  2:6752315/1‑150 (MQ=255)
            aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                              >  1:2226766/1‑150 (MQ=255)
                  gttGAAAGGGGTTTATCCCAAACTTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTACGTTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTtcacc                                                                                                                                        <  2:1471018/148‑1 (MQ=38)
                  gaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgggac                                                                                                                                        >  1:9009565/1‑145 (MQ=255)
                   aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACggggaccc                                                                                                                                       >  1:7029739/1‑143 (MQ=255)
                   aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgggccc                                                                                                                                       >  2:2840041/1‑144 (MQ=255)
                    tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcacccc                                                                                                                                      >  1:8669138/1‑149 (MQ=255)
                    tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcccccg                                                                                                                                      >  1:1382181/1‑149 (MQ=255)
                      aatgggCTTTATCCAAAACTTTTGGGGCGGGCCCACCCGGAAAATCCTGCCCCGCCGCGACGTACTAAAATTTATTCTTATGGTGGTTTTTCTCCCTCCCCCTGGGGTTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTtcaccctca                                                                                                                                    <  2:9140280/147‑1 (MQ=255)
                      aatgggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGTACGAAAATGTTTTCGTATGGTTGGTTTTCCCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctca                                                                                                                                    <  2:1390663/147‑1 (MQ=38)
                        tgggCTTTATCCCAAACTGGTCTGGCCGGACGACCCGGAAAATACTGCCCCGACGCGACGGACGAAAATGTTTTCGTATGTTTGTTTTTCTCCCTCCCCGTGTGGTTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtcaccctcacc                                                                                                                                  <  2:602528/149‑1 (MQ=255)
                                aTCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAAacca                                                                                                                          >  1:1350903/1‑149 (MQ=255)
                                  cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCGCCCCCACCCAAAaccaca                                                                                                                        >  1:1498456/1‑147 (MQ=38)
                                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGCCCCCCACCCAAACCCcctc                                                                                                                       >  2:1554108/1‑150 (MQ=38)
                                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCACACCCCAAaccaaaa                                                                                                                       >  1:5162349/1‑146 (MQ=37)
                                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCGCCCCCAAACCcacta                                                                                                                       >  1:6064634/1‑149 (MQ=37)
                                    cAAACTTGTGTGGCCGTCGGCCCCGCAAAATACGCCCCCGACGCGCCGGACTAAAATGTTTTCTTATGTTTGTTTTTCCCCCTCCCCCTGTGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                                      <  2:3000493/150‑1 (MQ=255)
                                    cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGAGGGACGAAAAAGTATTCGTATGTTTGTTTTTCTCCCTCTCCCTGTGGGTGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                                      <  2:7012915/150‑1 (MQ=38)
                                    cAAACATGTCGGGCCGGACGCCCCGGCAAATCCGGCCCCGCGGGGACGGACGAAAATGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                                      <  2:4628663/150‑1 (MQ=37)
                                         aTGTCTGCCCGCACCACCCGCCAAATCCTCCACCTACGCGAGTGACGAAAATTTATTCGTAAGGTTGTTTTACACCCTCTCCGTTTGGTTTCGGGCGGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                                 <  1:9177871/150‑1 (MQ=255)
                                          tttCTGGCCGGACGCCCCGCCAAATACTGCACCGACGCGACGTACGAAAAATTATTCGTATGGTTGGTTTTCCCCAGCCCCCTGGGGGATTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                <  2:3499879/148‑1 (MQ=255)
                                          tGTCTGGCGGGCCGACCCGTCAATTACTGCACCGACGCGACGGACTAAAAAGTATTGGAATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                <  1:8668393/150‑1 (MQ=37)
                                            tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCCCGGGCACCCGCCCCCCAACCCACGCCCCCAAggg                                                                                                              >  2:7958758/1‑150 (MQ=37)
                                            tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCGAAACCACACCCACAAcgg                                                                                                              >  2:1379358/1‑150 (MQ=37)
                                              tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCGCACCCCCACCCACACCCACACgggcg                                                                                                            >  1:3004856/1‑148 (MQ=37)
                                                gCCGGACGACCGGCCAAAGTCGGCCCCGCCGCGAGGCACTACAAAGCATTGGTACGGTTGGTTTCCGCCCTCGCCCGGGGGGCGTGGGCGGGGGTGCGGGCCTCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGagag                                                                                                          <  2:4431026/150‑1 (MQ=255)
                                                 ccGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGGGCGCCCGCACCCAAAACCACTCCCACAACGGAGAgg                                                                                                         >  1:9501135/1‑150 (MQ=37)
                                                  cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGGCCCCCACACCCCAAACCGCACCCCCACCGGAGAggg                                                                                                        >  1:9331439/1‑150 (MQ=25)
                                                          ccGGCAAATACTCCACCGCCGCGACTGACCAAAATTTATCGGTATGGTTGGTTTTCCCCCTCCCCGTGTGGGTGGGGGCGGGGGTGGGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                                <  1:8448254/150‑1 (MQ=255)
                                                           cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCCACCCCAAAACCACAACCACAACGGAGAGGGGCCCCATCg                                                                                               >  1:7774483/1‑150 (MQ=25)
                                                           cGGCAAAGACGGCCCCGCCGCGACGGACGAAAACGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGTGCGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCg                                                                                               <  2:647948/150‑1 (MQ=38)
                                                            ggCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCGCACCCAAAACCACACCCACCACGGAGAGGGGGccccaccc                                                                                              >  2:9564891/1‑144 (MQ=25)
                                                             gCAAATCCGGCACCGAGGCGCCGGCCTAAAAAGTTTTCGTATGGTTGGTTTGCTCCCTCCCCCTGTGGGAGTGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                             <  1:8038232/150‑1 (MQ=255)
                                                             gCAAATACTGCACCGCCGCGACGGACGAAAAAGTATTCTTATGGTGGGTTTTCTCCCTCTCCCTGTGGGCGCGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                             <  1:3909219/150‑1 (MQ=38)
                                                             caaaaaACTGCCCCGCCGCGACGGACCAAAATTTTTTCTTTCGGTTGTTTTTCCCCCTCGCCGTGTGGGTTTGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                             <  1:9213046/148‑1 (MQ=255)
                                                               aaaTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGCAACCACACCCAAACCCACACCCCCCAGGGGCGGGGGCCCGCTCgcccc                                                                                           >  2:5891263/1‑146 (MQ=25)
                                                               aaaTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGCCCGCCACCAAAACCACTCCCGCAAGGGTGGGGGGTCCGATCGGG‑‑cc                                                                                         >  1:1640876/1‑149 (MQ=25)
                                                               aaaGACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCCCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑ca                                                                                         <  2:3288669/150‑1 (MQ=255)
                                                                 aTACTGCACCGACGCGACGGACTAAAAATTATTCTTATGGTTGTTTTTCGCCCTCTCCCTGTGGGTGAGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑caca                                                                                       <  1:6254948/150‑1 (MQ=255)
                                                                       cACCGCCGCGAGGGACGAAAATGTTTTCGTACGGTTGGTTTTCCCCCTCCCCGTGTGGTCGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAAtttt                                                                                 <  1:1289881/150‑1 (MQ=255)
                                                                       cACCGACGGGCCGGACGAAAATGTATCGGTAAGTTTGGTTTTCGCCCCCTCCCTGGGGGAGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAAtttt                                                                                 <  2:8402080/150‑1 (MQ=255)
                                                                       cACCGACGCGACGGACGAAAAAGTTTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGTTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAAtttt                                                                                 <  2:4077002/150‑1 (MQ=255)
                                                                         cctccgcgacgGCCGAAAAAGTATTCGTATGTTTGGTTGTCGCCCTCTCCGTGTGGGTGTGGGCGGGGTTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTACCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTga                                                                               <  2:6069977/146‑1 (MQ=255)
                                                                         ccgccgCGCCGGACGAAACAGTATTCGTTCGGTGGGTTTTCTCCCTCTCCCGGTGGGCGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTga                                                                               <  2:1128826/150‑1 (MQ=38)
                                                                         ccGACGCGCCGGACGAAAAAGTATTCGTTAGGTGGGTTTTCTCCCTCGCCCGGTGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTga                                                                               <  1:6369568/150‑1 (MQ=255)
                                                                         ccGACGCGACGGACGAAAAAGTATGCGTAAGGTTGGTTTTCGCCCTCTCCCGGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTga                                                                               <  2:6142938/150‑1 (MQ=255)
                                                                           gacgcgacgGACGAAAAATTTTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTgaga                                                                             <  2:8721978/150‑1 (MQ=255)
                                                                            acgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCAAACCCTCTCCCCCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTtagag                                                                            >  1:4291444/1‑150 (MQ=255)
                                                                                    gacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATAtcttc                                                                    <  2:688931/150‑1 (MQ=255)
                                                                                      cgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCAAAACCACGCCCGCAACGGCGTGGGGCACCATCGAGCCCAAAAATTTTGAGatatattaat                                                                    >  2:9431263/1‑146 (MQ=16)
                                                                                         aaattttttCGTATGGTTGTTTTTCCCCCTCTCCGTGTGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTg                                                               <  1:593501/143‑1 (MQ=255)
                                                                                         aaaaTTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTg                                                               <  1:1230075/150‑1 (MQ=39)
                                                                                         aaaaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCCCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTg                                                               >  2:2285994/1‑150 (MQ=255)
                                                                                           acgTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGGGGGAGCGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTc                                                             <  2:9285968/148‑1 (MQ=255)
                                                                                           acgTATGCGTAAGGGGGGTGTGCGCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTGTGAGAGATATCTTCTTCTGTc                                                             <  2:7940896/148‑1 (MQ=38)
                                                                                              tATTCGTATGGTTGGTGTGCGCCCGCTCCCTGTGGGAGCGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                          <  1:5507498/150‑1 (MQ=39)
                                                                                              tATGCGTATGGTTGGTTTTCTCCCTCCCCCGGGGGGTGAGGGCGGGGGGGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                          <  2:2678192/150‑1 (MQ=255)
                                                                                               tttCGTAAGGTGGGTTGTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                         <  2:7282081/149‑1 (MQ=255)
                                                                                               aTTCTTAAGGTTGGATTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                         <  2:2245014/150‑1 (MQ=255)
                                                                                                  ggTACGGTTGGTTTGCTCCCCCCCCGGTTGGGAGCGGGCGGGGTTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                      <  2:3786579/149‑1 (MQ=255)
                                                                                                   gTATGGTTGTTTTTCTCCCTCTCCGTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                     <  2:5475732/150‑1 (MQ=255)
                                                                                                   gTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaaaa                                                     <  2:9610921/150‑3 (MQ=255)
                                                                                                   gTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                     <  1:7240554/150‑1 (MQ=255)
                                                                                                     atggttggttTGCTCCCTCTCCCGGGGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                   <  2:4999867/148‑1 (MQ=39)
                                                                                                      tggttggttTTCTCCCTCTCCCTGTGGTTGCGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                  <  1:4801186/149‑1 (MQ=255)
                                                                                                           gtttttCCCCATCTCCGTGTGGGTGAGGGGCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCTAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                             <  1:8219719/148‑1 (MQ=39)
                                                                                                           ggttTACCCCCTCTCCCTGTGGGTGAGGGCCGGGGTTCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                             <  1:4886780/150‑1 (MQ=255)
                                                                                                             ttttCTCCCCCTCCCGGTGGTTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGa                                           <  1:6001250/150‑1 (MQ=255)
                                                                                                             tttgcgccctcgcgcgGTGGGAGAGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGa                                           <  2:1237331/144‑1 (MQ=255)
                                                                                                              cccggccctctcccaCAGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCGCACCCAAACCCAAACCCGAAAGGGGGTGGGGGCCGATCGAG‑‑CACAAATTTTGAGTGATATATTATTAAGTCTTGTAAAATaagaaaagaa                                          >  2:2452632/6‑150 (MQ=16)
                                                                                                                 ctccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCCCGGTCACCCTCACCCGAACCCAAACCCTCAAGGGAGTGGGGACCGCTAGAG‑‑CCCGAATTTTGTGAGATATATTCTTCTGTCTTGTAACAGAAGAACAGaaaag                                       >  2:7285021/1‑149 (MQ=25)
                                                                                                                    cctctccggGTGGGCGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                   <  1:2099588/150‑1 (MQ=255)
                                                                                                                    cctctcccgGTGGGTGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                   <  1:9298717/150‑1 (MQ=255)
                                                                                                                    ccacaccctCTGGCAGAGGCCCGGGGTGAGGCCAACACCGCCCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGAACGAG‑‑CACAAATATTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                   <  2:3094224/145‑1 (MQ=255)
                                                                                                                      tctccctGGGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGc                                 <  1:3130269/150‑1 (MQ=255)
                                                                                                                           ctGTGGGAGCGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCACTCCCTCAAGGGAGAGGGGACCGCTCGAG‑‑CACAAATTTGGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAACTCGGGCTTTTg                            <  2:5755095/150‑1 (MQ=255)
                                                                                                                                ggtgagGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                       <  1:5716236/147‑1 (MQ=255)
                                                                                                                                 gcgcgggCGGGGGTGCGGGCCTCCGCGCGCACGTTCACCCTCACCCTACCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTg                      <  1:256531/146‑1 (MQ=255)
                                                                                                                                 gagaGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTg                      <  1:5754414/150‑1 (MQ=255)
                                                                                                                                  agagGGCCGGGGTGAGGGCATCAGCGCGCACGGGGACCCCCAACCCAACCCACAACCACAAGGGGGGGGGGGCCCCTCCCGCCCCAAATTTTGGAGATATAATAATAGGTGTTGTTAAAGAAGAACAGAAAATACGGCCTTTTggcccgg                      >  2:4846065/1‑144 (MQ=255)
                                                                                                                                        ccGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg               <  2:5540377/150‑1 (MQ=255)
                                                                                                                                         cGGGGTGCGGCCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc              <  1:7381562/150‑1 (MQ=255)
                                                                                                                                          gggtTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg             <  2:5790433/150‑1 (MQ=255)
                                                                                                                                          gggtTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg             <  2:5162349/150‑1 (MQ=255)
                                                                                                                                          ggggggTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg             <  1:9626094/150‑1 (MQ=255)
                                                                                                                                          ggggTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg             <  1:7921650/150‑1 (MQ=255)
                                                                                                                                          ggggTGAGGGCATCAGCGCGCACGTTCACCCTCCCCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg             <  2:7884337/150‑1 (MQ=255)
                                                                                                                                          ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg             <  1:4023986/150‑1 (MQ=255)
                                                                                                                                          ggggTGAGGGCATCAGCGCGCACGGTCACCCTCACCCTTACCCTCTCCCGCAAGGGAGAGGGGACCGATCGAG‑‑CACCAATTGTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGACAATCGGGCTTTTGCGCCTGATTATTGCg             <  2:5621513/150‑1 (MQ=255)
                                                                                                                                             gtgtGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGgag          <  2:3654171/150‑1 (MQ=255)
                                                                                                                                              tGAGGGCATCAGCGCGCACGTTCACCCACACCCAAAAACACAACCACAAGGGGGAGGGGACCCAAAGAG‑‑CACAAAATTTTTGTGATATATTATTATGTCTTGTAAAACAAGAAAAAAAAAACGGGCTTTTGCGCCTGAATATTTCggggg         >  2:1079480/1‑147 (MQ=255)
                                                                                                                                              ggcgggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGgaga         <  2:6583727/147‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                        
CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCA  >  NC_000913/164432‑164730

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: