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.

CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAG  >  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)
               agatATGAACGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGAAAATACGGCACCGACGCGACGGACTAAAAAGTATGCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑tc                                                                                                                                      <  2:15292201/150‑1 (MQ=38)
                gagaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAACAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgcg                                                                                                                                      >  2:1775914/1‑147 (MQ=255)
                 agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtggac                                                                                                                                     >  1:8065917/1‑146 (MQ=255)
                  ttttaatgggCGTTACCCAAAACCTGTCGGGCCGGACGACCCGGCAAATACTCCACCGACGCGCCGGACGAAAATTTATTGTTATGGTTGTTTTTCTCCCCCCCCCTGTGGGTGTGGGCGGGGGTGAGGCCATCAGCGCGCACGT‑tcacc                                                                                                                                   <  2:10689592/143‑1 (MQ=255)
                  tttGAAAGGGCGTTATCCCAAACTTGTCTGGCCGGACGACCGGGCAAATAGTGCCCCGAGGCGAGGTACTAAAAAGTTTTCTTATGTTTGTTTTTCTCCATCTCCGTGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGT‑tcacc                                                                                                                                   <  2:8321322/148‑1 (MQ=255)
                  tattaatggggGTTACCCCAAACTTGTCGGGCGGGACGCCCCGGAAAATACGGCCCCCCCGCGCGGCACTAAATTTTTTTGTTATGGTTTTTTTCCCCCCCCCCCCTGTGGGTGTGGGCGGGGGTGCGGCCATCAGCGCGCACGT‑tcacc                                                                                                                                   <  2:2930006/143‑1 (MQ=255)
                   aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑tcaccc                                                                                                                                  >  2:13552741/1‑150 (MQ=255)
                          ggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACTAAAAAGTATTGGTATGGTTGGTTTTCTCCCTCTCCCTGTGGTAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGT‑Tcaccctcaccct                                                                                                                           <  2:8848768/150‑1 (MQ=38)
                          ggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCACTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGTGCGCACGT‑gcacccacaccca                                                                                                                           >  1:11623633/1‑149 (MQ=255)
                            cGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTaa                                                                                                                         >  1:1384269/1‑150 (MQ=255)
                              ttATCCCAAACATGTGGGGCCGGACGACCCGGAAAAAACTGCACCGACGCGACGGACTAAAAATTTTTCTTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCGGGGTTGCGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAAcc                                                                                                                       <  2:1380478/150‑1 (MQ=255)
                              ttATCCCAAACATGTCTGGCCGGACGACCCGGCAATTACTGCACCGACGGGAGGGACGAAAAAGTATTCTTATGGTTGGTTTTCTCCCTCGCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAAcc                                                                                                                       <  2:15105557/150‑1 (MQ=38)
                                aTCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑GGGCCCACACCCAAAcccc                                                                                                                     >  1:13600115/1‑149 (MQ=38)
                                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCCCCCCCCCCCACACCctcac                                                                                                                  >  2:208192/1‑148 (MQ=38)
                                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCAACCTCACCCAAACCCactc                                                                                                                  >  2:3730347/1‑150 (MQ=255)
                                    cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGT‑TCACCCACCCCCCAAcccacgcc                                                                                                                 >  1:8773140/1‑145 (MQ=38)
                                    cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTAGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑GAAACCCCACCCAAAaccacaac                                                                                                                 >  1:10051377/1‑145 (MQ=37)
                                     aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCCCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTccc                                                                                                                <  2:7696848/150‑1 (MQ=255)
                                         aTGTCTGGCCGGCCGCCCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTACGGTTGTTTTTCTCCCGCTCCCTGTGGGAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                            <  2:2369925/150‑1 (MQ=255)
                                         aTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAATGTATTCGTACGTTTGGTTTTCTCCCTCGCCCTGTGGGAGCGGGCGGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                            <  2:3546248/150‑1 (MQ=38)
                                         aTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCGAACCCGCGcccgccc                                                                                                            >  2:10039306/1‑146 (MQ=255)
                                         aTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑GGCCCCCCCCCCAAACCCGCACCCCCaa                                                                                                            >  1:6132080/1‑150 (MQ=37)
                                          tgtggggCGGGACGCCCCGGAAAACAGTCCCCCGCCGCGACGGACCAAAATGTATCCGTACGTTTGTTTTTCCCCCTCTCCGTGTGGGTGTGGGCCGGGGTGTGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                           <  2:5464216/145‑1 (MQ=255)
                                          tGTCTGGCGGGACGACCCGGAAAATACTGCACCGCCGCGCGGGACTAAATTTTTTTGGTTTGTTGGTTTTTCTCCCGCCCCGGGTGGTTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                           <  2:4523037/150‑1 (MQ=255)
                                          tGTCTGGCCGGACTCCCCGGAAAATACGGCACCGACGCCCCGGACTAAAATTTTTTGTTTTGGTTGTTTTTCCCCCTCGCCGTGGGGGAGTGGGCCGGGGGGCGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                           <  1:3589608/150‑1 (MQ=255)
                                          ggtggggCGGGACGCCCCGGCTAATCCGGCCCCGCCGGGAGGGACGAAAATGTTTTGTTATGGTTGTTTTTCGCCCACCCGGTGTGGGTGTGGCGGGGGGTGCGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                           <  2:1539313/145‑1 (MQ=255)
                                          ggTCTGGCGGGACGACCCGGCATATACTGCACCGACGCGACGGCCGAAAATGTATTCGTAGGGTTGGTTGTCTCCCTCTCCGGGTGGGAGAGGGCGGGGGTGCGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                           <  2:679797/149‑1 (MQ=255)
                                            tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGT‑GGACCCACACCCAAACCCACTACCCCAAcgg                                                                                                         >  1:7559419/1‑150 (MQ=37)
                                             cTGGCCGGCCGACCCGGCAAATACTGCCCCGACGCGACGGACGAAAAAGTATACGTATGGTTGGTTTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGga                                                                                                        <  2:1100442/150‑1 (MQ=38)
                                             cTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCACCC‑CACCCCAAACCCCACCCTCACgggg                                                                                                        >  1:11174124/1‑149 (MQ=37)
                                              tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑GGACCCGCAACCAAAACCACACCCACCAgggcg                                                                                                       >  2:1143428/1‑148 (MQ=25)
                                              tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑GCCCCCACACCCAAACCCACACCCCCAAGGgag                                                                                                       >  2:2632735/1‑150 (MQ=37)
                                                gcgggACCACCCGGAAAATACTGCACCGACGCGCCGGACGAAAATGTATCCTTATGGTTGTTTTTCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                     <  1:11109396/147‑1 (MQ=255)
                                                gCCGGAGGACCCGCAAAATCCTGCCCCGCCGCGCCGCCCTAAAATTTTTTGGTTCGGTGGGTTTTCTCCCCCCCCGGGTGGGTGCGGGCGGGGGGGCGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTTGAGGGagag                                                                                                     <  2:1384269/150‑1 (MQ=255)
                                                gCCGGACGCCCCGGAAAATACTGCACCGACGCGCCGGACTAAAATTTTTTCGTACGTTTGTTTTTCGCCCCCGCCGTGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                     <  2:4635590/150‑1 (MQ=255)
                                                       gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGG‑TGACCCGCCCCCGAAACCACAACCACAAGGGCGTGGGGGCCg                                                                                              >  1:10257671/1‑150 (MQ=25)
                                                       gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGG‑TCACCCACACCCAAAACCACACCCACAAGGGGGAGGGGGccc                                                                                              >  1:5044091/1‑149 (MQ=37)
                                                             gaaaaTACTGCCCCGACGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCGCCCTGGGGGAGAGGGCGGGGGTGTGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                        <  2:6278142/148‑1 (MQ=255)
                                                                   aCTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGT‑GGCCCCGCACCCGAACCCACACCCACAAGGGGGGGGGGCCCGCTCGGGCACaaa                                                                                  >  1:1759936/1‑150 (MQ=25)
                                                                    gtGCACCTACGCTACGTACTAAAAAGTTTTCGTATGGTTGTTTTTCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAt                                                                                 <  1:12828432/149‑1 (MQ=255)
                                                                     tGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCACACCCCAAACCAAACCCCAAACGGCGAGGGGCCCCATCGAGCACcaaat                                                                                >  1:15270681/1‑148 (MQ=25)
                                                                      gCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAAc                                                                                                                        <  2:9614089/109‑1 (MQ=255)
                                                                       ccccgccgCGACGGCCGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGCGGGCCGGGGTGCGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                              <  2:7506104/148‑1 (MQ=38)
                                                                       cACCGCCGCGCCGCAAGAAAATGTATTCGTATGGTTGGTTTTCTCCCTCTCCGTGTGGGAGGGGGCGGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                              <  1:3049147/150‑1 (MQ=38)
                                                                       cACCGACGCGACGGACGAAAAATTTTTCGTCTGGTTGGTTTTCTCCCTCTCCGTGTGGGTGAGGGCGGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                              <  1:10117301/150‑1 (MQ=38)
                                                                         ccGACGCGCCGGACTAAATTGTTTTCGTATGGTTGTTTTTCCCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                            <  2:13259715/150‑1 (MQ=38)
                                                                                  cggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtct                                                                   <  2:2213046/150‑1 (MQ=255)
                                                                                        aaaatttATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                                             <  1:11377300/144‑1 (MQ=38)
                                                                                         aaaaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑GGCCCCACCCCCAAACCCAATCCCACAATAGCGGGGGGGACGATCGCGCACAAAAATTGATTGATATATTATTATg                                                            >  1:3451362/1‑150 (MQ=12)
                                                                                           attttttCGTATGGTTGGTTTTCTCCCTCTCCGTTTGGGATTGGGCGGGGGTGTGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                          <  2:1179741/145‑1 (MQ=255)
                                                                                           atgTATTCGTTGGGTTGTTTTTCCCCCGCCCGGGGTGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                          <  1:4543404/148‑1 (MQ=255)
                                                                                           actttttCGTCTGGTTGGTTTCCTCCCCCTCCGTGTGGGCGTGGGCCGGGGTGCGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                          <  2:3084833/145‑1 (MQ=255)
                                                                                           acgTAGTCGTAAGGTTGGTTTTCTCCCTCGCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGG‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                          <  2:14252524/148‑1 (MQ=255)
                                                                                              tttgggTAGGGTGGGTTTCCGCCCTCGCCGGGGGGGCGGGGGCGGGGGGGCGGGCCTCAGCGCGCACGG‑GCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                       <  1:12008081/145‑1 (MQ=255)
                                                                                              tagtCGTAAGGTTGGTTTCCGCCCCCCCCCGGTGGGAGCGGCCCGGGGGGCGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                       <  2:382965/147‑1 (MQ=255)
                                                                                               aTTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                      <  2:8103885/150‑1 (MQ=255)
                                                                                                ttCGTATGGTTGTTTTTCTCCCTATCCCTGGGGGAGCGGGCGGGGGTGTGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa                                                     <  1:9173553/150‑1 (MQ=255)
                                                                                                ttCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa                                                     <  2:1235319/150‑1 (MQ=255)
                                                                                                 tcttttggttggttTTCTCCCTCTCCGTGTGGTTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                    <  2:5832386/144‑1 (MQ=255)
                                                                                                   ttagggttggttTCCGCCCTCTCCGTGGGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                  <  1:210041/146‑1 (MQ=255)
                                                                                                   gtatgtttgtttTTCTCCCTCGCCGTGTGGGTGGGGGCCGGGGTGCGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                  <  1:8187721/144‑1 (MQ=255)
                                                                                                   gTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGACCGGGGTGAGGGCATCAGCGCACACGT‑TCACCCACACCCGAACCCACTCCCACAAGGGGGAGGGGACCGATCGAGCACAAATTTTTTGAGATATATTCTTCTGTATTGTAACa                                                  >  1:5273760/1‑150 (MQ=37)
                                                                                                     aaGGTTGGTTTTCTCCCTCTCCCTGGGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                <  1:15389350/150‑1 (MQ=255)
                                                                                                        gttgtttTTCTCCCTCTCCGTGTGGTCGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                             <  1:5503924/150‑1 (MQ=255)
                                                                                                         tgggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                            <  1:9474047/148‑1 (MQ=255)
                                                                                                         gtggttTTCGCCCTCCCCCTGGGGGTGTGGGCCGGGGTGCGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                            <  1:645128/149‑1 (MQ=255)
                                                                                                                 ggccctctcccaCAGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGT‑TCACCCACACCCTAACCCACACCCTCAAGGGGGAGGGGACCCATCCAGCCCAAATTTTGTGAGTTATATTCTTCTGTGTTGTTACAGAAGAACAGaaaaa                                    >  2:11745862/3‑149 (MQ=255)
                                                                                                                  gccctctcccgGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                                   <  1:14880499/149‑1 (MQ=255)
                                                                                                                   ccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑GGACCCGCACCCGAACACTCTACCAAAAGGGAGAGGGGTCCCGTCGAGCACAAATTTTTAGATTTATATTATTCTTTGTATTAACAGAAGAACAGAAAAtag                                  >  2:3621674/1‑148 (MQ=16)
                                                                                                                   ccctccccctGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                  <  2:3349574/150‑1 (MQ=255)
                                                                                                                    cctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGT‑GCCCCCGCACCCGAAACCGCTCCCGCCAGGGGGGGGGGGCCGCTCGGGCACCAATTTTGAGGGATATCTTATTATGTCTTGTTACAGTAGGACAGAAAATCgg                                 >  2:5170107/1‑150 (MQ=19)
                                                                                                                     ctccccctGGGGGAGCGGGCCGGGGGGCGGGCCTCAGCGCGCCCGT‑TCACCCTCACCCTAACCCTCGCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                                <  2:12411959/150‑1 (MQ=255)
                                                                                                                     cactccctGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                                <  1:2632735/148‑1 (MQ=255)
                                                                                                                           ctGGGGGAGCGGGCGGGGGTGCGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTg                          <  2:1891687/150‑1 (MQ=255)
                                                                                                                            tgtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑GGACCCACACCCAAACCCACACCCACAAGGGGGAGGGGCCCGAGCCCGCACAAAGTTTGTGAGATATATTATTGTGTGGTTTAACAGAAGAACATAAAAACggggtgtttg                         >  2:12777529/1‑143 (MQ=255)
                                                                                                                              ggggtgcgggCGGGGGTGCGCGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                       <  1:6677105/143‑1 (MQ=255)
                                                                                                                              ggggcgtgggCGGGGGTGCGGCCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                       <  2:10491560/143‑1 (MQ=255)
                                                                                                                                ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCCCACCCAAAACCACAACCACAAAGGAGAGGGGGCCCATACAGCACAAAAGTTTAGAGATATATTCTAATGTCTTTTAACAGAAAAAAAGAAAATAGGGCTTTTTCGccg                     >  2:12877215/1‑149 (MQ=255)
                                                                                                                                       gCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCGCACCCCAAACCTCGCCCACCCTGGCGCGGGGACCCATAGAGCACAAAATTTGGGAGATATATTCATCGGTTTTGTTAGAGAAGAAAAGCAAAACAGGCTTTTGGGACTGTATAtt              >  2:7372347/1‑150 (MQ=255)
                                                                                                                                        ccGGGGTGAGGCCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg             <  2:6145551/150‑1 (MQ=255)
                                                                                                                                        ccGGGGGGCGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg             <  2:14620384/150‑1 (MQ=255)
                                                                                                                                         cGGGGAGAGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc            <  1:12461192/150‑1 (MQ=255)
                                                                                                                                          ggggTGCGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCTGGCTTTTGCGCCTGAATATTGCg           <  2:4171764/150‑1 (MQ=255)
                                                                                                                                          ggggTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg           <  2:11302830/150‑1 (MQ=255)
                                                                                                                                            ggggcgggCCTCAGCGCGCACGG‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga         <  1:14056666/145‑1 (MQ=255)
                                                                                                                                            ggTGCGGCCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga         <  1:11750620/150‑1 (MQ=255)
                                                                                                                                            ggTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga         <  2:10191159/150‑1 (MQ=255)
                                                                                                                                            ggTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga         <  1:12275009/150‑1 (MQ=255)
                                                                                                                                            cttccacgCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga         <  1:14225456/143‑1 (MQ=255)
                                                                                                                                                  ggCATCCGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTCCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa   <  1:251097/150‑1 (MQ=255)
                                                                                                                                                  ggCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa   <  2:14847711/150‑1 (MQ=255)
                                                                                                                                                  gccATCCGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa   <  1:3980254/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAG  >  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: ‑