Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,606,642 C→T 5.9% intergenic (‑241/+27) leuQ ← / ← rsmC tRNA‑Leu/16S rRNA m(2)G1207 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,606,6420CT 5.9% 271.0 / 4.3 85intergenic (‑241/+27)leuQ/rsmCtRNA‑Leu/16S rRNA m(2)G1207 methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (42/38);  new base T (2/3);  total (44/41)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCG  >  NC_000913/4606494‑4606789
                                                                                                                                                    |                                                                                                                                                   
aaTTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCa                                                                                                                                                    <  1:690303/150‑1 (MQ=255)
aaTTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTACACGCCGCa                                                                                                                                                    >  2:1420352/1‑150 (MQ=255)
 aTTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCAt                                                                                                                                                   >  1:285530/1‑150 (MQ=255)
   tGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                                 <  1:929674/150‑1 (MQ=37)
   tGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATcc                                                                                                                                                 <  2:1083916/150‑1 (MQ=255)
   tGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATcc                                                                                                                                                 <  1:931233/150‑1 (MQ=255)
   tGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTGTAGGTCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                                <  2:1111504/150‑1 (MQ=37)
    gCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCg                                                                                                                                                >  1:142069/1‑150 (MQ=255)
     cTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGa                                                                                                                                               <  1:1242319/150‑1 (MQ=255)
      tGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGac                                                                                                                                              <  1:431154/150‑1 (MQ=255)
        gAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATGCGGCGTGAACGCCTTATCCGacct                                                                                                                                            <  2:141152/150‑3 (MQ=37)
        gAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGacac                                                                                                                                            >  1:578395/1‑150 (MQ=255)
        gAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGacac                                                                                                                                            <  2:254301/150‑1 (MQ=255)
        gAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGacac                                                                                                                                            >  1:1524659/1‑150 (MQ=255)
            aaTGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTa                                                                                                                                        <  2:310361/150‑1 (MQ=255)
               gTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTc                                                                                                                                     <  1:1390377/150‑1 (MQ=255)
                 ccATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCgg                                                                                                                                   >  2:405437/1‑150 (MQ=255)
                  cATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGc                                                                                                                                  >  1:919996/1‑150 (MQ=255)
                  cATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGc                                                                                                                                  <  2:1458430/150‑1 (MQ=255)
                  cATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGc                                                                                                                                  >  2:1088981/1‑150 (MQ=255)
                      aGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACTCCGCATCCGACACGGTATTCGGCgaga                                                                                                                              >  1:337015/1‑150 (MQ=255)
                      aGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCgaga                                                                                                                              >  1:1014592/1‑150 (MQ=255)
                            cTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTa                                                                                                                        <  2:919996/150‑1 (MQ=255)
                              aGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAAc                                                                                                                      >  2:335997/1‑150 (MQ=255)
                                 gTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCtt                                                                                                                   <  2:345813/150‑1 (MQ=255)
                                  tGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCttt                                                                                                                  >  1:924220/1‑150 (MQ=255)
                                      ggTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTttcttc                                                                                                              >  2:1237596/1‑150 (MQ=255)
                                      ggTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTttcttc                                                                                                              >  1:789065/1‑150 (MQ=255)
                                        tttGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGc                                                                                                            <  1:1019405/150‑1 (MQ=255)
                                              atgcatgcGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCg                                                                                                      >  2:1278270/1‑150 (MQ=255)
                                              atgcatgcGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCg                                                                                                      >  1:1491141/1‑150 (MQ=255)
                                                gcatgcGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCggg                                                                                                    <  1:35422/150‑1 (MQ=255)
                                                   tgcGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCa                                                                                                 <  2:142290/150‑1 (MQ=255)
                                                         ggCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATgg                                                                                           <  1:1278270/150‑1 (MQ=255)
                                                         ggCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATgg                                                                                           >  1:819354/1‑150 (MQ=255)
                                                            agatgccagatgcGACGCTGGCGGGTCTTATCTGGCCTACGATGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGgcgc                                                                                        <  2:378602/150‑1 (MQ=255)
                                                            agatgccagatgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGgcgc                                                                                        <  1:722437/150‑1 (MQ=255)
                                                                gccagatgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATa                                                                                    >  2:298346/1‑150 (MQ=255)
                                                                  cagatgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATaca                                                                                  >  1:170133/1‑150 (MQ=255)
                                                                  cagatgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATaca                                                                                  <  1:855686/150‑1 (MQ=255)
                                                                    tatgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACAcc                                                                                <  2:1055412/149‑1 (MQ=255)
                                                                    gatgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACAcc                                                                                <  1:971533/150‑1 (MQ=255)
                                                                      tgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCtt                                                                              <  1:723664/150‑1 (MQ=255)
                                                                       gcGACGCTGGCGCGTCTAATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGTATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATAAACCTTTCTTCGCCTGGCGGGTCATAATTGCGCGATACACCTTg                                                                             <  2:103098/150‑1 (MQ=255)
                                                                        cGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGa                                                                            <  2:924220/150‑1 (MQ=255)
                                                                                gcgcGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCt                                                                    >  1:1462435/1‑150 (MQ=255)
                                                                                 cgcgTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTg                                                                   >  1:702053/1‑150 (MQ=255)
                                                                                 cgcgTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTg                                                                   <  1:1401997/150‑1 (MQ=255)
                                                                                        tATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTc                                                                                                  <  2:280141/112‑1 (MQ=255)
                                                                                        tATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTc                                                                                                  >  1:280141/1‑112 (MQ=255)
                                                                                           cTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGAt                                                         >  1:1251753/1‑150 (MQ=255)
                                                                                            tGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATc                                                        >  1:1224200/1‑150 (MQ=255)
                                                                                             ggCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGTCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCa                                                       >  2:865400/1‑150 (MQ=255)
                                                                                             ggCCTACGAAGGGCTAACGTGCAGGTTTTGTAGCTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCa                                                       >  2:613395/1‑150 (MQ=255)
                                                                                               ccTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACt                                                     >  1:415280/1‑150 (MQ=255)
                                                                                                   cGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGt                                                 >  1:1419490/1‑150 (MQ=255)
                                                                                                   cGAAGGGCTAACGTGCAGGTTTTGTACGTCGGATAAGGCGTACACGCCGCATACCACACCGTATACGACGACATAAATAACCTTTCTTCGACTGGCGGGTAATAATGTCCCGATACACCTTCAATCCACCTCTTTGAGCGATCACTTcct                                                 >  2:1084944/1‑148 (MQ=255)
                                                                                                    gAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTg                                                <  1:1537860/150‑1 (MQ=255)
                                                                                                     aaGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTgg                                               >  2:651725/1‑150 (MQ=255)
                                                                                                     aaGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTgg                                               <  1:274934/150‑1 (MQ=255)
                                                                                                     aaGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTgg                                               <  2:578395/150‑1 (MQ=255)
                                                                                                           tAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCGTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCa                                         <  2:518327/150‑1 (MQ=255)
                                                                                                            aaCGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCaa                                        <  2:789065/150‑1 (MQ=255)
                                                                                                              cGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAAt                                      >  1:843263/1‑150 (MQ=255)
                                                                                                                tGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGt                                    >  1:576814/1‑150 (MQ=255)
                                                                                                                tGCAGGTTTTGTAGGTCGGATAAGGCCTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATACTGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGt                                    >  2:171249/1‑150 (MQ=255)
                                                                                                                   aGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGtctc                                 <  2:285530/150‑1 (MQ=255)
                                                                                                                   aGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGtctc                                 <  1:1021661/150‑1 (MQ=255)
                                                                                                                    ggTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCa                                >  1:1031768/1‑150 (MQ=255)
                                                                                                                        ttGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCa                            <  1:1088981/150‑1 (MQ=255)
                                                                                                                         tGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAGGCCAAATGTCTCATCCAg                           >  2:609841/1‑150 (MQ=255)
                                                                                                                          gTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGc                          >  1:1225179/1‑150 (MQ=255)
                                                                                                                          ccgtgTCGGATGCGGCGTGAACGCCTTATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGAGGAAGCCAAATGTCTCATCCAGc                          >  1:1310676/5‑150 (MQ=255)
                                                                                                                           tAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCa                         <  2:170133/150‑1 (MQ=255)
                                                                                                                            gtgTCGGATGCGGCGTGAACGCCTTATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCAc                        >  2:267734/3‑150 (MQ=255)
                                                                                                                                cGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTcc                    >  2:870083/1‑150 (MQ=255)
                                                                                                                                  gATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCgg                  >  1:324659/1‑150 (MQ=255)
                                                                                                                                    tAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGAt                <  1:569309/150‑1 (MQ=255)
                                                                                                                                      aGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATaa              >  2:729863/1‑150 (MQ=255)
                                                                                                                                       ggCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAg             >  1:427605/1‑150 (MQ=255)
                                                                                                                                        gCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAgg            <  1:60921/150‑1 (MQ=255)
                                                                                                                                          gTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCa          <  2:1200956/150‑1 (MQ=255)
                                                                                                                                            tCACGCCGCATCGGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAgg        <  1:613395/150‑1 (MQ=255)
                                                                                                                                            tCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAgg        >  1:553292/1‑150 (MQ=255)
                                                                                                                                            tCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAgg        >  1:1048244/1‑150 (MQ=255)
                                                                                                                                               cgccgcATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAg     >  1:1208193/1‑150 (MQ=255)
                                                                                                                                                  cgcATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCg  <  2:1014592/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCG  >  NC_000913/4606494‑4606789

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: