Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,581 A→C 21.5% 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,5810AC21.5% 197.2 / 18.9 121intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (44/48);  new base C (5/23);  total (49/72)
Fisher's exact test for biased strand distribution p-value = 4.62e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 122 total aligned reads displayed.

CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAA  >  NC_000913/164432‑164726
                                                                                                                                                     |                                                                                                                                                 
ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                   >  2:4119438/1‑150 (MQ=255)
   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCACACATGTCTGGCCGGACGACCCGGCAAATACGGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                <  2:577266/150‑1 (MQ=255)
   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                <  2:283497/150‑1 (MQ=255)
   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                >  1:3311776/1‑150 (MQ=255)
       gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCCgcgc                                                                                                                                            >  2:241582/1‑150 (MQ=255)
       gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCGCAgcgc                                                                                                                                            >  1:943765/1‑150 (MQ=255)
       gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgc                                                                                                                                            >  2:2522529/1‑150 (MQ=255)
           aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCCAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTACGGGTGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCGgcgcgcgg                                                                                                                                        >  2:2769001/1‑148 (MQ=38)
           aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCCGCGCGCAc                                                                                                                                        >  1:2506457/1‑150 (MQ=255)
           aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                        >  1:4051553/1‑150 (MQ=255)
           aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                        >  2:2769033/1‑150 (MQ=255)
           aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                        >  2:1576400/1‑150 (MQ=255)
           aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                        >  1:2429102/1‑150 (MQ=255)
            aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGTTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACg                                                                                                                                       >  2:653419/1‑150 (MQ=255)
            aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                       >  2:652861/1‑150 (MQ=255)
               agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACgggc                                                                                                                                    >  1:4072574/1‑147 (MQ=255)
               agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCAGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgg                                                                                                                                    >  1:2805426/1‑148 (MQ=255)
                 agttGAAAGGGCGTTTTCCCAAACTTGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACTAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGTAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcac                                                                                                                                  <  2:3073395/147‑1 (MQ=38)
                   aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccc                                                                                                                                >  1:3399928/1‑150 (MQ=255)
                    tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcaccca                                                                                                                               >  2:2026446/1‑149 (MQ=255)
                      aaaGGGCGTTATCCCAAACTTGTCTGGCCGGACGCCCCGTAAAATACGGCTCCTAGGGGCGGGTCTAAAATTTATTCTTATGGTTGTTTTTGGCCCGCGCCCGGTGGTTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctca                                                                                                                             <  2:3280830/150‑1 (MQ=255)
                         gggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtgacccccaccc                                                                                                                          >  1:558502/1‑150 (MQ=255)
                            cGTTATCCCAACCATGTCTGGCCGGCCGCCCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCACTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTaa                                                                                                                       <  2:2189466/150‑1 (MQ=255)
                            cGTTATCCCAAAGTTGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGTCGGAGTAAAAATTTTTCTTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTaa                                                                                                                       <  2:2189502/150‑1 (MQ=37)
                            cGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCaaa                                                                                                                       >  1:2377201/1‑150 (MQ=255)
                             gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCCaaa                                                                                                                      >  1:953701/1‑149 (MQ=255)
                             gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCAAAc                                                                                                                      >  1:338824/1‑150 (MQ=255)
                                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCCCCCCAAcccgcgc                                                                                                                >  2:907184/1‑146 (MQ=38)
                                    cAAACTTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTTTTCGTAAGTTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                               <  2:3782164/150‑1 (MQ=255)
                                    cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCCCCCGAACCCACTcc                                                                                                               >  2:2511300/1‑150 (MQ=255)
                                    cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACCCCCGAAcccacgac                                                                                                               >  2:2511256/1‑145 (MQ=255)
                                         aTGTCGGGCCGGACGACCCGGCAAATCCTGCACCGACGCGACGGACGAAAAAGTTTTCGTACGGTGGGTTTTCCCCCCCGCCGGGTGGGTGAGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                          <  2:188782/150‑1 (MQ=255)
                                          tGTCTGCCCGGACGCCCCGCCACATACGGCACCGCCGCGCGGGACGAAAATTTATTCGAACGTTTGGTTTTCTCCCTCTCCGTGTGGGCGTGGCCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCttgag                                                                                                         <  2:1105208/150‑5 (MQ=255)
                                          ggTCTGGCCGTACTACCCGGAAAATCGTGCCCGGACGCGCGGGACGAAATAGTTTTGTTATGGTTGTTTTTCGCCCCCTCCGGGTGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCttgag                                                                                                         <  2:1105099/149‑5 (MQ=255)
                                             gtGGCCGGACGACCGGGCAAATCCCGCACCGACGCGCCGGACGAAAATGTATGCGTATGGTTGGTTTGCGCCCTCTCCCTGTGGGTGCGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGga                                                                                                      <  1:3962434/149‑1 (MQ=255)
                                             cgggCCGGACGCCCCGGAAAATACTGCACCGACGCGACGGACTAAAATGTATTCGTATGGTTGGTTTGCCCCCTCTCCGGGGGGGAGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGga                                                                                                      <  1:3962516/148‑1 (MQ=255)
                                             cTGGCCGGACGACCCGGCAAATACTGCACCGAGGCGACGGACGAAAACTTTTTGGTTCGGTTGGTTTTCGCCCTCTCCCTGTGGTAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGga                                                                                                      <  2:2802809/150‑1 (MQ=37)
                                                gCCGGACGCCCCTCAAAAACCTCCACCGCCGCGCCGGCCTAAAATTTTTTCGTAGGTTTGTTTTCCCCCCACGCGGTGGGGGCGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                   <  1:585107/150‑1 (MQ=255)
                                                 ccGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGGCCCCGCCCCCGAACCCACTCCCACCAgggggggg                                                                                                  >  1:3928528/1‑145 (MQ=25)
                                                 ccGGACGACCCGGAAAATACCGACCCTACGCGACGCACGAAAATTTTTTCGATAGGTTGTTTTGCTCCCCATCCCTGTGGGAGGGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                  <  2:322189/150‑1 (MQ=255)
                                                  gggACGACCCGGCAAATCCTGCACCGACGCGCCGGACGAAAATTTTTTCGTATGGTTGGTTTCCTCCCTCTCCCTGGGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                                 <  1:3054329/149‑1 (MQ=38)
                                                  cgcacgacCGGGCCAATACTGCACCGACGGGAGGGACTTAAAAGTATTCTTACTGTGGGTTTTCTCCAACTCCGTGTGGGAGTGGGCCGGGGGGAGGGCATCAGCGCGCCCGGTCTCCCTCACCCTAACCCTCTCCCTTGAGGGAGAggg                                                                                                 <  2:1686575/147‑1 (MQ=255)
                                                  cGGCCGACCCGGCAAATACTGCCCCGACGCGAGGGACTAAAAAGTTTTCTTACGGTTGGTTTTCGCCCCCCCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                                 <  2:3254692/150‑1 (MQ=255)
                                                  cGGACGACCCGGCAAATCCTGCACCGCCGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                                 <  2:592033/150‑1 (MQ=255)
                                                  cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGGCCCCACACCCCGACCCCCACCCACCAGGGGGAggg                                                                                                 >  1:1127625/1‑150 (MQ=25)
                                                     acgacCGGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCCCCCGGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgggtta                                                                                              <  2:3715098/150‑4 (MQ=38)
                                                     acgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGGAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCCGCGCGCACGGGCACCCCCCCCCGACCCCGCGCCCCCCAGGGCGCggggcc                                                                                              >  2:4036839/1‑148 (MQ=25)
                                                        acCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCCCACCCACAACCACACCCCCAAAGGAGAGGGGGcccc                                                                                           >  2:2055698/1‑148 (MQ=25)
                                                        acCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGGAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGCACCCGCACCCACAACCGCACCCGCCAGGGGGGGGGGGCCgc                                                                                           >  2:4089293/1‑149 (MQ=25)
                                                         cccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCGCACCCGAACCCACAACCACACGGGTGGGGGGACCgct                                                                                          >  1:2626794/1‑148 (MQ=25)
                                                            ggCAATCACGGCACCGCCGGGCCGGACTAAAATGTATTCGTATGGTTGGTTTCCGCCCTCTCCGTGGGGGTGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                       <  1:2749863/150‑1 (MQ=255)
                                                            ggCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCGGCGCGCACGTTCACCCGCACCCGAAACCAATACCTCAAGGGGGTGGGGAGCggtggg                                                                                       >  2:3312620/1‑145 (MQ=25)
                                                            ggCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGGCCCCGCACCCGAACCCCCACCCACACGGGGGGGGGGGcccggcgg                                                                                       >  1:572707/1‑144 (MQ=25)
                                                            ggCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCGCACCCGAACCCACTCCCACCAGGGGGCGGGGGCCgcggga                                                                                       >  2:3312590/1‑145 (MQ=25)
                                                            ggCAAAGACTGCACCGACGCGACGGACGAACAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                       <  2:3016192/150‑1 (MQ=255)
                                                               aaaaACTCCACCGACGCGACGGACTAAAATGTATTCGTACGGTTGGTTTACCCCCTCTCCGTGTGGGTGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                                                    <  1:3000004/150‑1 (MQ=255)
                                                               aaaTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTGACCCACACCCGAACCCTCACCCACAACGGGGAGGGGGCCCAtacagcc                                                                                    >  2:4050333/1‑144 (MQ=25)
                                                                aaTACTGCACCGACGCGACGGACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGcac                                                                                   <  2:4119482/150‑1 (MQ=255)
                                                                  tACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCAACCACACCCCAACCCAAACCCACAAGGGAGGGGGGACCGCGCGAGCCCaa                                                                                 >  2:1326277/1‑150 (MQ=25)
                                                                    cTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCGAACCCACACCCTCCACGGGGAGGGGACCCCTCGAGCCCaaaa                                                                               >  2:4108338/1‑149 (MQ=25)
                                                                    cTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCAAACCCGCACCCACACGGGGGAGGGGACCCATAGAGCCCaaaa                                                                               >  2:2074834/1‑149 (MQ=25)
                                                                    cTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCGCACCCCAACCCACACCCACAACGGCGGGGGGCCCCATAGCGCACaaaa                                                                               >  1:3012772/1‑149 (MQ=21)
                                                                     tGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGGGCGCCCCCCCCCAACCCCCCAACCCCAACGGGGAGGGGGCCCCGCGAGCACcaaat                                                                              >  2:96788/1‑148 (MQ=16)
                                                                       cACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAAACCACAACCCCAAGGGAGATGGGACCCAAACAGCACAAAattt                                                                            >  2:3337472/1‑150 (MQ=25)
                                                                         ccgccgCGCCGGACTAAAACGTATTCGTACGGTTGGTTTTCTCCCTCGCCGTGTGGTTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                          <  2:3281316/150‑1 (MQ=255)
                                                                         ccGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                          <  2:3391177/150‑1 (MQ=255)
                                                                                 acggactAAAAAGTTTTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                                                  <  2:285867/150‑1 (MQ=255)
                                                                                        aaaaaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCACACCCCAAACCGCACCCGCAAGGGCGCGGGGGCCCAGCGAGCCCCAATTTTGTGGGATGtcttgttct                                                           >  2:2025801/1‑150 (MQ=16)
                                                                                         aaaaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCACACCCCAAACCACACACACACACGAGAGGGGCCCCCACCAGCACAAAAATTGTGATATAtattattaag                                                          >  1:983641/1‑147 (MQ=12)
                                                                                           attttttCGTATGTTTGTTTTTCCCCCTCGCCCGGTGGTTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                        <  2:811252/145‑1 (MQ=255)
                                                                                           aaGTATTCGTACGGTTGTTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                        <  2:953701/150‑1 (MQ=255)
                                                                                            tgttttCGTCTGGTTGTTTTTCCCCCCCGCCCGGTGGGCGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCt                                                       <  1:382673/146‑1 (MQ=255)
                                                                                             ttATTCTTAAGTTTGGTTTACCCCCTCTCCGGGTGGGAGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                      <  1:139790/149‑1 (MQ=255)
                                                                                                  gtttcggttggttTACCCCCTCCCCGTGTGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                 <  1:2021792/145‑1 (MQ=255)
                                                                                                   gTAAGGTTGTTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                <  2:3012568/150‑1 (MQ=255)
                                                                                                    tAAGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACag                                               <  1:2385231/150‑1 (MQ=255)
                                                                                                     aaGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                              <  2:1366492/150‑1 (MQ=255)
                                                                                                     aaGGTTGGTTTTCTCCCTCCCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                              <  1:2976384/150‑1 (MQ=255)
                                                                                                     aaGGTTGGTTTTCGCCCTCTCCCTGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                              <  2:1366520/150‑1 (MQ=255)
                                                                                                         tttgttTACTCCCCCTCGATGGGGGAGAGGGCGGGGTTGATGGCTTCAGCGCGCACGTTCACCCTCACCCTATCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCTCAAATTTTGAGAGATATCTTATTCTGTCTTGTAACAgaagaa                                          <  1:1352225/147‑1 (MQ=255)
                                                                                                         tgggttTGCCCCCTCGCCGTGGGGGTGCGGGGGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                          <  1:1576400/148‑1 (MQ=255)
                                                                                                          tggttTTCTCCCGCGCCCTGTGGGAGAGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                         <  1:4119438/150‑1 (MQ=255)
                                                                                                           ggttTTCTCCCGCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCGAACCCTCTCCCTCAAGGGAGAGGGGACCGAACGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTAGTAACAGAAGAACa                                        <  2:297121/150‑1 (MQ=255)
                                                                                                              tttCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCGCCCCCGAACCCACGCCCACAATGGGGAGGGGCCCGATCCAGCCCAAATTTTGAGAGATATCTTGTTCTGTCTTGTAACAGAAGAACAGaa                                     >  2:2059319/1‑150 (MQ=25)
                                                                                                                    cctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCAAAACCTCACCCTCAAGGGGGTGGGGGCCGATCGAGCACAAATTTTGAGGGATATATTCTTCTGTCTTGTGAAAGAAGTAAAGAAAATCgg                               >  2:3094714/1‑150 (MQ=37)
                                                                                                                         ccctGGGGGAGCGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                          <  1:3307000/150‑1 (MQ=255)
                                                                                                                          cctGGGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCCCGTGCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                         <  1:1562785/150‑1 (MQ=255)
                                                                                                                             gtgGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcg                      <  1:1557801/150‑1 (MQ=255)
                                                                                                                              tgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCAAAACCACAACCACCAGGGGGAGGGGGCCCATCGAGCCCAAAATTTGAGAGATAAAATATTATGTATTGTAACAGAAAAACATCAAAAAGGGCGTTTtcgc                     >  2:1376930/1‑150 (MQ=255)
                                                                                                                              ggggtgcgggCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                     <  1:1437939/143‑1 (MQ=255)
                                                                                                                              ggggggtgggCGGGGGGGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                     <  2:377491/143‑1 (MQ=255)
                                                                                                                              ggggAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                     <  1:3884677/149‑1 (MQ=255)
                                                                                                                                   gagGGCCGGGGTGAGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGaa                <  1:3051556/150‑1 (MQ=255)
                                                                                                                                        ccGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCGAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg           <  1:1325402/150‑1 (MQ=255)
                                                                                                                                         cGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc          <  1:433476/150‑1 (MQ=255)
                                                                                                                                         cGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc          <  1:4108338/150‑1 (MQ=255)
                                                                                                                                          ggggTGAGGGCATCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg         <  1:2026446/150‑1 (MQ=255)
                                                                                                                                          ggggTGAAGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg         <  2:610582/150‑1 (MQ=255)
                                                                                                                                            ggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga       <  1:2711821/150‑1 (MQ=255)
                                                                                                                                                 ggccATCAGCGCCCACGTTCACCCTCACCCTACCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAATTCGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa  <  1:458293/147‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAA  >  NC_000913/164432‑164726

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: