Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,139,632 A→G 8.7% intergenic (+39/+84) yghU → / ← hybG disulfide reductase/organic hydroperoxide reductase/hydrogenase 2 accessory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,139,6320AG 8.7% 435.1 / 22.7 151intergenic (+39/+84)yghU/hybGdisulfide reductase/organic hydroperoxide reductase/hydrogenase 2 accessory protein
Reads supporting (aligned to +/- strand):  ref base A (77/60);  new base G (8/5);  total (85/66)
Fisher's exact test for biased strand distribution p-value = 7.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

Only 100 of 152 total aligned reads displayed.

AACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTT  >  NC_000913/3139484‑3139780
                                                                                                                                                    |                                                                                                                                                    
aaCGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATaa                                                                                                                                                     >  2:79330/1‑150 (MQ=255)
    ggggAGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCg                                                                                                                                                 <  2:2465651/149‑1 (MQ=255)
        cGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCa                                                                                                                                             >  2:1137273/1‑150 (MQ=255)
             tGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCAcgccg                                                                                                                                        >  2:2197999/1‑150 (MQ=255)
              gTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCAcgccgc                                                                                                                                       >  1:2190668/1‑150 (MQ=255)
              gTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCAcgccgc                                                                                                                                       >  2:2393927/1‑150 (MQ=255)
              gTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCAcgccgc                                                                                                                                       >  1:1634000/1‑150 (MQ=255)
                  aCCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                                                   <  2:455122/150‑1 (MQ=255)
                   ccGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCAATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGCTTGGTGTTAGTCGAAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                                                  >  2:1814766/1‑150 (MQ=255)
                   ccGCACCAACGGACCGATGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                                                  >  1:750927/1‑150 (MQ=255)
                    cGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                 >  2:1514338/1‑150 (MQ=255)
                    cGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTATTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                 <  1:678958/150‑1 (MQ=255)
                    cGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGGCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                 <  1:2363529/150‑1 (MQ=255)
                     gCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                <  1:344081/150‑1 (MQ=255)
                      cACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                                                                               >  2:1287035/1‑150 (MQ=255)
                         cAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATGCGGCGTGAACGCCTTATCCGGcctgg                                                                                                                            <  2:1213076/150‑5 (MQ=255)
                             ggACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGc                                                                                                                        >  2:1946451/1‑150 (MQ=255)
                               aCCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATTCGGAATATAAGCGTCAGGGGTAAGGGTTGTTGTTCGTCGAAGCAAGCCATCCAGGACGGATAATGCGTTCAAGCCGCATCCGGCAATGGTGcac                                                                                                                      <  2:974807/150‑1 (MQ=255)
                                  gCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAAt                                                                                                                   >  2:565431/1‑150 (MQ=255)
                                     gAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCTTTCACGACGCATCCGGCAATCGTGCACAATGcc                                                                                                                <  2:1449904/150‑1 (MQ=255)
                                     gAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGcc                                                                                                                <  1:262/150‑1 (MQ=255)
                                      aaTGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCt                                                                                                               <  1:463507/150‑1 (MQ=255)
                                                 gCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCt                                                                                                    >  1:426425/1‑150 (MQ=255)
                                                 gCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCt                                                                                                    >  2:2394725/1‑150 (MQ=255)
                                                  cATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTg                                                                                                   >  2:1271461/1‑150 (MQ=255)
                                                  cATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTg                                                                                                   >  1:490461/1‑150 (MQ=255)
                                                  cATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTg                                                                                                   >  1:1110165/1‑150 (MQ=255)
                                                   aTGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGc                                                                                                  >  2:872676/1‑150 (MQ=255)
                                                    tGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGcc                                                                                                 >  1:2360983/1‑150 (MQ=255)
                                                     gAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGCTAAGCCTCAGGGGGAAGGGTTGGGGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCcg                                                                                                <  2:2103849/150‑1 (MQ=255)
                                                        cgcCATGACGCCAGTGATTTCGAGACGAATACGGAAGATCAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGt                                                                                             <  2:2363346/150‑1 (MQ=255)
                                                        cgcCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGt                                                                                             >  1:678584/1‑150 (MQ=255)
                                                            aTGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTa                                                                                         <  2:1687648/150‑1 (MQ=255)
                                                                cGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAg                                                                                     <  2:1113549/150‑1 (MQ=255)
                                                                cGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAg                                                                                     <  1:746427/150‑1 (MQ=255)
                                                                cGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAg                                                                                     >  2:2114855/1‑150 (MQ=255)
                                                                 gCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAgg                                                                                    <  1:1896584/150‑1 (MQ=255)
                                                                 gCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAgg                                                                                    >  1:87121/1‑150 (MQ=255)
                                                                 gCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTGGGGGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATAAGACGCGGCAGCGTCGCATCAgg                                                                                    <  1:1154515/150‑1 (MQ=38)
                                                                  ccAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGTATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATAAGACGCGGCAGCGTCGCATCAGGc                                                                                   <  2:2222394/150‑1 (MQ=38)
                                                                  ccAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGc                                                                                   >  2:856451/1‑150 (MQ=255)
                                                                  ccAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGc                                                                                   >  2:2155868/1‑150 (MQ=255)
                                                                    aGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGGGCTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATAAGACGCGGCAGCGTCGCATCAGGcat                                                                                 <  1:538896/150‑3 (MQ=37)
                                                                     gTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTa                                                                                >  1:415146/1‑150 (MQ=255)
                                                                          ttCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGaaa                                                                           <  2:965680/150‑1 (MQ=255)
                                                                           tCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAg                                                                          <  1:813808/150‑1 (MQ=255)
                                                                            cGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGc                                                                         >  1:785088/1‑150 (MQ=255)
                                                                             gagaCGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCt                                                                        <  2:564299/150‑1 (MQ=255)
                                                                             gagaCGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCt                                                                        <  2:2360983/150‑1 (MQ=255)
                                                                             gagaCGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCt                                                                        >  2:2480003/1‑150 (MQ=255)
                                                                             gagaCGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCt                                                                        >  1:2401378/1‑150 (MQ=255)
                                                                                aCGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAat                                                                     <  1:462278/150‑1 (MQ=255)
                                                                                aCGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAat                                                                     <  2:1826942/150‑1 (MQ=255)
                                                                                 cGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAatc                                                                    <  1:854790/150‑1 (MQ=255)
                                                                                     tACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAAtcatca                                                                <  1:1514338/150‑1 (MQ=255)
                                                                                             aTAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTgg                                                        <  1:856451/150‑1 (MQ=255)
                                                                                              tAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGt                                                       >  1:1839490/1‑150 (MQ=255)
                                                                                                   gTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGt                                                  >  2:2014523/1‑150 (MQ=255)
                                                                                                     cAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCg                                                >  1:233168/1‑150 (MQ=255)
                                                                                                     cAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCg                                                >  2:2167228/1‑150 (MQ=255)
                                                                                                         ggTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAAt                                            >  1:2224236/1‑150 (MQ=255)
                                                                                                         ggTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAAt                                            >  1:1174160/1‑150 (MQ=255)
                                                                                                          gTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATc                                           >  2:904394/1‑150 (MQ=255)
                                                                                                             aGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCAt                                        <  2:2037179/150‑1 (MQ=255)
                                                                                                              gggTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCAtt                                       >  2:1714214/1‑150 (MQ=255)
                                                                                                              gggTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCAtt                                       <  1:1271461/150‑1 (MQ=255)
                                                                                                              gggTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCAtt                                       >  2:495352/1‑150 (MQ=255)
                                                                                                                   ggTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGgcg                                  >  1:2296833/1‑150 (MQ=255)
                                                                                                                   ggTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGACGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGgcg                                  >  2:704606/1‑150 (MQ=255)
                                                                                                                      gtTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAg                               >  1:351168/1‑150 (MQ=255)
                                                                                                                          gtcgCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCg                           <  1:2314910/150‑1 (MQ=255)
                                                                                                                            cgCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTc                         >  2:323218/1‑150 (MQ=255)
                                                                                                                            cgCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTc                         <  1:2175613/150‑1 (MQ=255)
                                                                                                                            cgCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTc                         >  2:2248519/1‑150 (MQ=255)
                                                                                                                            cgCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTc                         >  2:1331438/1‑150 (MQ=255)
                                                                                                                            cgCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTc                         <  1:565431/150‑1 (MQ=255)
                                                                                                                             gcagcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCt                        >  1:1957167/1‑150 (MQ=255)
                                                                                                                                gcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAAt                     <  2:731126/150‑1 (MQ=255)
                                                                                                                                  aaGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAAtgt                   <  1:2393927/150‑1 (MQ=255)
                                                                                                                                   aGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAAtgtg                  <  1:1692138/150‑1 (MQ=255)
                                                                                                                                     ccATGCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtggc                >  1:1130596/1‑150 (MQ=255)
                                                                                                                                     ccATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtggc                >  1:1699666/1‑150 (MQ=255)
                                                                                                                                      cATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtggct               >  2:1643248/1‑150 (MQ=255)
                                                                                                                                       acgattgCCGGATGCGGCGTGAACGCCTTATCCGGCAATCGTGCACAATGCCTGATGCGATGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtggctt              >  2:2280739/7‑150 (MQ=38)
                                                                                                                                        cgattgCCGGATGCGGCGTGAACGCCTTATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtggcttt             >  1:1128757/6‑150 (MQ=38)
                                                                                                                                         ccAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtggctttg            <  1:1004575/150‑1 (MQ=255)
                                                                                                                                          attgCCGGATGCGGCGTGAACGCCTTATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtggctttgg           >  1:975711/4‑150 (MQ=38)
                                                                                                                                           aGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtggctttggc          >  1:1805988/1‑150 (MQ=255)
                                                                                                                                            tgCCGGATGCGGCGTGGACGGCGTAGCCGGCAATCGTGCACAAGGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAACGCTAATCATCATGCGCTGGGAAGGTCGGAGTCGCTTTGGCGCAGTGCGTCTAATGtggctttggct         >  2:2115025/2‑150 (MQ=255)
                                                                                                                                             gCCGGATGCGGCGTGAACGCCTTATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtggctttggctt        >  2:113987/1‑150 (MQ=38)
                                                                                                                                             gCCGGATGCGGCGTGAACGCCTTATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtggctttggct         <  1:485233/149‑1 (MQ=39)
                                                                                                                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGtggctttggctt        >  1:839679/1‑150 (MQ=255)
                                                                                                                                               cGGATAAGGCTTTCACGCCGCATCCGGCAATCTTTCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACTCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCa      <  2:495289/150‑1 (MQ=255)
                                                                                                                                               cGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCa      <  1:909004/150‑1 (MQ=255)
                                                                                                                                               cGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCa      <  1:2248519/150‑1 (MQ=255)
                                                                                                                                               cGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCa      <  1:323218/150‑1 (MQ=255)
                                                                                                                                                ggATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCAt     >  1:122688/1‑150 (MQ=255)
                                                                                                                                                 gATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATc    <  2:785088/150‑1 (MQ=255)
                                                                                                                                                   tAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCtt  <  2:2224236/150‑1 (MQ=255)
                                                                                                                                                   tAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCtt  <  2:159699/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGCAGTGCGTCTAATGTGGCTTTGGCTTCATCTT  >  NC_000913/3139484‑3139780

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: