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.

CGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑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)
              agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGCAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCCCGGtc                                                                                                                                       >  2:3827422/1‑150 (MQ=255)
                agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcac                                                                                                                                     >  2:2400052/1‑150 (MQ=255)
                    gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcacccac                                                                                                                                 >  2:4128449/1‑148 (MQ=255)
                    gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGgcacccac                                                                                                                                 >  1:3057933/1‑148 (MQ=255)
                     aaaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcacccaca                                                                                                                                >  1:683817/1‑150 (MQ=255)
                        gggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcacccacaccc                                                                                                                             >  1:977471/1‑150 (MQ=255)
                        gggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcgcccgcaccc                                                                                                                             >  1:2797349/1‑150 (MQ=255)
                         ggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtgccccccaccca                                                                                                                            >  2:301235/1‑149 (MQ=255)
                           cGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCcccaacccaa                                                                                                                          >  2:466972/1‑150 (MQ=255)
                           cGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCACTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCaaa                                                                                                                          >  1:4188783/1‑150 (MQ=255)
                            gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCCCACCCCAAc                                                                                                                         >  2:2619388/1‑150 (MQ=255)
                             ttACCCCAAACTTGTCGGGCGGGACGCCCCGGAAACTACGGAACCTCCGCGCCGGACTAAAACGTATTGTTATGTTTGGTTTTCTCCCCCTCCGTGTGGTAGAGGGCGGGGTTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                        <  2:1634472/150‑1 (MQ=255)
                                  ccACACCTGTGTGGCCGGCCGACCCGGCAAATACTGCCCCGACGCGACGGACTAAAATTTATTCTTATGGTTGGTTTTCCCCCTCCCCCTGTGGGTGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                                   <  2:1958027/150‑1 (MQ=37)
                                  ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGCCACCACCCCAAcacacaa                                                                                                                   >  2:3175504/1‑144 (MQ=37)
                                  ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCCCCCCCCAAcccacaa                                                                                                                   >  1:2304671/1‑146 (MQ=37)
                                  ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTGCCCCGCACCCAAAcccacac                                                                                                                   >  2:4144503/1‑146 (MQ=38)
                                  ccAAACATGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACTAAAAATTATTGTTACGGTTGGTTTTCTCCCTCCCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                                   <  2:3027375/150‑1 (MQ=38)
                                   cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGCCCCCCCACCCGAAcccacgcc                                                                                                                  >  1:3072775/1‑145 (MQ=38)
                                   cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCGCACCCAAAacccaacc                                                                                                                  >  1:2613620/1‑145 (MQ=38)
                                   cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCCCCCCCAAAcccacccc                                                                                                                  >  1:729982/1‑145 (MQ=37)
                                      ccttGTCTGGCCGGACGACGCGGCAATTAGTCCACCGAGGCGCCGTCCGACAAATATTTGGTAGGGTGGTTTTTCTCCCGCTCCGTGTGGTAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTc                                                                                                               <  2:2833880/147‑1 (MQ=255)
                                      actggTCTGGCCGGACGACCCGGAAAATACTGCCCCGACGCGCCGGACGAAAAAGTTTTCTTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTc                                                                                                               <  1:3831493/146‑1 (MQ=37)
                                      acccgccTGGCCGGACCACCCGGAAAATACCGCCCCGCCGCCCCGGACTAAAAAGTATTCGTACGGTGGGTTTTCCCCCTCCCCCGGGGGGAGTGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTc                                                                                                               <  1:3085644/144‑1 (MQ=255)
                                      aCAGGTCTGGCCGGACGCCCCGGCAAATACTGCCCCGCCGCGACGGACTAAAACGTATTGGTAAGGTTGGTTTTCCCCCTCCCCCTGTGGGAGCGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTc                                                                                                               <  2:2910537/150‑1 (MQ=255)
                                        ctGTCGGGCCGGACGACCGGCAAAATCCGGAACCGCCGCGACGTACTAAATTTTTTTGTTTTGTTTGTTTTTCTCCCACCCCGTGTGGGAGCGGGCGGTGTTTCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                             <  2:1086768/149‑1 (MQ=255)
                                        agctCTGGCCGGACGACCCGCAAAATACTGCACCGACGCGACGGACGAAAAATTATTCGTAAGGTTGGTTTTCTCCCTCTCCCGGTGGGAGTGGGCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                             <  2:59404/147‑1 (MQ=255)
                                        aTGTCTGGCCGGACGACCCGGCAAATACTGCACCGAGGCGCCGGACGAAAAAGTATTGGTATGGTTGTTTTTCTCCCTCTCCCTGTGGTAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                             <  2:1687721/150‑1 (MQ=38)
                                        aTGTCTGGCCGGACGACCCGGAAATTCCTCCACCGACGCGCCGGACTAAAATGTATTCGTACGGTTGGTTTTCCCCCTCTCCCGGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                             <  2:3255373/150‑1 (MQ=255)
                                        aTGTCTGGCCGGACCACCCGGCAAATACTCCACCCACGCGACGGACGAAAAAGTATTCGTAAGGTTGTTTTTCTCCCTCTCCCTGGGGGAGAGGCCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                             <  2:2922286/150‑1 (MQ=255)
                                        aTGTCGGGCCGGACGCCCCGGCAAATACGGCACCGACGCGCCGGACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                             <  2:438547/150‑1 (MQ=38)
                                        aTGTCGGGCCGCACGACCCGCCAAATCCCCCACCGCCGCGCCGCCCCAAAAAGTATTCCTACGGTTGGTTTTCTCCCACCCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                             <  2:814044/150‑1 (MQ=25)
                                         ggTCGGGCGGGACCACCCGGCAAATTCTGCACGGACGCGCCGGACGAAAATTTATTCGTACGGTTGGTTTTCTCCCCCCCCGTGTGGGAGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                            <  2:745408/149‑1 (MQ=255)
                                          gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                           <  2:3145354/150‑1 (MQ=255)
                                          gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCGCCCCCCAACCCACAccccccacg                                                                                                           >  1:1801300/1‑144 (MQ=37)
                                             tGGCGGGACGACCCGGAAAATACTGCACCGCCGCGACGGACGAAAAAGTATTCGTACGGTGGGTTTTCTCCCTCGCCCTGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                                                                                        <  2:1335041/150‑1 (MQ=255)
                                             tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCCCACCCAAAACCAATCCCACAAgggcg                                                                                                        >  1:1677964/1‑148 (MQ=37)
                                               gccggccgCCCCGGAAAATACGCCACCGCCGCGCGGTACTAAAATTTTTTGGTATGTTTGGTTGTCCCCCCCTCCGTGTGGGCGGGGGCCGGGGGGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                      <  1:1614152/150‑1 (MQ=255)
                                               gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCACACCCAAACACACACCCACAAcggggcg                                                                                                      >  2:1544559/1‑146 (MQ=37)
                                               gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAACCCACACCCACCAGGGagag                                                                                                      >  2:3438385/1‑150 (MQ=38)
                                               gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCACCCCAACCCCCACCCCCCAggggggg                                                                                                      >  2:3635484/1‑146 (MQ=37)
                                               gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACACTCACCCAATCACACACCCGCAAggggggg                                                                                                      >  2:1544541/1‑146 (MQ=37)
                                                ccGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAATGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                     <  2:3503784/150‑1 (MQ=255)
                                                ccGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                     <  2:2700728/150‑1 (MQ=255)
                                                ccGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                     <  2:1409280/150‑1 (MQ=255)
                                                ccGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACACCCACCCAAACCCACACACCCAAgggggggg                                                                                                     >  1:2872566/1‑145 (MQ=25)
                                                ccGGACGACCCGCCAAATACGGCACCGACGCGCCGGACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCACTCCCTGTGGGATAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                     <  2:3262506/150‑1 (MQ=38)
                                                 cGGACGCCCCGGCAACGACTGCCCCGAGCCGCCGGACGAAAACGTATTCGTACGGTTGGTTTTCCCCCACTCCCGGTGGGAGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                                    <  1:3010141/150‑1 (MQ=255)
                                                   gacgccCCGGCAAATACTGCACCGACGCGCGGGACTAAAAATTTTTCGTACGTTTGTTTTCCTCCCTCGCCCTGTGGTTGTGGCCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                                  <  2:3989231/150‑1 (MQ=255)
                                                   gacgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                                  <  2:1034608/150‑1 (MQ=255)
                                                      gacCCGGCAAATCCTGCACCGACGCGCCGGACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGCGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgggttagg                                                                                               <  2:2343394/150‑6 (MQ=38)
                                                         ccGGCAAATCCTCCCCCGCCGCGCCGGCCGAAAACGTTTTCGTATGTTGGTTTTTCGCCCGCTCCGTGTGGGAGTGGGCCGGGGGGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                            <  1:4128449/150‑1 (MQ=255)
                                                         ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGGCACCCCACCCCAACCCACACCCACAAGGGCGAGGGGACCGCTc                                                                                            >  1:957116/1‑150 (MQ=25)
                                                           ggCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                          <  2:2736759/150‑1 (MQ=255)
                                                                      ccccGACGCGACGTACTAAAATTTTTCCGTCTGTTTGTTTTACACCCTCTCCCTGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAAtttt                                                                             <  1:424201/148‑1 (MQ=255)
                                                                        ccgccgCGAGGCACTAATATGTATTGGTATGGTGGGTTGGCGCCCTCCCCGGGTGGGTGTGGGCGGGGGTGGGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTga                                                                           <  1:57849/150‑1 (MQ=255)
                                                                        ccGACGCTACGCACGAAAAATTATTCTTAAGGTTGGTTTTCTCCCTCTCCCTGTGGTAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTga                                                                           <  1:957421/150‑1 (MQ=38)
                                                                        ccGACGCGACGGCCGAAAAAGTATTCGTACGGTTGGTGTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTga                                                                           <  2:669317/150‑1 (MQ=39)
                                                                          gacgcgacgGACTAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCGCCCTGTGGGCGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTgaga                                                                         <  1:952935/150‑1 (MQ=255)
                                                                          gacgcgacgGACGAAAAAGTTTTCGTACGGTTGGTTTTCTCCCTCCCCCTGGGGGAGCGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTgaga                                                                         <  2:1720248/150‑1 (MQ=39)
                                                                           acgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGGCCCCACCCCCCAACCCACACCCACAAGCGCGGGGGGGCCGCGCCAG‑‑CCCCAAATTTgtgag                                                                        >  1:4098263/1‑150 (MQ=16)
                                                                               gacggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGCACCCACCCCCCCACCAAAAACCCCAAAGGAGAGGGGGCCCCACGAGCCCAAAAATTTTgagaaaa                                                                      >  1:1317486/1‑147 (MQ=12)
                                                                               gacggacgAAAAAGTATTCGTAAGGGTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTCACCCACACCCTAACCCTCTCCCGCCAGGGAGCGGGGACCGATCGCG‑‑CACCAAATTTGAGAGatat                                                                    >  2:3128375/1‑150 (MQ=37)
                                                                                  ggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCAAACCCACTACCTACAGGGGGAGGGGCCCCCTCGCG‑‑CACAAAATTTGGGGGATAtctt                                                                 >  2:311300/1‑150 (MQ=21)
                                                                                          atttATTCGTTTGGTTGGTTTCCCCCCTCCCGGGGGGGGTGGGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTc                                                         <  1:1219389/147‑1 (MQ=255)
                                                                                            gttttCTTATGGTTGGTTTTCTCCCTCTCCGTGGGGGCGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                       <  1:2400052/147‑1 (MQ=255)
                                                                                             tATTCGTATGGTTGGTTTTCGCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                      <  2:1710035/150‑1 (MQ=255)
                                                                                              aTTCTTAAGTTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                     <  1:2045120/150‑1 (MQ=255)
                                                                                                 cGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                  <  2:1774761/150‑1 (MQ=255)
                                                                                                  gTACGGTTGGTTGCCCCCCCCCCCCGGGGGGAGCGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                 <  1:3418600/150‑1 (MQ=255)
                                                                                                   tatggttggttTTCTCCCGCTCCCTGGGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACag                                                <  2:152232/147‑1 (MQ=255)
                                                                                                      gtttggttTTCCCCCTCGCCCTGTGGGAGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaag                                             <  1:3901874/148‑1 (MQ=255)
                                                                                                       gttgtttTTCTCCCTCCCCCTGGGGTAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                            <  2:2482990/150‑1 (MQ=255)
                                                                                                        ttggttTTCTCCCTCTCGGTGTGGGTGTGGGCCGTGGTGATGGCCTCAGCGCTCTCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                           <  1:3393039/150‑1 (MQ=255)
                                                                                                        ttggttTTCTCCCTCTCCGTGTGGGCGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                           <  2:3072775/150‑1 (MQ=255)
                                                                                                        tgggttTCCTCGCTCGCCGGGTGGTAGCGGGCGGGGGTTCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                           <  1:3186946/148‑1 (MQ=255)
                                                                                                          gtttttCTCCCTCCCCCTGTGGGTGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                         <  1:2813077/148‑1 (MQ=255)
                                                                                                           gttTTCTCCCTCTCCGTGGGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAg                                        <  1:780141/150‑1 (MQ=255)
                                                                                                                  ccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                 <  1:2774371/150‑1 (MQ=255)
                                                                                                                    ctctccgtGGGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                               <  1:2939502/150‑1 (MQ=255)
                                                                                                                        ccgggtgGGTGTGGGCGGGGGTGTGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                           <  1:3273145/146‑1 (MQ=255)
                                                                                                                        ccctGTGGGTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                           <  1:301235/150‑1 (MQ=255)
                                                                                                                         cctGGGGGAGAGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                          <  1:1522174/150‑1 (MQ=255)
                                                                                                                             ggggtgtgggCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                      <  1:3633608/143‑1 (MQ=255)
                                                                                                                               ggtgagGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                    <  2:2304671/147‑1 (MQ=255)
                                                                                                                               ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                    <  1:2793826/150‑1 (MQ=255)
                                                                                                                               ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACAACCCAACCCAAAACCACAACGGAGAGGGGACCGATCCAG‑‑CCCAAAATTTTAGATATATATTAATATGTATTGTAAAATAAGAACAGAAAATCGGGCGTTTGAGccg                    >  2:2046446/1‑149 (MQ=255)
                                                                                                                                 cgcgggCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGa                  <  1:3181261/147‑1 (MQ=255)
                                                                                                                                 agtgggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGa                  <  1:2039659/147‑1 (MQ=255)
                                                                                                                                       cgggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg            <  1:2269054/148‑1 (MQ=255)
                                                                                                                                       ccGGGGTGAGGGCATCAGCGCGCACGTTGACCCACAACCAAAAACAAAAACACAAGCGAGAGGGGGCCGATCGAG‑‑CACAAATATTTAGAGATATAATGTTCTTTGTTGTAAAATAAGAAAAGAAAATAGGGCGTTTTGGCCGGAATATTg            >  2:4005530/1‑150 (MQ=255)
                                                                                                                                       ccGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg            <  2:2797349/150‑1 (MQ=255)
                                                                                                                                        cGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc           <  2:3276145/150‑1 (MQ=255)
                                                                                                                                         ggggTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg          <  2:3375369/150‑1 (MQ=255)
                                                                                                                                          gggttcgggCCTCAGCGAGCACGTTCACCCTCACCCTAACCCTCTACCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACCGAAGAACAGTAAATCGGGCTTTTGCGCCTGAATATTGCgg         <  2:3479439/144‑1 (MQ=255)
                                                                                                                                           ggTTAGGTCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga        <  2:2811213/150‑1 (MQ=255)
                                                                                                                                           ggTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga        <  1:2514564/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑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: