Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,080,041 T→G 7.5% intergenic (+72/+258) yiiF → / ← fdhE putative thymol sensitivity protein, CopG family putative transcriptional regulator/formate dehydrogenase formation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,080,0410TG 7.5% 427.4 / 19.3 147intergenic (+72/+258)yiiF/fdhEputative thymol sensitivity protein, CopG family putative transcriptional regulator/formate dehydrogenase formation protein
Reads supporting (aligned to +/- strand):  ref base T (62/74);  new base G (6/5);  total (68/79)
Fisher's exact test for biased strand distribution p-value = 7.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

Only 100 of 149 total aligned reads displayed.

ACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTGAGGCCAAACC  >  NC_000913/4079892‑4080189
                                                                                                                                                     |                                                                                                                                                    
      gtTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                                >  2:133607/1‑150 (MQ=255)
       ttCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                               >  1:557738/1‑150 (MQ=255)
       ttCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                               <  2:763167/150‑1 (MQ=39)
       ttCCTGGCAGCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                               <  2:436466/150‑1 (MQ=255)
         gcTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTGGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGca                                                                                                                                             <  2:472153/149‑2 (MQ=38)
         ccTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                             <  2:1725340/150‑1 (MQ=255)
         ccTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGca                                                                                                                                             <  2:134935/150‑2 (MQ=38)
         ccTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATAGCAGCGTCAGCTGCGCGAGGAATGGTACTGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                             <  1:1124112/150‑1 (MQ=255)
          cTGGCATCTGGCAAGGGGGTGAGGAAGATGGTGTCGAATATCACCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCt                                                                                                                                            <  2:911907/150‑1 (MQ=255)
           tGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                                                                                                           <  1:452228/150‑1 (MQ=255)
            ggCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                                                                                                          >  1:441810/1‑150 (MQ=255)
             gCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACa                                                                                                                                         <  2:1308734/150‑1 (MQ=255)
                 cTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGt                                                                                                                                     >  2:19900/1‑150 (MQ=255)
                 cTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGt                                                                                                                                     >  2:267403/1‑150 (MQ=255)
                 cTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGt                                                                                                                                     <  1:1274557/150‑1 (MQ=255)
                 cTGGCAAGGGTGTGAGGAAGATGGTGGCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGt                                                                                                                                     <  1:1281367/150‑1 (MQ=255)
                   ggCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCg                                                                                                                                   >  2:1688019/1‑150 (MQ=255)
                   ggCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCg                                                                                                                                   >  2:66127/1‑150 (MQ=255)
                      aggggTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGtat                                                                                                                                <  2:531769/148‑1 (MQ=255)
                        gggTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATaa                                                                                                                              >  1:774128/1‑150 (MQ=255)
                                ggAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATg                                                                                                                      >  2:149058/1‑150 (MQ=255)
                                    gATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTAtt                                                                                                                  <  1:1106648/150‑1 (MQ=255)
                                      tggtggAGAATATCAGCGTACGCTGCGCGAGGCATGGTAATGGGAAAAAATGATTTCAATCAGATTGCTGACACCGTGCGCGTTGTTCATGCCGGAGGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGc                                                                                                                <  2:89832/150‑1 (MQ=25)
                                      tGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGc                                                                                                                <  1:6204/150‑1 (MQ=255)
                                      gggTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGc                                                                                                                <  2:1453231/149‑1 (MQ=255)
                                            cGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTa                                                                                                          >  1:343054/1‑150 (MQ=255)
                                             gAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTAt                                                                                                         <  1:669906/150‑1 (MQ=255)
                                             gAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTAt                                                                                                         >  2:1673418/1‑150 (MQ=255)
                                                 atCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGt                                                                                                     >  1:794224/1‑150 (MQ=255)
                                                  tCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTa                                                                                                    <  1:1713429/150‑1 (MQ=255)
                                                  tCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTa                                                                                                    >  2:1075283/1‑150 (MQ=255)
                                                   cAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCGTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAg                                                                                                   >  2:1382836/1‑150 (MQ=255)
                                                    aGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCGGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAgg                                                                                                  <  2:410793/150‑1 (MQ=255)
                                                    aGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAgg                                                                                                  >  1:1522723/1‑150 (MQ=255)
                                                     gCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGCCAACGTGCGAGTTGTTCATGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGGTTATTGGTAGGc                                                                                                 >  2:1073696/1‑150 (MQ=255)
                                                     gCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGc                                                                                                 <  1:1609132/150‑1 (MQ=255)
                                                      cGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGcc                                                                                                <  1:65800/150‑1 (MQ=255)
                                                       gTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCt                                                                                               <  2:1312058/150‑1 (MQ=255)
                                                        tAAGATGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTg                                                                                              <  2:589524/150‑1 (MQ=255)
                                                           tcTGCGCGAGGAATGTTAATGGGTATAAATTATTTCAATCAGATTGCTGAGAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGACTACAAAGTAGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATa                                                                                           <  1:302746/149‑1 (MQ=25)
                                                            cTGCGCGAGGAATGGTAATGGGAATAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATaa                                                                                          <  2:859894/150‑1 (MQ=255)
                                                               cgcgAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACACCGTGCGCGTTGTTCATGCCGGATGCGGCGTGCACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAagcg                                                                                       >  2:277514/1‑150 (MQ=255)
                                                               cgcgAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAagcg                                                                                       >  1:1343389/1‑150 (MQ=255)
                                                                gcgAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAagcgc                                                                                      >  1:1012328/1‑150 (MQ=255)
                                                                gcgAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGACTACAAAGTCGAATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAaacgc                                                                                      >  1:1106329/1‑150 (MQ=255)
                                                                   aGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAagcgcagc                                                                                   >  1:445576/1‑150 (MQ=255)
                                                                   aGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAagcgcagc                                                                                   >  1:1596165/1‑150 (MQ=255)
                                                                    ggtattgTATTGGGAAAAAATGTTGTCAATCAGATTGCTTACAATGTGCGCGTGGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAGAAAGTCGTTAAATTTCAATGTATTGCTGTTTATTTGGAGGCCTGATAagcgcagcg                                                                                  <  2:1141393/144‑1 (MQ=21)
                                                                    ggAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAagcgcagcg                                                                                  >  1:1500054/1‑150 (MQ=255)
                                                                     gAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAagcgcagcgc                                                                                 <  2:155886/150‑1 (MQ=255)
                                                                     gAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAagcgcagcgc                                                                                 <  1:1041648/150‑1 (MQ=255)
                                                                      aaTGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAgcgcagcgca                                                                                <  1:891659/150‑1 (MQ=255)
                                                                         ggTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCa                                                                             <  2:74241/150‑1 (MQ=255)
                                                                            aaTGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGc                                                                          <  1:1019717/150‑1 (MQ=255)
                                                                            aaTGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGc                                                                          >  2:1068132/1‑150 (MQ=255)
                                                                               gggTAAAAATGATGTCAATCAGATTGCAGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATGCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGGAGCGCATCAGGCAAt                                                                       <  2:1503732/150‑1 (MQ=255)
                                                                                 gAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCt                                                                     <  1:1395012/150‑1 (MQ=255)
                                                                                 gAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCt                                                                     >  1:470298/1‑150 (MQ=255)
                                                                                    aaaaTGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTAc                                                                  >  1:1452383/1‑150 (MQ=255)
                                                                                         gatgTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTg                                                             <  2:441810/150‑1 (MQ=255)
                                                                                            gTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGtca                                                          >  2:316621/1‑150 (MQ=255)
                                                                                                aTCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAg                                                      <  2:1216204/150‑1 (MQ=255)
                                                                                                aTCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAg                                                      <  1:267403/150‑1 (MQ=255)
                                                                                                 tCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGt                                                     <  1:133607/150‑1 (MQ=255)
                                                                                                  cAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCAGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGtc                                                    <  1:1499262/150‑1 (MQ=255)
                                                                                                    gATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGtctc                                                  >  1:900809/1‑150 (MQ=255)
                                                                                                      ttGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCAt                                                <  1:1102921/150‑1 (MQ=255)
                                                                                                       catgaaCAACGCGCACGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATg                                               >  1:1522677/6‑150 (MQ=255)
                                                                                                         cTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGcc                                             >  2:396601/1‑150 (MQ=255)
                                                                                                         cTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGcc                                             <  1:293899/150‑1 (MQ=255)
                                                                                                         cTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTAATGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGTCAATCTTACGCTTGTCATCAGTCTCATGcc                                             >  2:1577574/1‑150 (MQ=255)
                                                                                                           gACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGGAGTTTATTGGTAGGCCTGATAAGCGCCGCGCATCAGGCATTCTTACGCTTGTCTTCAGTCTCATGCCCa                                           >  2:1196854/1‑150 (MQ=255)
                                                                                                           gACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCa                                           >  1:1114690/1‑150 (MQ=255)
                                                                                                            aCAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCAt                                          >  1:1071739/1‑150 (MQ=255)
                                                                                                                cGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAAcc                                      >  2:637172/1‑150 (MQ=255)
                                                                                                                 gTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCt                                     <  1:6228/150‑1 (MQ=255)
                                                                                                                   gcgcGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTa                                   >  1:920782/1‑150 (MQ=255)
                                                                                                                    cgcgTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTAt                                  <  2:586273/150‑1 (MQ=255)
                                                                                                                    cgcgTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTATATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTAt                                  <  1:982745/150‑1 (MQ=255)
                                                                                                                      cgTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTTTTGGTAGGCCTGTTAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTg                                <  2:1371756/150‑1 (MQ=255)
                                                                                                                      cgTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTg                                >  1:513914/1‑150 (MQ=255)
                                                                                                                         tgttCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGGCGTATAAATTAAATGTGTTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCAGCAGGCAATCTTACGCTTGTCATCAGGCTCATGCCCATAACCTTATTGAtg                             <  2:720536/150‑1 (MQ=37)
                                                                                                                         tgttCATGCCGGAGGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTGTTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGAtg                             <  2:714184/150‑1 (MQ=255)
                                                                                                                             taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTTTTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATgtggt                         >  1:1167441/4‑150 (MQ=38)
                                                                                                                             taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAAGTCAATGTATTGCAGTTTAGTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATttggt                         >  1:1240922/4‑150 (MQ=37)
                                                                                                                             cATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTAtt                                 >  1:360535/1‑142 (MQ=255)
                                                                                                                             cATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATgtggt                         <  1:1601800/150‑1 (MQ=255)
                                                                                                                              aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATgtggtg                        >  2:1144958/3‑150 (MQ=38)
                                                                                                                              aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTTAAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATATTACGCTTGTCATCAGTCTCATGCCCATATCCTTATTGATgtggtg                        >  2:354078/1‑150 (MQ=255)
                                                                                                                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTTCGCTGGTCATCAGTCTCATGCCCATAACCTTAGTGATgtggtgg                       >  2:121995/2‑150 (MQ=37)
                                                                                                                                 ccGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCCGTTTATTGGTAGGCCTTATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATgtggtggtg                     <  2:554394/150‑1 (MQ=255)
                                                                                                                                        cggcgTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAAt              >  1:38701/1‑150 (MQ=255)
                                                                                                                                           cgTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTGa           >  2:935247/1‑150 (MQ=255)
                                                                                                                                           cgTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGAAATCTTACGCTTGTCATCAGTCTCATGCCCATACCCTTAGTGATGTGGTGGGGGGTAAATTGa           >  2:1727333/1‑150 (MQ=255)
                                                                                                                                             tGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTGAgg         >  2:193644/1‑150 (MQ=255)
                                                                                                                                             tGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTGAgg         >  2:1546200/1‑150 (MQ=255)
                                                                                                                                              gAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTGAGGc        <  1:1275750/150‑1 (MQ=255)
                                                                                                                                               aaCGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGACGTGGTGGTGGGTAAATTGAGGcc       <  2:451734/150‑1 (MQ=255)
                                                                                                                                               aaCGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATGGATGTGGTGGTGGGTAAATTGAGGcc       >  2:1273163/1‑150 (MQ=255)
                                                                                                                                                  gCCTTATCCGGCCTACAAAGTCGTATAAATTCAAGGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTGAGGCCaaa    >  1:1292726/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTGAGGCCAAACC  >  NC_000913/4079892‑4080189

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