Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,581 A→C 20.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,5810AC20.5% 185.9 / 19.0 112intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (34/53);  new base C (4/21);  total (38/74)
Fisher's exact test for biased strand distribution p-value = 3.37e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 113 total aligned reads displayed.

CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGC  >  NC_000913/164432‑164729
                                                                                                                                                     |                                                                                                                                                    
ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                      >  1:2593872/1‑150 (MQ=255)
  gTGAGGTGCAAAAAGAGATGAAAGGGCTTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGTTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATc                                                                                                                                                    <  2:71963/150‑1 (MQ=255)
  gTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATc                                                                                                                                                    <  2:2283210/150‑1 (MQ=255)
   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                   >  2:1325102/1‑150 (MQ=255)
        tGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGGGGGCCGGGGTGAGGGCATCAgcgcg                                                                                                                                              >  2:3003315/1‑150 (MQ=255)
         gCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGGAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCgcgcgc                                                                                                                                             >  1:745790/1‑150 (MQ=255)
          cAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCa                                                                                                                                            >  2:506899/1‑150 (MQ=255)
          cAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTAGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcgcc                                                                                                                                            >  2:1931315/1‑149 (MQ=255)
            aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                          >  1:3400426/1‑150 (MQ=255)
                   aTGAAAGGGCGTTATCCCAAACCTGTCTGGCCGGACGACCCGGCAAATCCGGCACCGCCGCGACGGACGACAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGGGAGGGCCGCCGCGCGCACgggggccc                                                                                                                                   >  2:3664517/1‑143 (MQ=37)
                   aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcgccc                                                                                                                                   >  1:699288/1‑150 (MQ=255)
                   aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACgggggccc                                                                                                                                   >  1:2352665/1‑143 (MQ=255)
                   aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcaccc                                                                                                                                   >  2:3974175/1‑150 (MQ=255)
                     gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTggaccccc                                                                                                                                 >  1:34187/1‑148 (MQ=255)
                     gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTggacccac                                                                                                                                 >  1:669919/1‑148 (MQ=255)
                     gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgccccccc                                                                                                                                 >  1:4344904/1‑148 (MQ=255)
                     gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcacccac                                                                                                                                 >  1:4344553/1‑148 (MQ=255)
                             gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCACCCCAAc                                                                                                                         >  2:3729248/1‑150 (MQ=255)
                             gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACCCCCCAAc                                                                                                                         >  2:1265371/1‑150 (MQ=38)
                                  cccAAACCTGTCTGCGGGTACGACCCGGAAAATACTGCACGGCCGGGACGCACTAAAATGTATTCGTAAGGTTGGTTTTCTCCCCCTCCCTGTGGGCGTGGGCCGGGTTGGGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctct                                                                                                                    <  2:2745467/150‑1 (MQ=255)
                                   ccAAACTTGTCTGGCCGGACGACCCGGAAAATACTGCCCCGACGCGACGGCCGAAAAATTTTTGGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                                   <  2:1402197/150‑1 (MQ=37)
                                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                                   <  2:3604364/150‑1 (MQ=255)
                                   ccAAACATGACCGCCCGGACGACCCGCCAAAAACTCCACCGACGCGACGCACCAAAAAGTATTCGAAAGGTTGGTTTTCACCCACTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                                   <  2:3503004/150‑1 (MQ=37)
                                    cAAACTTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTTTTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                                  <  2:684667/150‑1 (MQ=255)
                                         gtGTCGGGCGGGACGACCGGGCTAAGACGGCCCCGCGGGGACGGAGGAAAAAGTTTTCGTACGGTTGGTTTGCTCCCTGTCCCTGGGGGAGCGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTACCCCTCTCCCTCaa                                                                                                             <  2:2376622/149‑1 (MQ=255)
                                         ctGTCGGGCCGGACGCCCCGCCAAATCCTCCACCCCCGCGCCGGCCCAAAAAGTTTTCGTACGGTGGGTTTCCCCCCCCACCCTGTGGGAGCGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                             <  2:2376637/149‑1 (MQ=255)
                                          ggccgggCCGGACGCCCCGCCAACCCCGCCACCGCCGCGCCGTCCGAAAAAGTATTCGTACGGTTGGTTTACTCCCGCGCCGTGTGGCCGAGGGCCGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                            <  2:825186/145‑1 (MQ=255)
                                             cgggCCGGACGACCCGGAAAATACTGCACCGCCGCGACGTACTAAAAAGTATTCGTATGGTTGTTTTACTCCCTCTCCCGGTGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGga                                                                                                         <  1:1772575/148‑1 (MQ=255)
                                                gCCGCACGACCCGCCAAATCCTGCACCGCCGCGACGGACGAAAATGTATTGGTTCGGTTGGTTTGCTCCCTCCCCCTGTGGGTGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                      <  2:3784362/150‑1 (MQ=37)
                                                 ccGGCCGACCCGGCACCCCCTGCCCCGCCGCGCCGGCCGCAAACGTAGTCGTCCGGTTGGTTTTCGCCCTCTCCCGGGGGGAGCGGGCCGGGGGGCGGGCATCAGCGCGCCCGGTCACCCTCACCCTAACCCTCCCCCTCAAGGGAGAgg                                                                                                     <  2:2181704/150‑1 (MQ=255)
                                                 ccGGACGCCCCGGAACATACTCAACCGACGCGCGGGACTAAAATTTATTCGTATGTTTGTTTTTCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgg                                                                                                     <  1:1743424/150‑1 (MQ=255)
                                                     acgacCCGGCAAATCGTGCCCCGACGCGACGGACGAAAAAGTATTGGTAGGTTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTACCCCTCTCCCTCAAGGGAGAGGGGAc                                                                                                 <  2:1387975/150‑1 (MQ=37)
                                                       tacCCGGCAAATACTGCACCGACGCGACGGACGAAAATGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCAACCGAAACCGAAAACACTAACACAAGGGAGAgggggcgg                                                                                               >  1:1043967/2‑146 (MQ=21)
                                                       gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGTGGAACCCCACCCAAAAACACACCCACAAGCGGGTGGGGGCCg                                                                                               >  2:1382402/1‑150 (MQ=25)
                                                          ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGGGCACGGTGAGCCGCGCCCTCACCCGCTACCACAAGGGGGGGGGGGCCgcgg                                                                                            >  1:3019887/1‑147 (MQ=25)
                                                          ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGCCCCACACCCGAAACCACTACCAACAGGGGGAGGGGGCCGagc                                                                                            >  1:670004/1‑148 (MQ=25)
                                                          ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGCCCCCCCCCCCAAACCCGCTCCCCCCAGGGGGGGGGGCCCGagc                                                                                            >  2:3901284/1‑148 (MQ=25)
                                                           cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGCACCCGCACCCCAACCCACACCCCCCACGGCGCGGGGGCCGAGCg                                                                                           >  2:4106148/1‑150 (MQ=25)
                                                           cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAGGTGCACCCACCCCCAAACCCACACCCTCAGGGGCGCGGGGCCCgcgcg                                                                                           >  2:3555320/1‑146 (MQ=25)
                                                           cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCTCAACCAAACCCACCCCCGCCAGGGGGGGGGGGCCCATCg                                                                                           >  1:285333/1‑150 (MQ=25)
                                                            ggaaaaTCCGGCACCGACGCGACGGACGAAAATGTTTTGGTACGGTTGGTTTTCGCCCACACCCGGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                          <  2:454427/147‑1 (MQ=255)
                                                            ggCAAATACGGCCCCGCCGCGACGGACGAAAACGTATTGGTATGGTTGGTTTTCCCCCTCCCCCTGGGGGAGCGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                          <  2:1698865/150‑1 (MQ=255)
                                                             gCAAATACGGCACCGACGCGACGGACGACAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                         <  2:3259764/150‑1 (MQ=255)
                                                              cAAATACGGCACCGCCGCGACGGACGAAAAAGTATGCGTAAGGTTGGTTTTCTCCCTCGCCCTGTGGTCGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGc                                                                                        <  2:2488079/150‑1 (MQ=255)
                                                                   aCTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCCCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACaaa                                                                                   >  2:1898816/1‑150 (MQ=255)
                                                                   aCTGCACCGACGCGACGGACGAAAAAGTATGCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCCCCCCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCCCaaa                                                                                   >  2:402767/1‑150 (MQ=255)
                                                                     tGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCCGCGCGCACGTGCCACCGCACCCGAACCCACTCCCGCAAGGGGGGGGGGACCGCTAGAGCCCaaaat                                                                                 >  2:2891126/1‑148 (MQ=25)
                                                                     agCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCGCACCCCAACCCACACCCACAAAGGCGACGGGCCCGCTCCAGcaccacat                                                                                 >  2:2335189/2‑145 (MQ=25)
                                                                       ccccGACGCGACGGACGAAAAATTATTCGTACGGTTGGTTTTCTCCCACTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                               <  1:3586140/148‑1 (MQ=255)
                                                                       cACCGCCGCGACGGACGAAAATGTATTCGTATGGTTGGTTTTCTCCCTCCCCCTGTGGGTGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAATtttt                                                                               <  2:2310358/150‑1 (MQ=255)
                                                                       cACCGACGCGAGGGAGTAAAATTTTTTGTTATGGTTGTTTTTCTCCCTGTCCGGGTGGTTGTGGGCGGGGTTTTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                               <  2:1238883/150‑1 (MQ=255)
                                                                       cACCGACGCGAGGGACTAAAATGTTTTCTTATGGTTGGTTTTCTCCCTCCCCGTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                               <  1:1361873/150‑1 (MQ=255)
                                                                         ccGACGCGCCGGACTAAAATTTTTTCTTATGTTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                             <  1:2734164/150‑1 (MQ=37)
                                                                              gcgacgGACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat                                                                        <  2:299897/150‑1 (MQ=255)
                                                                              gcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCACAAACAAAACCACAACCGCAAGGGAGAGGGGGCCCCAAGAGCACAAAATATGAGTGat                                                                        >  2:493242/1‑150 (MQ=16)
                                                                                    gacgAAAATGTATTCGTACGGTTGGTTTCCCCCCCCGCCCTGGGGGTGCGGGGCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttc                                                                  <  1:3342574/150‑1 (MQ=255)
                                                                                           tttttttCGTTGGGTGGGTTTCCCCCCTGCCGGTGTGGGGGTGGGCGGGGGTGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                           <  1:144562/145‑1 (MQ=255)
                                                                                           aaGTATCCTTATGGTTGGTTTTCTCCCCCCCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                           <  1:672957/150‑1 (MQ=38)
                                                                                            gttttacGTAGGGTTGGTTGTCCCCAACTCCCTGGGGGTGTGGGCGGGGTTGAGTGCATCAGCGCGCACGTTCACCCTCACCCTTACCCTCTCCCTCAAGTGAGATGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCt                                                          <  1:4002091/144‑1 (MQ=255)
                                                                                            agtagtCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCt                                                          <  2:700314/150‑1 (MQ=255)
                                                                                            aGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCt                                                          <  2:1082003/150‑1 (MQ=255)
                                                                                             ttATTCGTATGGTTGTTTTTCTCCCGCTCCGGGTGGGTGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                         <  1:1645209/149‑1 (MQ=255)
                                                                                             gTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                         <  2:1384084/150‑1 (MQ=255)
                                                                                             gTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                         <  2:1384526/150‑1 (MQ=255)
                                                                                             gTATTCGTAAGGTTGGTTTTCTCCCTCTCCCGGTGGGAGAGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                         <  2:2194751/150‑1 (MQ=255)
                                                                                             gTATACGAACGGTGGGTTTCCCCCCACGCCCGGTGGTCGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                         <  1:849101/150‑1 (MQ=255)
                                                                                                 tcgtctggttggttTTCTCCCTCTCCGTGTGGGTGTGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                     <  1:339820/144‑1 (MQ=255)
                                                                                                  cGTATGGTTGTTTTACTCCCTCTCCCTGTGGGTGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                    <  1:3739507/150‑1 (MQ=39)
                                                                                                    tatgtttgtttTACTCCCCCCCCGGGTGGGAGCGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACag                                                  <  2:1624059/145‑1 (MQ=255)
                                                                                                      aGGTTGGTTTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                <  2:3598086/150‑1 (MQ=255)
                                                                                                         ttggttTTCGCCCACTCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                             <  2:3460705/150‑1 (MQ=255)
                                                                                                         ttggttTTCCCCCTCTCCGTGGGGGAGTGGGCGGGGGTGGGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                             <  2:4026770/150‑1 (MQ=255)
                                                                                                          tggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                            <  2:2421969/150‑1 (MQ=255)
                                                                                                           ggttTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                           <  1:4283046/150‑1 (MQ=255)
                                                                                                           ggttTTCTCCCTCACCCGGGGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                           <  1:4283030/150‑1 (MQ=255)
                                                                                                            gttTTCTCCCTCGCCGGGTGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAg                                          <  1:1719084/150‑1 (MQ=255)
                                                                                                            gttTTCACCCTCTCCCTGTGGGTGTGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAg                                          <  1:841908/150‑1 (MQ=255)
                                                                                                                  tccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCGCACCCGAACCCACGCCCCCAAGGGCGCGGGGCCCGCGCGCGCCCCCATGTTGTGAGATATCTTCGTCTGTCGTGTAACAGACGGACAGaaaaac                                    >  2:3900462/1‑148 (MQ=16)
                                                                                                                    cctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCCCCCGCACCCAAACCCGCGCCCGCCAGGGCGGGGGGCCCGATCGGGCACAAAATTTGAGAGCTATATTAATCGGTCGTGGAAACGAAGAACACAAAAACgg                                  >  2:2882496/1‑150 (MQ=16)
                                                                                                                     ctctcccgGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTTACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                                 <  2:304268/150‑1 (MQ=255)
                                                                                                                              tgGGAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                        <  2:2593872/150‑1 (MQ=255)
                                                                                                                              ggggAGCGGGCCGGGGGGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                        <  1:1325102/149‑1 (MQ=255)
                                                                                                                                ggAGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                      <  1:479139/150‑1 (MQ=255)
                                                                                                                                ggAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                      <  2:3479750/150‑1 (MQ=255)
                                                                                                                                  agagGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCCCACCCAAAACCACACCCACAAGGGCGAGGGGCCCCAGCGAGCACAAATTTTTGGAGATAAAATATTATGTCTTGTAAAAGAAGAAAAGAAAAAACGGCCTTTTCGccagc                    >  2:665373/1‑147 (MQ=255)
                                                                                                                                       gcgggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATAtt               <  2:1717567/147‑1 (MQ=255)
                                                                                                                                       gCCGGGGTGAGGGCATCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGACAATCGGGCTTTTGCGCCTGAATAtt               <  2:2763974/150‑1 (MQ=255)
                                                                                                                                          ggggTGCGGGACTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGTTCGAGCACTAATATTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGTAAAGAAAATCGTGCTTTTGGGCCTTAATATTGCg            <  1:2587058/150‑1 (MQ=38)
                                                                                                                                          ggggTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCACTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg            <  2:847500/150‑1 (MQ=255)
                                                                                                                                          ggggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg            <  2:544449/150‑1 (MQ=255)
                                                                                                                                          ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAACTTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg            <  2:3221886/150‑1 (MQ=255)
                                                                                                                                          ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCNGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg            <  2:2780016/150‑1 (MQ=255)
                                                                                                                                            ggTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga          <  1:497657/150‑1 (MQ=255)
                                                                                                                                            ggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga          <  1:497656/150‑1 (MQ=255)
                                                                                                                                            ggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga          <  1:2061802/150‑1 (MQ=255)
                                                                                                                                               gcgggCCTCAGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGCTCGAGCACAAATTTGGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaa       <  1:3098685/148‑1 (MQ=255)
                                                                                                                                               gAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaa       <  2:3790083/150‑1 (MQ=255)
                                                                                                                                               gAGGGAATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaa       <  2:3014715/150‑1 (MQ=255)
                                                                                                                                                 gggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa     <  1:346152/150‑1 (MQ=255)
                                                                                                                                                  ggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa    <  2:597062/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGC  >  NC_000913/164432‑164729

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: