Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,581 A→C 17.1% 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,5810AC17.1% 237.2 / 18.9 133intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (40/68);  new base C (4/21);  total (44/89)
Fisher's exact test for biased strand distribution p-value = 5.84e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 133 total aligned reads displayed.

CGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAA  >  NC_000913/164433‑164727
                                                                                                                                                    |                                                                                                                                                     
cGTGAGGTGCAAAAAGAGATGCAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACGGCACCGACGCTACGGACGACAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                      <  2:3899856/150‑1 (MQ=255)
 gTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCACACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTc                                                                                                                                                     >  1:3963338/1‑150 (MQ=255)
 gTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATc                                                                                                                                                     >  1:1676174/1‑150 (MQ=255)
  tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAACAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                    >  1:4102360/1‑150 (MQ=255)
              agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGCAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCCCGG‑tc                                                                                                                                       >  2:3827422/1‑150 (MQ=255)
                agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑tcac                                                                                                                                     >  2:2400052/1‑150 (MQ=255)
                    gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑tcacccac                                                                                                                                 >  2:4128449/1‑148 (MQ=255)
                    gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGG‑gcacccac                                                                                                                                 >  1:3057933/1‑148 (MQ=255)
                     aaaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑tcacccaca                                                                                                                                >  1:683817/1‑150 (MQ=255)
                        gggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑tcacccacaccc                                                                                                                             >  1:977471/1‑150 (MQ=255)
                        gggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑gcgcccgcaccc                                                                                                                             >  1:2797349/1‑150 (MQ=255)
                         ggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑tgccccccaccca                                                                                                                            >  2:301235/1‑149 (MQ=255)
                           cGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑GCACCcccaacccaa                                                                                                                          >  2:466972/1‑150 (MQ=255)
                           cGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCACTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCaaa                                                                                                                          >  1:4188783/1‑150 (MQ=255)
                            gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCCCACCCCAAc                                                                                                                         >  2:2619388/1‑150 (MQ=255)
                             ttACCCCAAACTTGTCGGGCGGGACGCCCCGGAAACTACGGAACCTCCGCGCCGGACTAAAACGTATTGTTATGTTTGGTTTTCTCCCCCTCCGTGTGGTAGAGGGCGGGGTTGCGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAAcc                                                                                                                        <  2:1634472/150‑1 (MQ=255)
                                  ccACACCTGTGTGGCCGGCCGACCCGGCAAATACTGCCCCGACGCGACGGACTAAAATTTATTCTTATGGTTGGTTTTCCCCCTCCCCCTGTGGGTGAGGGCGGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCtctc                                                                                                                   <  2:1958027/150‑1 (MQ=37)
                                  ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TGGCCACCACCCCAAcacacaa                                                                                                                   >  2:3175504/1‑144 (MQ=37)
                                  ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑GGGCCCCCCCCCCAAcccacaa                                                                                                                   >  1:2304671/1‑146 (MQ=37)
                                  ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGG‑TGCCCCGCACCCAAAcccacac                                                                                                                   >  2:4144503/1‑146 (MQ=38)
                                  ccAAACATGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACTAAAAATTATTGTTACGGTTGGTTTTCTCCCTCCCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCtctc                                                                                                                   <  2:3027375/150‑1 (MQ=38)
                                   cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGT‑GCCCCCCCACCCGAAcccacgcc                                                                                                                  >  1:3072775/1‑145 (MQ=38)
                                   cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TGACCCGCACCCAAAacccaacc                                                                                                                  >  1:2613620/1‑145 (MQ=38)
                                   cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑GGGCCCCCCCCCAAAcccacccc                                                                                                                  >  1:729982/1‑145 (MQ=37)
                                      ccttGTCTGGCCGGACGACGCGGCAATTAGTCCACCGAGGCGCCGTCCGACAAATATTTGGTAGGGTGGTTTTTCTCCCGCTCCGTGTGGTAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTc                                                                                                               <  2:2833880/147‑1 (MQ=255)
                                      actggTCTGGCCGGACGACCCGGAAAATACTGCCCCGACGCGCCGGACGAAAAAGTTTTCTTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCGGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTc                                                                                                               <  1:3831493/146‑1 (MQ=37)
                                      acccgccTGGCCGGACCACCCGGAAAATACCGCCCCGCCGCCCCGGACTAAAAAGTATTCGTACGGTGGGTTTTCCCCCTCCCCCGGGGGGAGTGGGCCGGGGTGCGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTc                                                                                                               <  1:3085644/144‑1 (MQ=255)
                                      aCAGGTCTGGCCGGACGCCCCGGCAAATACTGCCCCGCCGCGACGGACTAAAACGTATTGGTAAGGTTGGTTTTCCCCCTCCCCCTGTGGGAGCGGGCCGGGGTGCGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTc                                                                                                               <  2:2910537/150‑1 (MQ=255)
                                        ctGTCGGGCCGGACGACCGGCAAAATCCGGAACCGCCGCGACGTACTAAATTTTTTTGTTTTGTTTGTTTTTCTCCCACCCCGTGTGGGAGCGGGCGGTGTTTCGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                             <  2:1086768/149‑1 (MQ=255)
                                        agctCTGGCCGGACGACCCGCAAAATACTGCACCGACGCGACGGACGAAAAATTATTCGTAAGGTTGGTTTTCTCCCTCTCCCGGTGGGAGTGGGCCGGGGGGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                             <  2:59404/147‑1 (MQ=255)
                                        aTGTCTGGCCGGACGACCCGGCAAATACTGCACCGAGGCGCCGGACGAAAAAGTATTGGTATGGTTGTTTTTCTCCCTCTCCCTGTGGTAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                             <  2:1687721/150‑1 (MQ=38)
                                        aTGTCTGGCCGGACGACCCGGAAATTCCTCCACCGACGCGCCGGACTAAAATGTATTCGTACGGTTGGTTTTCCCCCTCTCCCGGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                             <  2:3255373/150‑1 (MQ=255)
                                        aTGTCTGGCCGGACCACCCGGCAAATACTCCACCCACGCGACGGACGAAAAAGTATTCGTAAGGTTGTTTTTCTCCCTCTCCCTGGGGGAGAGGCCCGGGGTGAGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                             <  2:2922286/150‑1 (MQ=255)
                                        aTGTCGGGCCGGACGCCCCGGCAAATACGGCACCGACGCGCCGGACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                             <  2:438547/150‑1 (MQ=38)
                                        aTGTCGGGCCGCACGACCCGCCAAATCCCCCACCGCCGCGCCGCCCCAAAAAGTATTCCTACGGTTGGTTTTCTCCCACCCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                             <  2:814044/150‑1 (MQ=25)
                                         ggTCGGGCGGGACCACCCGGCAAATTCTGCACGGACGCGCCGGACGAAAATTTATTCGTACGGTTGGTTTTCTCCCCCCCCGTGTGGGAGCGGGCCGGGGTGCGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                            <  2:745408/149‑1 (MQ=255)
                                          gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                           <  2:3145354/150‑1 (MQ=255)
                                          gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑GGCCCCGCCCCCCAACCCACAccccccacg                                                                                                           >  1:1801300/1‑144 (MQ=37)
                                             tGGCGGGACGACCCGGAAAATACTGCACCGCCGCGACGGACGAAAAAGTATTCGTACGGTGGGTTTTCTCCCTCGCCCTGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                                                                                        <  2:1335041/150‑1 (MQ=255)
                                             tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑GGGCCCCCACCCAAAACCAATCCCACAAgggcg                                                                                                        >  1:1677964/1‑148 (MQ=37)
                                               gccggccgCCCCGGAAAATACGCCACCGCCGCGCGGTACTAAAATTTTTTGGTATGTTTGGTTGTCCCCCCCTCCGTGTGGGCGGGGGCCGGGGGGCGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                      <  1:1614152/150‑1 (MQ=255)
                                               gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TGACCCACACCCAAACACACACCCACAAcggggcg                                                                                                      >  2:1544559/1‑146 (MQ=37)
                                               gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCACACCCAAACCCACACCCACCAGGGagag                                                                                                      >  2:3438385/1‑150 (MQ=38)
                                               gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑GGACCCCCACCCCAACCCCCACCCCCCAggggggg                                                                                                      >  2:3635484/1‑146 (MQ=37)
                                               gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑GGACACTCACCCAATCACACACCCGCAAggggggg                                                                                                      >  2:1544541/1‑146 (MQ=37)
                                                ccGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAATGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGG‑TCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                     <  2:3503784/150‑1 (MQ=255)
                                                ccGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                     <  2:2700728/150‑1 (MQ=255)
                                                ccGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                     <  2:1409280/150‑1 (MQ=255)
                                                ccGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑GGACACCCACCCAAACCCACACACCCAAgggggggg                                                                                                     >  1:2872566/1‑145 (MQ=25)
                                                ccGGACGACCCGCCAAATACGGCACCGACGCGCCGGACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCACTCCCTGTGGGATAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                     <  2:3262506/150‑1 (MQ=38)
                                                 cGGACGCCCCGGCAACGACTGCCCCGAGCCGCCGGACGAAAACGTATTCGTACGGTTGGTTTTCCCCCACTCCCGGTGGGAGCGGGCCGGGGTGCGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                                    <  1:3010141/150‑1 (MQ=255)
                                                   gacgccCCGGCAAATACTGCACCGACGCGCGGGACTAAAAATTTTTCGTACGTTTGTTTTCCTCCCTCGCCCTGTGGTTGTGGCCGGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                                  <  2:3989231/150‑1 (MQ=255)
                                                   gacgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                                  <  2:1034608/150‑1 (MQ=255)
                                                      gacCCGGCAAATCCTGCACCGACGCGCCGGACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGCGGGCCGGGGTGCGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgggttagg                                                                                               <  2:2343394/150‑6 (MQ=38)
                                                         ccGGCAAATCCTCCCCCGCCGCGCCGGCCGAAAACGTTTTCGTATGTTGGTTTTTCGCCCGCTCCGTGTGGGAGTGGGCCGGGGGGCGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                            <  1:4128449/150‑1 (MQ=255)
                                                         ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGGCACCC‑CACCCCAACCCACACCCACAAGGGCGAGGGGACCGCTc                                                                                            >  1:957116/1‑150 (MQ=25)
                                                           ggCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                          <  2:2736759/150‑1 (MQ=255)
                                                                      ccccGACGCGACGTACTAAAATTTTTCCGTCTGTTTGTTTTACACCCTCTCCCTGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAAtttt                                                                             <  1:424201/148‑1 (MQ=255)
                                                                        ccgccgCGAGGCACTAATATGTATTGGTATGGTGGGTTGGCGCCCTCCCCGGGTGGGTGTGGGCGGGGGTGGGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTga                                                                           <  1:57849/150‑1 (MQ=255)
                                                                        ccGACGCTACGCACGAAAAATTATTCTTAAGGTTGGTTTTCTCCCTCTCCCTGTGGTAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTga                                                                           <  1:957421/150‑1 (MQ=38)
                                                                        ccGACGCGACGGCCGAAAAAGTATTCGTACGGTTGGTGTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTga                                                                           <  2:669317/150‑1 (MQ=39)
                                                                          gacgcgacgGACTAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCGCCCTGTGGGCGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTgaga                                                                         <  1:952935/150‑1 (MQ=255)
                                                                          gacgcgacgGACGAAAAAGTTTTCGTACGGTTGGTTTTCTCCCTCCCCCTGGGGGAGCGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTgaga                                                                         <  2:1720248/150‑1 (MQ=39)
                                                                           acgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGT‑GGCCCCACCCCCCAACCCACACCCACAAGCGCGGGGGGGCCGCGCCAG‑‑CCCCAAATTTgtgag                                                                        >  1:4098263/1‑150 (MQ=16)
                                                                               gacggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGG‑GCACCCACCCCCCCACCAAAAACCCCAAAGGAGAGGGGGCCCCACGAGCCCAAAAATTTTgagaaaa                                                                      >  1:1317486/1‑147 (MQ=12)
                                                                               gacggacgAAAAAGTATTCGTAAGGGTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGT‑TCACCCACACCCTAACCCTCTCCCGCCAGGGAGCGGGGACCGATCGCG‑‑CACCAAATTTGAGAGatat                                                                    >  2:3128375/1‑150 (MQ=37)
                                                                                  ggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑GCACCCACACCCAAACCCACTACCTACAGGGGGAGGGGCCCCCTCGCG‑‑CACAAAATTTGGGGGATAtctt                                                                 >  2:311300/1‑150 (MQ=21)
                                                                                          atttATTCGTTTGGTTGGTTTCCCCCCTCCCGGGGGGGGTGGGGGCCGGGGTGCGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTc                                                         <  1:1219389/147‑1 (MQ=255)
                                                                                            gttttCTTATGGTTGGTTTTCTCCCTCTCCGTGGGGGCGTGGGCGGGGGTGCGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                       <  1:2400052/147‑1 (MQ=255)
                                                                                             tATTCGTATGGTTGGTTTTCGCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                      <  2:1710035/150‑1 (MQ=255)
                                                                                              aTTCTTAAGTTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                     <  1:2045120/150‑1 (MQ=255)
                                                                                                 cGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                  <  2:1774761/150‑1 (MQ=255)
                                                                                                  gTACGGTTGGTTGCCCCCCCCCCCCGGGGGGAGCGGGCCGGGGTGCGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                 <  1:3418600/150‑1 (MQ=255)
                                                                                                   tatggttggttTTCTCCCGCTCCCTGGGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACag                                                <  2:152232/147‑1 (MQ=255)
                                                                                                      gtttggttTTCCCCCTCGCCCTGTGGGAGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaag                                             <  1:3901874/148‑1 (MQ=255)
                                                                                                       gttgtttTTCTCCCTCCCCCTGGGGTAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                            <  2:2482990/150‑1 (MQ=255)
                                                                                                        ttggttTTCTCCCTCTCGGTGTGGGTGTGGGCCGTGGTGATGGCCTCAGCGCTCTCGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                           <  1:3393039/150‑1 (MQ=255)
                                                                                                        ttggttTTCTCCCTCTCCGTGTGGGCGTGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                           <  2:3072775/150‑1 (MQ=255)
                                                                                                        tgggttTCCTCGCTCGCCGGGTGGTAGCGGGCGGGGGTTCGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                           <  1:3186946/148‑1 (MQ=255)
                                                                                                          gtttttCTCCCTCCCCCTGTGGGTGAGGGCGGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                         <  1:2813077/148‑1 (MQ=255)
                                                                                                           gttTTCTCCCTCTCCGTGGGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAg                                        <  1:780141/150‑1 (MQ=255)
                                                                                                                  ccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                 <  1:2774371/150‑1 (MQ=255)
                                                                                                                    ctctccgtGGGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                               <  1:2939502/150‑1 (MQ=255)
                                                                                                                        ccgggtgGGTGTGGGCGGGGGTGTGGCCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                           <  1:3273145/146‑1 (MQ=255)
                                                                                                                        ccctGTGGGTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                           <  1:301235/150‑1 (MQ=255)
                                                                                                                         cctGGGGGAGAGGGCGGGGGTGCGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                          <  1:1522174/150‑1 (MQ=255)
                                                                                                                             ggggtgtgggCGGGGGTGCGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                      <  1:3633608/143‑1 (MQ=255)
                                                                                                                               ggtgagGGCCGGGGTGCGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                    <  2:2304671/147‑1 (MQ=255)
                                                                                                                               ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                    <  1:2793826/150‑1 (MQ=255)
                                                                                                                               ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑GCACCCACAACCCAACCCAAAACCACAACGGAGAGGGGACCGATCCAG‑‑CCCAAAATTTTAGATATATATTAATATGTATTGTAAAATAAGAACAGAAAATCGGGCGTTTGAGccg                    >  2:2046446/1‑149 (MQ=255)
                                                                                                                                 cgcgggCCGGGGTGCGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGa                  <  1:3181261/147‑1 (MQ=255)
                                                                                                                                 agtgggCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGa                  <  1:2039659/147‑1 (MQ=255)
                                                                                                                                       cgggggTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg            <  1:2269054/148‑1 (MQ=255)
                                                                                                                                       ccGGGGTGAGGGCATCAGCGCGCACGT‑TGACCCACAACCAAAAACAAAAACACAAGCGAGAGGGGGCCGATCGAG‑‑CACAAATATTTAGAGATATAATGTTCTTTGTTGTAAAATAAGAAAAGAAAATAGGGCGTTTTGGCCGGAATATTg            >  2:4005530/1‑150 (MQ=255)
                                                                                                                                       ccGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg            <  2:2797349/150‑1 (MQ=255)
                                                                                                                                        cGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc           <  2:3276145/150‑1 (MQ=255)
                                                                                                                                         ggggTGTGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg          <  2:3375369/150‑1 (MQ=255)
                                                                                                                                          gggttcgggCCTCAGCGAGCACGT‑TCACCCTCACCCTAACCCTCTACCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACCGAAGAACAGTAAATCGGGCTTTTGCGCCTGAATATTGCgg         <  2:3479439/144‑1 (MQ=255)
                                                                                                                                           ggTTAGGTCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga        <  2:2811213/150‑1 (MQ=255)
                                                                                                                                           ggTGCGGGCCTCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga        <  1:2514564/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGT‑TCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAA  >  NC_000913/164433‑164727

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