Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,787,619 Δ1 bp 5.4% coding (782/1260 nt) envC → murein hydrolase activator EnvC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,787,6190T. 5.4% 501.4 / 21.9 166coding (782/1260 nt)envCmurein hydrolase activator EnvC
Reads supporting (aligned to +/- strand):  ref base T (81/76);  new base . (4/5);  total (85/81)
Fisher's exact test for biased strand distribution p-value = 7.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-01

Only 100 of 169 total aligned reads displayed.

CTGGCAGGGCTGGAGTCTTCCATCCAGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACGCTGCATCGCTATG  >  NC_000913/3787471‑3787768
                                                                                                                                                    |                                                                                                                                                      
cTGGCAGGGCTGGAGTCTTCCATCCAGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGtt                                                                                                                                                       >  2:1443046/1‑150 (MQ=255)
  ggCAGGGCTGGAGTCTTCCATCCAGCAAGGTCAGCAACATTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTcg                                                                                                                                                     >  2:2101352/1‑150 (MQ=255)
  ggCAGGGCTGGAGTCTTCCATCCAGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTcg                                                                                                                                                     <  2:1859469/150‑1 (MQ=255)
    cAGGGCTGGAGTCTTCCATCCAGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTcgcg                                                                                                                                                   <  2:1085036/150‑1 (MQ=255)
    cAGGGCTGGAGTCTTCCATCCAGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTcgcg                                                                                                                                                   >  1:648174/1‑150 (MQ=255)
     aGGGCTGGAGTCTTCCATCCAGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGa                                                                                                                                                 <  1:337689/150‑1 (MQ=255)
       ggCTGGAGTCTTCCATCCAGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGAcc                                                                                                                                               >  2:889384/1‑150 (MQ=255)
       ggCTGGAGTCTTCCATCCAGCAAGCTCAGCAACACTTGAGGGAGCTGCGCGCCAACGAATCCGGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGTcc                                                                                                                                               <  2:1399440/150‑1 (MQ=255)
        gCTGGAGTCTTCCATCCAGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCg                                                                                                                                              >  2:2201034/1‑150 (MQ=255)
        gCTGGAGTCTTCCATCCAGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCg                                                                                                                                              >  2:1193370/1‑150 (MQ=255)
         cTGGAGTCTTCCATCCAGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGc                                                                                                                                             >  2:1020163/1‑150 (MQ=255)
          tGGAGTCTTCCATCCAGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGcc                                                                                                                                            <  2:766221/150‑1 (MQ=255)
           ggAGTCTTCCATCCAGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCa                                                                                                                                           >  2:235264/1‑150 (MQ=255)
            gAGTCTTCCATCCAGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTCGCGAACCGcctg                                                                                                                                          <  1:1968009/150‑3 (MQ=255)
              gTCTTCCATCCAGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCagaa                                                                                                                                        >  1:1735302/1‑150 (MQ=255)
              gTCTTCCATCCAGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTCGCGAACCGcctggg                                                                                                                                        <  1:71183/150‑5 (MQ=255)
                cTTCCATCCAGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCagaaag                                                                                                                                      >  1:1449371/1‑150 (MQ=255)
                 ttCCATCCAGCAAGGTCAGCAACAGTTGAGCGAGCTGAGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCGAGGCGGTTCGCGACCTCCagaaagt                                                                                                                                     <  2:1703301/150‑2 (MQ=255)
                    catccaGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGc                                                                                                                                  >  2:995180/1‑150 (MQ=255)
                    catccaGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGc                                                                                                                                  >  2:517112/1‑150 (MQ=255)
                    catccaGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGc                                                                                                                                  >  1:1616541/1‑150 (MQ=255)
                     atccaGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCg                                                                                                                                 <  2:1830214/150‑1 (MQ=255)
                     atccaGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCg                                                                                                                                 <  1:447529/150‑1 (MQ=255)
                        caGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGAcg                                                                                                                              >  2:263216/1‑150 (MQ=255)
                              ggTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCaaa                                                                                                                        >  1:2104191/1‑150 (MQ=255)
                               ctCACCAACAGTAGAGACAGATTCTCGCCAACGAATAACGTCTGCGAAACACCATTGAACGTTAGAAAGCATAGGCCAAAGCGCGTCCAGAACGAGAAGCACGTGAGGCCCAGGCGCTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAg                                                                                                                       >  2:866856/2‑150 (MQ=255)
                                  aGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGCCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGca                                                                                                                    <  1:79448/150‑1 (MQ=255)
                                   gCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGcac                                                                                                                   >  2:1902355/1‑150 (MQ=255)
                                         gTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGGGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTa                                                                                                             <  1:1976065/150‑1 (MQ=255)
                                           tGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCATGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACa                                                                                                           <  2:1883995/150‑1 (MQ=255)
                                               cgagcTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAacc                                                                                                       <  2:764407/150‑1 (MQ=255)
                                                  gcTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAaccgac                                                                                                    >  1:1320537/1‑150 (MQ=255)
                                                  gcTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAaccgac                                                                                                    <  2:543926/150‑1 (MQ=255)
                                                   cTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAaccgacc                                                                                                   <  1:2201034/150‑1 (MQ=255)
                                                    tGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAaccgaccg                                                                                                  >  1:872540/1‑150 (MQ=255)
                                                        cgcCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGACACCACCTACAAACCGACCGAAAg                                                                                              >  2:348468/1‑150 (MQ=255)
                                                         gcCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGc                                                                                             <  1:104111/150‑1 (MQ=255)
                                                            aaCGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGATGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGaa                                                                                          <  1:1784748/150‑1 (MQ=255)
                                                                aaTCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAAt                                                                                      >  2:1752178/1‑150 (MQ=255)
                                                                aaTCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAAt                                                                                      <  2:648174/150‑1 (MQ=255)
                                                                   cccGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGc                                                                                   <  2:570465/150‑1 (MQ=255)
                                                                    ccGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCt                                                                                  <  2:2197573/150‑1 (MQ=255)
                                                                    ccGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCt                                                                                  <  2:1611751/150‑1 (MQ=255)
                                                                    ccGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCt                                                                                  >  2:303117/1‑150 (MQ=255)
                                                                     cGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTg                                                                                 <  1:1666686/150‑1 (MQ=255)
                                                                     cGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGGGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTg                                                                                 <  2:1585387/150‑1 (MQ=255)
                                                                         tGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGt                                                                             >  2:1763133/1‑150 (MQ=255)
                                                                          gCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTc                                                                            >  1:1066649/1‑150 (MQ=255)
                                                                           cGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTcc                                                                           <  2:2041056/150‑1 (MQ=255)
                                                                           cGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTcc                                                                           <  2:1066649/150‑1 (MQ=255)
                                                                            gTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTccc                                                                          <  2:954371/150‑1 (MQ=255)
                                                                            gTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTccc                                                                          >  2:1225359/1‑150 (MQ=255)
                                                                                cAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTAc                                                                      >  2:46161/1‑150 (MQ=255)
                                                                                cAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTAc                                                                      >  1:1009188/1‑150 (MQ=255)
                                                                                    aTTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCggt                                                                  <  2:1403436/150‑1 (MQ=255)
                                                                                      tGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCggtgg                                                                <  1:178089/150‑1 (MQ=255)
                                                                                       gCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGCAGCGACGCGCAAAGGCCCCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTCCCggtggg                                                               >  1:448304/1‑149 (MQ=255)
                                                                                          cGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTg                                                            >  1:1317640/1‑150 (MQ=255)
                                                                                           gTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTgg                                                           >  1:1306196/1‑150 (MQ=255)
                                                                                           gTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTgg                                                           <  1:547482/150‑1 (MQ=255)
                                                                                               ggAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGgcgc                                                       >  2:2114622/1‑150 (MQ=255)
                                                                                                  aGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGcgcc                                                    >  1:497216/1‑150 (MQ=255)
                                                                                                  aGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGcgcc                                                    >  1:1075895/1‑150 (MQ=255)
                                                                                                   gCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGcgccg                                                   >  1:619756/1‑150 (MQ=255)
                                                                                                     cgcgGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCCCCACCTACAAACCGACCGAACGCGAAAACTCGCTGATGTCGCGTACCGGTGGTCTGGGCGCGCcgcg                                                 >  2:1294680/1‑150 (MQ=255)
                                                                                                     cgcgGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCcgcg                                                 <  1:1193370/150‑1 (MQ=255)
                                                                                                     cgcgGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCcgcg                                                 >  1:2460490/1‑150 (MQ=255)
                                                                                                     cgcgGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCcgcg                                                 >  2:60414/1‑150 (MQ=255)
                                                                                                        ggcgAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTgg                                              <  2:1614061/150‑1 (MQ=255)
                                                                                                           gAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCa                                           <  1:995180/150‑1 (MQ=255)
                                                                                                             aaGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAgg                                         <  1:2027969/150‑1 (MQ=255)
                                                                                                              aGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGc                                        <  1:833149/150‑1 (MQ=255)
                                                                                                                cgcgTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCAt                                      >  1:1009098/1‑150 (MQ=255)
                                                                                                                      cAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGc                                >  1:1078563/1‑150 (MQ=255)
                                                                                                                      cAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGc                                >  2:382770/1‑150 (MQ=255)
                                                                                                                       aGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGcc                               <  2:32819/150‑1 (MQ=255)
                                                                                                                       aGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGcc                               >  2:2158416/1‑150 (MQ=255)
                                                                                                                        gAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCg                              <  2:1320537/150‑1 (MQ=255)
                                                                                                                        gAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCg                              >  2:1926400/1‑150 (MQ=255)
                                                                                                                         aaCGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCgg                             <  2:1009098/150‑1 (MQ=255)
                                                                                                                          aCGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGt                            <  1:1752178/150‑1 (MQ=255)
                                                                                                                          aCGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGt                            >  1:1507853/1‑150 (MQ=255)
                                                                                                                               aaGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTcgcg                       <  2:935559/150‑1 (MQ=255)
                                                                                                                                aGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCgg                      <  2:2380236/150‑1 (MQ=255)
                                                                                                                                 gCACGTGAGTCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCTAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCATGAATTCTGGCCGGTTCGCggg                     <  2:684055/150‑1 (MQ=255)
                                                                                                                                 gCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCggg                     >  2:1366133/1‑150 (MQ=255)
                                                                                                                                 gCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCggg                     <  1:517112/150‑1 (MQ=255)
                                                                                                                                   aCGGGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGcc                   <  1:1571843/150‑1 (MQ=255)
                                                                                                                                     gTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGa                 >  1:1915053/1‑150 (MQ=255)
                                                                                                                                     gTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGa                 >  1:50677/1‑150 (MQ=255)
                                                                                                                                        tctttctggCGGTCGCGAACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACGc              >  2:708703/8‑150 (MQ=255)
                                                                                                                                         ggCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAACCGAAAAATCGCTGATGTCCCGTACCGCTCGTCTGGGCGCGCCGCGAGGACAGGCATTCTGGCCGGTTCGCGGGCCGACGCt             >  2:747586/1‑150 (MQ=255)
                                                                                                                                           ttctggCGGTCGCGAACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACGCTGc           >  2:2175553/5‑150 (MQ=255)
                                                                                                                                           ttctggCGGTCGCGAACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCCAAAAATCGCTCATGTCCCGTACCGGTGGTCTGGCAGCGCCGCGTGGACAGCCATTCTGGCCGGTTCGCGGGCCGACGCTGc           >  2:2214966/5‑150 (MQ=255)
                                                                                                                                           cccAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACGCTGc           >  1:211501/1‑150 (MQ=255)
                                                                                                                                              aGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACGCTGCATc        <  1:1725953/150‑1 (MQ=255)
                                                                                                                                              aGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACGCTGCATc        <  1:1225359/150‑1 (MQ=255)
                                                                                                                                                gCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACGCTGCATCGc      <  2:1965381/150‑1 (MQ=255)
                                                                                                                                                  ggTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACGCTGCATCGCTa    <  1:263216/150‑1 (MQ=255)
                                                                                                                                                    ttCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACGCTGCATCGCTATg  >  1:38412/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CTGGCAGGGCTGGAGTCTTCCATCCAGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACGCTGCATCGCTATG  >  NC_000913/3787471‑3787768

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: