Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,183,345:1 +T 5.8% intergenic (+22/+468) zupT → / ← ribB zinc transporter/3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,183,3451.T 5.8% 472.1 / 23.8 120intergenic (+22/+468)zupT/ribBzinc transporter/3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base . (60/53);  new base T (4/3);  total (64/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

Only 100 of 125 total aligned reads displayed.

CGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAATAACGTTTCAAATGCAG  >  NC_000913/3183196‑3183494
                                                                                                                                                      |                                                                                                                                                     
cGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCa                                                                                                                                                        >  2:1716941/1‑150 (MQ=255)
 ggTTGATGAATTAATGCCGCTCGCCAAAGTGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAg                                                                                                                                                      <  2:1796978/150‑1 (MQ=255)
    tgatgaATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAgcgc                                                                                                                                                   >  2:715939/1‑150 (MQ=255)
     gattaATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCa                                                                                                                                                  <  1:912389/150‑1 (MQ=255)
      atgaATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTTCAAACGGCGGGAATTGGTTATATACAACATCAGAAGTTGTCCAGCGCAt                                                                                                                                                 <  2:1097373/150‑1 (MQ=255)
       tgaATTAATGCCGCTCGCCAAAGAGATTGACACTAATAATAACCCCAGCTATGGAGTACTATGTGTAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATg                                                                                                                                                <  2:1589159/150‑1 (MQ=255)
         aaTTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGa                                                                                                                                              >  1:647553/1‑150 (MQ=255)
         aaTTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGAGATGTGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGa                                                                                                                                              >  2:1373684/1‑150 (MQ=255)
          aTTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGTATTGGTTAAAAACAACATCAGAAGTTGTCAAGCGCATGGAc                                                                                                                                             >  1:1195524/1‑150 (MQ=255)
          aTTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGAc                                                                                                                                             >  2:347606/1‑150 (MQ=255)
           ttAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACTAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACa                                                                                                                                            <  1:372095/150‑1 (MQ=255)
            tAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACaa                                                                                                                                           >  2:928756/1‑150 (MQ=255)
            tAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACaa                                                                                                                                           >  1:1002934/1‑150 (MQ=255)
             aaTGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAAc                                                                                                                                          >  1:13009/1‑150 (MQ=255)
              aTGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTGTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACt                                                                                                                                         <  1:784194/150‑1 (MQ=255)
               tGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACtt                                                                                                                                        >  2:1535740/1‑150 (MQ=255)
                gCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGAGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTACAGCGCATGGACAACTTc                                                                                                                                       >  2:84367/1‑150 (MQ=255)
                  cGCTCGCCAAAGAGATTGCCCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTAAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTg                                                                                                                                     <  1:1216913/150‑1 (MQ=255)
                   gCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACTTCTga                                                                                                                                    <  2:899550/150‑2 (MQ=255)
                    cTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGCTGGACAACTTCTgat                                                                                                                                   <  2:1360251/150‑3 (MQ=255)
                    cTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGcc                                                                                                                                   <  2:414066/150‑1 (MQ=255)
                     tCGCCAAAGAGATTGACCCTAATAATAAACCAAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGTTATTGGTTAAAAACAACAGCAGAAGTTGTCCATGCGCTGGACAACTTCTgatg                                                                                                                                  <  1:1281090/150‑4 (MQ=255)
                      cGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAg                                                                                                                                 >  2:1744908/1‑150 (MQ=255)
                       gCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCACAAGTTGTCCAGCGCATGGACAACTTCTGCCAGc                                                                                                                                >  2:297835/1‑150 (MQ=255)
                        ccAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCt                                                                                                                               >  2:500356/1‑150 (MQ=255)
                          aaaGAGATTGACCCTAATCATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGGGATGGGAGTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGGTAAACACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTc                                                                                                                             <  1:469249/150‑1 (MQ=255)
                            agagATTGAACCTAATAATAACCCCCGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAACCAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCAc                                                                                                                           <  2:363268/150‑1 (MQ=255)
                                   gACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGcc                                                                                                                    >  1:957866/1‑150 (MQ=255)
                                   gACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGcc                                                                                                                    >  1:510347/1‑150 (MQ=255)
                                   gACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGcc                                                                                                                    <  1:1151689/150‑1 (MQ=255)
                                     cccTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCtt                                                                                                                  >  2:1391985/1‑150 (MQ=255)
                                     cccTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCtt                                                                                                                  <  2:28201/150‑1 (MQ=255)
                                          ataataaCCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCt                                                                                                             <  2:1536337/150‑1 (MQ=255)
                                            aataaCCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGATTTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTAc                                                                                                           >  2:1065108/1‑150 (MQ=255)
                                              taaCCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCt                                                                                                         >  1:1404046/1‑150 (MQ=255)
                                                  cccAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTTCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACATTTTc                                                                                                     >  2:1622349/1‑150 (MQ=255)
                                                  cccAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTc                                                                                                     <  2:1448181/150‑1 (MQ=255)
                                                     aGCTATGGCGTACTATGTTGAATGTCAGTTATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGc                                                                                                  <  1:362467/150‑1 (MQ=255)
                                                     aGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGc                                                                                                  <  2:8810/150‑1 (MQ=255)
                                                      gCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCt                                                                                                 >  2:1599006/1‑150 (MQ=255)
                                                      gCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCt                                                                                                 >  2:74594/1‑150 (MQ=255)
                                                       cTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTg                                                                                                <  2:792233/150‑1 (MQ=255)
                                                          tGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACATTTTCGGCTGtct                                                                                             >  2:1510940/1‑150 (MQ=255)
                                                              gTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctcttc                                                                                         <  2:647553/150‑1 (MQ=255)
                                                                   aTGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCg                                                                                    <  2:1220155/150‑1 (MQ=255)
                                                                    tgtgGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGt                                                                                   >  2:836928/1‑150 (MQ=255)
                                                                       ggAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACt                                                                                >  1:537821/1‑150 (MQ=255)
                                                                            gtgaGTGATGGGATTCAGTTTAGTGCTGCTACAAACTGAGGGAATTGGTTAAAAACATCATCAGTAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTa                                                                           <  2:480381/147‑1 (MQ=255)
                                                                            gTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTa                                                                           <  1:1397714/150‑1 (MQ=255)
                                                                             tCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTaa                                                                          <  1:624775/150‑1 (MQ=255)
                                                                              cAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTaat                                                                         >  2:162554/1‑149 (MQ=255)
                                                                              cAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAg                                                                         <  2:1512741/150‑1 (MQ=255)
                                                                                 tGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTAt                                                                      <  1:846072/150‑1 (MQ=255)
                                                                                    tGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTg                                                                   <  1:514991/150‑1 (MQ=255)
                                                                                        aTTCAGTTTAGTGCTGCTACAACCGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCGAATTCTACCGTTTAGGCTGTCTCTACTCTCGTACTGTTTAAGGATTTGTTcc                                                               >  1:351170/1‑150 (MQ=255)
                                                                                         ttCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCcg                                                              <  2:1210053/150‑1 (MQ=255)
                                                                                            aGTTTAGTGCTGCTACAATCGGCGGTAATTGGTTAAAAACAACATCATAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCgcgc                                                           <  1:157632/150‑1 (MQ=255)
                                                                                            aGTTAAGTGCGGCAACCAACGGCGGGAATTGGCTAAAAACAACATCAGACGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTAATCTCCCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCgcgc                                                           <  2:248544/150‑1 (MQ=255)
                                                                                             gTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCTCGCt                                                          >  2:808841/1‑150 (MQ=255)
                                                                                             gTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCATCACGAAAGCCTTATTCTACCATTACGGCTGTCTCTGCTCTCGTACTGTTTAAGTATATGTTCCgagca                                                          >  1:1633393/1‑149 (MQ=255)
                                                                                                tAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCa                                                       <  1:1535740/150‑1 (MQ=255)
                                                                                                   tgctgcTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGc                                                    <  1:1316681/150‑1 (MQ=255)
                                                                                                         tACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCg                                              <  1:1744908/150‑1 (MQ=255)
                                                                                                             aaCGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGtctc                                          <  1:1396873/150‑1 (MQ=255)
                                                                                                                  cggGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTTTTCCGCGATTCAAGCAGTGCGTCTCTGTcc                                     >  2:1338279/1‑150 (MQ=255)
                                                                                                                  cggGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTcc                                     >  1:1070637/1‑150 (MQ=255)
                                                                                                                   gggAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCt                                    <  1:1297840/150‑1 (MQ=255)
                                                                                                                    ggAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTg                                   >  1:1660913/1‑150 (MQ=255)
                                                                                                                     gAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTgg                                  >  1:1742340/1‑150 (MQ=255)
                                                                                                                       aTTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCg                                >  2:663800/1‑150 (MQ=255)
                                                                                                                       aTTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCg                                <  2:1801694/150‑1 (MQ=255)
                                                                                                                       aTTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGAAGTGCGTCTCTGTCCTGGCg                                <  1:1505169/150‑1 (MQ=255)
                                                                                                                          ggTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTGTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGaaa                             >  2:1444881/1‑150 (MQ=255)
                                                                                                                          ggTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTATGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCACTTCTCTCGTACTGTTTAAGTATTTGTTGCGCGCTTCAAGCAGTGCGTCTCAGTCCTGGCGaaa                             >  2:785177/1‑150 (MQ=255)
                                                                                                                           gTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAAc                            <  1:1622349/150‑1 (MQ=255)
                                                                                                                           gTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAAc                            >  2:1355358/1‑150 (MQ=255)
                                                                                                                           gTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACATTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAAc                            >  1:1015321/1‑150 (MQ=255)
                                                                                                                           gTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGACAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTATTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTATCTGTCCTGGCGAAAc                            >  1:882749/1‑150 (MQ=255)
                                                                                                                               aaaaCAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCGTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGAGCGTCTCTGTCCTGGCGAAACGccc                        >  1:1714354/1‑150 (MQ=255)
                                                                                                                                aaaCAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCg                       >  1:1324312/1‑150 (MQ=255)
                                                                                                                                aaaCAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCg                       <  1:505868/150‑1 (MQ=255)
                                                                                                                                   caacaTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACa                    >  1:1070298/1‑150 (MQ=255)
                                                                                                                                     acaTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCTAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAAt                  <  1:1800226/150‑1 (MQ=255)
                                                                                                                                      caTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGAAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAATa                 >  1:969618/1‑150 (MQ=255)
                                                                                                                                       aTCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAATaa                >  2:562033/1‑150 (MQ=255)
                                                                                                                                        tCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTTCTGTTTAAGTTTTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAATAAc               >  1:1737182/1‑150 (MQ=255)
                                                                                                                                        tCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAATAAc               >  2:894596/1‑150 (MQ=255)
                                                                                                                                        tCAGAAGTTGTCCAGCGAATGGACAAATTCTGCAAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAATAAc               <  2:351170/150‑1 (MQ=255)
                                                                                                                                        gcAGAAGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAATAAc               >  2:1754071/2‑150 (MQ=255)
                                                                                                                                             aGTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAATAACGTTTc          >  2:873111/1‑150 (MQ=255)
                                                                                                                                             aGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAATAACGTTTc          >  1:1737581/1‑150 (MQ=255)
                                                                                                                                              gTTGTCCATGCGCTGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGGAACGCCCGACAATAACGTTTCa         >  1:925834/1‑150 (MQ=255)
                                                                                                                                                tGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAATAACGTTTCaaa       >  1:1673870/1‑150 (MQ=255)
                                                                                                                                                tGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAATAACGTTTCaaa       <  2:1735343/150‑1 (MQ=255)
                                                                                                                                                 gTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAATAACGTTTCAAAt      <  2:537821/150‑1 (MQ=255)
                                                                                                                                                  tCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAATAACGTTTCAAATg     <  2:1404046/150‑1 (MQ=255)
                                                                                                                                                   ccAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAATAACGTTTCAAATGc    <  1:1374833/150‑1 (MQ=255)
                                                                                                                                                   ccAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAATAACGTTTCAAATGc    >  1:637970/1‑150 (MQ=255)
                                                                                                                                                    cAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAATAACGTTTCAAATGCa   <  1:84367/150‑1 (MQ=255)
                                                                                                                                                     aGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAATAACGTTTCAAATGCAg  <  2:411663/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CGGTTGATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCAGCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAATAACGTTTCAAATGCAG  >  NC_000913/3183196‑3183494

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