Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,581 A→C 19.0% 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,5810AC19.0% 543.6 / 59.4 336intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (101/166);  new base C (11/58);  total (112/224)
Fisher's exact test for biased strand distribution p-value = 5.19e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 337 total aligned reads displayed.

CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAG  >  NC_000913/164432‑164728
                                                                                                                                                     |                                                                                                                                                    
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                     >  2:9939794/1‑150 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                     >  2:11109396/1‑150 (MQ=255)
  gTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATc                                                                                                                                                    <  2:1301240/150‑1 (MQ=255)
   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                   >  2:1152799/1‑150 (MQ=255)
   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                   >  1:13294809/1‑150 (MQ=255)
   tGAGGTGAAAAAAGATATGAAAGGGCGTTATCCCAAACATGTATGGCCTGACGACCCGGCAAATACTGCACCTACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTACATCTCCCTGTGGGAGAGGGCCGGGGTGGGGGCATCa                                                                                                                                                   >  1:4866924/1‑150 (MQ=255)
   agAGGCGCAAACAGAGATGAAAGGGCGTTCTCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACACGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                   <  2:8558376/149‑1 (MQ=255)
           aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                           >  2:5064163/1‑150 (MQ=255)
               agatATGAACGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGAAAATACGGCACCGACGCGACGGACTAAAAAGTATGCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTtc                                                                                                                                      <  2:15292201/150‑1 (MQ=38)
                gagaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAACAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgcg                                                                                                                                      >  2:1775914/1‑147 (MQ=255)
                 agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtggac                                                                                                                                     >  1:8065917/1‑146 (MQ=255)
                  ttttaatgggCGTTACCCAAAACCTGTCGGGCCGGACGACCCGGCAAATACTCCACCGACGCGCCGGACGAAAATTTATTGTTATGGTTGTTTTTCTCCCCCCCCCTGTGGGTGTGGGCGGGGGTGAGGCCATCAGCGCGCACGTtcacc                                                                                                                                   <  2:10689592/143‑1 (MQ=255)
                  tttGAAAGGGCGTTATCCCAAACTTGTCTGGCCGGACGACCGGGCAAATAGTGCCCCGAGGCGAGGTACTAAAAAGTTTTCTTATGTTTGTTTTTCTCCATCTCCGTGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTtcacc                                                                                                                                   <  2:8321322/148‑1 (MQ=255)
                  tattaatggggGTTACCCCAAACTTGTCGGGCGGGACGCCCCGGAAAATACGGCCCCCCCGCGCGGCACTAAATTTTTTTGTTATGGTTTTTTTCCCCCCCCCCCCTGTGGGTGTGGGCGGGGGTGCGGCCATCAGCGCGCACGTtcacc                                                                                                                                   <  2:2930006/143‑1 (MQ=255)
                   aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccc                                                                                                                                  >  2:13552741/1‑150 (MQ=255)
                          ggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACTAAAAAGTATTGGTATGGTTGGTTTTCTCCCTCTCCCTGTGGTAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTcaccctcaccct                                                                                                                           <  2:8848768/150‑1 (MQ=38)
                          ggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCACTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGTGCGCACGTgcacccacaccca                                                                                                                           >  1:11623633/1‑149 (MQ=255)
                            cGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTaa                                                                                                                         >  1:1384269/1‑150 (MQ=255)
                              ttATCCCAAACATGTGGGGCCGGACGACCCGGAAAAAACTGCACCGACGCGACGGACTAAAAATTTTTCTTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCGGGGTTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                       <  2:1380478/150‑1 (MQ=255)
                              ttATCCCAAACATGTCTGGCCGGACGACCCGGCAATTACTGCACCGACGGGAGGGACGAAAAAGTATTCTTATGGTTGGTTTTCTCCCTCGCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                       <  2:15105557/150‑1 (MQ=38)
                                aTCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCACACCCAAAcccc                                                                                                                     >  1:13600115/1‑149 (MQ=38)
                                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCCCCCCCACACCctcac                                                                                                                  >  2:208192/1‑148 (MQ=38)
                                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCAACCTCACCCAAACCCactc                                                                                                                  >  2:3730347/1‑150 (MQ=255)
                                    cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTCACCCACCCCCCAAcccacgcc                                                                                                                 >  1:8773140/1‑145 (MQ=38)
                                    cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTAGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGAAACCCCACCCAAAaccacaac                                                                                                                 >  1:10051377/1‑145 (MQ=37)
                                     aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCCCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTccc                                                                                                                <  2:7696848/150‑1 (MQ=255)
                                         aTGTCTGGCCGGCCGCCCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTACGGTTGTTTTTCTCCCGCTCCCTGTGGGAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                            <  2:2369925/150‑1 (MQ=255)
                                         aTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAATGTATTCGTACGTTTGGTTTTCTCCCTCGCCCTGTGGGAGCGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                            <  2:3546248/150‑1 (MQ=38)
                                         aTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCGAACCCGCGcccgccc                                                                                                            >  2:10039306/1‑146 (MQ=255)
                                         aTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCCCCCCCAAACCCGCACCCCCaa                                                                                                            >  1:6132080/1‑150 (MQ=37)
                                          tgtggggCGGGACGCCCCGGAAAACAGTCCCCCGCCGCGACGGACCAAAATGTATCCGTACGTTTGTTTTTCCCCCTCTCCGTGTGGGTGTGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                           <  2:5464216/145‑1 (MQ=255)
                                          tGTCTGGCGGGACGACCCGGAAAATACTGCACCGCCGCGCGGGACTAAATTTTTTTGGTTTGTTGGTTTTTCTCCCGCCCCGGGTGGTTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                           <  2:4523037/150‑1 (MQ=255)
                                          tGTCTGGCCGGACTCCCCGGAAAATACGGCACCGACGCCCCGGACTAAAATTTTTTGTTTTGGTTGTTTTTCCCCCTCGCCGTGGGGGAGTGGGCCGGGGGGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                           <  1:3589608/150‑1 (MQ=255)
                                          ggtggggCGGGACGCCCCGGCTAATCCGGCCCCGCCGGGAGGGACGAAAATGTTTTGTTATGGTTGTTTTTCGCCCACCCGGTGTGGGTGTGGCGGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                           <  2:1539313/145‑1 (MQ=255)
                                          ggTCTGGCGGGACGACCCGGCATATACTGCACCGACGCGACGGCCGAAAATGTATTCGTAGGGTTGGTTGTCTCCCTCTCCGGGTGGGAGAGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                           <  2:679797/149‑1 (MQ=255)
                                            tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGGACCCACACCCAAACCCACTACCCCAAcgg                                                                                                         >  1:7559419/1‑150 (MQ=37)
                                             cTGGCCGGCCGACCCGGCAAATACTGCCCCGACGCGACGGACGAAAAAGTATACGTATGGTTGGTTTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGga                                                                                                        <  2:1100442/150‑1 (MQ=38)
                                             cTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCACCCCACCCCAAACCCCACCCTCACgggg                                                                                                        >  1:11174124/1‑149 (MQ=37)
                                              tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCGCAACCAAAACCACACCCACCAgggcg                                                                                                       >  2:1143428/1‑148 (MQ=25)
                                              tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCACACCCAAACCCACACCCCCAAGGgag                                                                                                       >  2:2632735/1‑150 (MQ=37)
                                                gcgggACCACCCGGAAAATACTGCACCGACGCGCCGGACGAAAATGTATCCTTATGGTTGTTTTTCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                     <  1:11109396/147‑1 (MQ=255)
                                                gCCGGAGGACCCGCAAAATCCTGCCCCGCCGCGCCGCCCTAAAATTTTTTGGTTCGGTGGGTTTTCTCCCCCCCCGGGTGGGTGCGGGCGGGGGGGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGagag                                                                                                     <  2:1384269/150‑1 (MQ=255)
                                                gCCGGACGCCCCGGAAAATACTGCACCGACGCGCCGGACTAAAATTTTTTCGTACGTTTGTTTTTCGCCCCCGCCGTGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                     <  2:4635590/150‑1 (MQ=255)
                                                       gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTGACCCGCCCCCGAAACCACAACCACAAGGGCGTGGGGGCCg                                                                                              >  1:10257671/1‑150 (MQ=25)
                                                       gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCACACCCAAAACCACACCCACAAGGGGGAGGGGGccc                                                                                              >  1:5044091/1‑149 (MQ=37)
                                                             gaaaaTACTGCCCCGACGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCGCCCTGGGGGAGAGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                        <  2:6278142/148‑1 (MQ=255)
                                                                   aCTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGGCCCCGCACCCGAACCCACACCCACAAGGGGGGGGGGCCCGCTCGGGCACaaa                                                                                  >  1:1759936/1‑150 (MQ=25)
                                                                    gtGCACCTACGCTACGTACTAAAAAGTTTTCGTATGGTTGTTTTTCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAt                                                                                 <  1:12828432/149‑1 (MQ=255)
                                                                     tGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCCAAACCAAACCCCAAACGGCGAGGGGCCCCATCGAGCACcaaat                                                                                >  1:15270681/1‑148 (MQ=25)
                                                                      gCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAc                                                                                                                        <  2:9614089/109‑1 (MQ=255)
                                                                       ccccgccgCGACGGCCGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                              <  2:7506104/148‑1 (MQ=38)
                                                                       cACCGCCGCGCCGCAAGAAAATGTATTCGTATGGTTGGTTTTCTCCCTCTCCGTGTGGGAGGGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                              <  1:3049147/150‑1 (MQ=38)
                                                                       cACCGACGCGACGGACGAAAAATTTTTCGTCTGGTTGGTTTTCTCCCTCTCCGTGTGGGTGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                              <  1:10117301/150‑1 (MQ=38)
                                                                         ccGACGCGCCGGACTAAATTGTTTTCGTATGGTTGTTTTTCCCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                            <  2:13259715/150‑1 (MQ=38)
                                                                                  cggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtct                                                                   <  2:2213046/150‑1 (MQ=255)
                                                                                        aaaatttATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                                             <  1:11377300/144‑1 (MQ=38)
                                                                                         aaaaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCACCCCCAAACCCAATCCCACAATAGCGGGGGGGACGATCGCGCACAAAAATTGATTGATATATTATTATg                                                            >  1:3451362/1‑150 (MQ=12)
                                                                                           attttttCGTATGGTTGGTTTTCTCCCTCTCCGTTTGGGATTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                          <  2:1179741/145‑1 (MQ=255)
                                                                                           atgTATTCGTTGGGTTGTTTTTCCCCCGCCCGGGGTGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                          <  1:4543404/148‑1 (MQ=255)
                                                                                           actttttCGTCTGGTTGGTTTCCTCCCCCTCCGTGTGGGCGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                          <  2:3084833/145‑1 (MQ=255)
                                                                                           acgTAGTCGTAAGGTTGGTTTTCTCCCTCGCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                          <  2:14252524/148‑1 (MQ=255)
                                                                                              tttgggTAGGGTGGGTTTCCGCCCTCGCCGGGGGGGCGGGGGCGGGGGGGCGGGCCTCAGCGCGCACGGGCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                       <  1:12008081/145‑1 (MQ=255)
                                                                                              tagtCGTAAGGTTGGTTTCCGCCCCCCCCCGGTGGGAGCGGCCCGGGGGGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                       <  2:382965/147‑1 (MQ=255)
                                                                                               aTTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                      <  2:8103885/150‑1 (MQ=255)
                                                                                                ttCGTATGGTTGTTTTTCTCCCTATCCCTGGGGGAGCGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa                                                     <  1:9173553/150‑1 (MQ=255)
                                                                                                ttCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa                                                     <  2:1235319/150‑1 (MQ=255)
                                                                                                 tcttttggttggttTTCTCCCTCTCCGTGTGGTTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                    <  2:5832386/144‑1 (MQ=255)
                                                                                                   ttagggttggttTCCGCCCTCTCCGTGGGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                  <  1:210041/146‑1 (MQ=255)
                                                                                                   gtatgtttgtttTTCTCCCTCGCCGTGTGGGTGGGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                  <  1:8187721/144‑1 (MQ=255)
                                                                                                   gTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGACCGGGGTGAGGGCATCAGCGCACACGTTCACCCACACCCGAACCCACTCCCACAAGGGGGAGGGGACCGATCGAGCACAAATTTTTTGAGATATATTCTTCTGTATTGTAACa                                                  >  1:5273760/1‑150 (MQ=37)
                                                                                                     aaGGTTGGTTTTCTCCCTCTCCCTGGGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                <  1:15389350/150‑1 (MQ=255)
                                                                                                        gttgtttTTCTCCCTCTCCGTGTGGTCGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                             <  1:5503924/150‑1 (MQ=255)
                                                                                                         tgggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                            <  1:9474047/148‑1 (MQ=255)
                                                                                                         gtggttTTCGCCCTCCCCCTGGGGGTGTGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                            <  1:645128/149‑1 (MQ=255)
                                                                                                                 ggccctctcccaCAGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCACACCCTAACCCACACCCTCAAGGGGGAGGGGACCCATCCAGCCCAAATTTTGTGAGTTATATTCTTCTGTGTTGTTACAGAAGAACAGaaaaa                                    >  2:11745862/3‑149 (MQ=255)
                                                                                                                  gccctctcccgGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                                   <  1:14880499/149‑1 (MQ=255)
                                                                                                                   ccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCGCACCCGAACACTCTACCAAAAGGGAGAGGGGTCCCGTCGAGCACAAATTTTTAGATTTATATTATTCTTTGTATTAACAGAAGAACAGAAAAtag                                  >  2:3621674/1‑148 (MQ=16)
                                                                                                                   ccctccccctGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                  <  2:3349574/150‑1 (MQ=255)
                                                                                                                    cctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGCCCCCGCACCCGAAACCGCTCCCGCCAGGGGGGGGGGGCCGCTCGGGCACCAATTTTGAGGGATATCTTATTATGTCTTGTTACAGTAGGACAGAAAATCgg                                 >  2:5170107/1‑150 (MQ=19)
                                                                                                                     ctccccctGGGGGAGCGGGCCGGGGGGCGGGCCTCAGCGCGCCCGTTCACCCTCACCCTAACCCTCGCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                                <  2:12411959/150‑1 (MQ=255)
                                                                                                                     cactccctGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                                <  1:2632735/148‑1 (MQ=255)
                                                                                                                           ctGGGGGAGCGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTg                          <  2:1891687/150‑1 (MQ=255)
                                                                                                                            tgtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCAAACCCACACCCACAAGGGGGAGGGGCCCGAGCCCGCACAAAGTTTGTGAGATATATTATTGTGTGGTTTAACAGAAGAACATAAAAACggggtgtttg                         >  2:12777529/1‑143 (MQ=255)
                                                                                                                              ggggtgcgggCGGGGGTGCGCGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                       <  1:6677105/143‑1 (MQ=255)
                                                                                                                              ggggcgtgggCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                       <  2:10491560/143‑1 (MQ=255)
                                                                                                                                ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCCCACCCAAAACCACAACCACAAAGGAGAGGGGGCCCATACAGCACAAAAGTTTAGAGATATATTCTAATGTCTTTTAACAGAAAAAAAGAAAATAGGGCTTTTTCGccg                     >  2:12877215/1‑149 (MQ=255)
                                                                                                                                       gCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCCAAACCTCGCCCACCCTGGCGCGGGGACCCATAGAGCACAAAATTTGGGAGATATATTCATCGGTTTTGTTAGAGAAGAAAAGCAAAACAGGCTTTTGGGACTGTATAtt              >  2:7372347/1‑150 (MQ=255)
                                                                                                                                        ccGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg             <  2:6145551/150‑1 (MQ=255)
                                                                                                                                        ccGGGGGGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg             <  2:14620384/150‑1 (MQ=255)
                                                                                                                                         cGGGGAGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc            <  1:12461192/150‑1 (MQ=255)
                                                                                                                                          ggggTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCTGGCTTTTGCGCCTGAATATTGCg           <  2:4171764/150‑1 (MQ=255)
                                                                                                                                          ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg           <  2:11302830/150‑1 (MQ=255)
                                                                                                                                            ggggcgggCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga         <  1:14056666/145‑1 (MQ=255)
                                                                                                                                            ggTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga         <  1:11750620/150‑1 (MQ=255)
                                                                                                                                            ggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga         <  2:10191159/150‑1 (MQ=255)
                                                                                                                                            ggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga         <  1:12275009/150‑1 (MQ=255)
                                                                                                                                            cttccacgCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga         <  1:14225456/143‑1 (MQ=255)
                                                                                                                                                  ggCATCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTCCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa   <  1:251097/150‑1 (MQ=255)
                                                                                                                                                  ggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa   <  2:14847711/150‑1 (MQ=255)
                                                                                                                                                  gccATCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa   <  1:3980254/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAG  >  NC_000913/164432‑164728

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: