Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,080,042 T→C 6.8% intergenic (+73/+257) 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,0420TC 6.8% 435.1 / 16.3 147intergenic (+73/+257)yiiF/fdhEputative thymol sensitivity protein, CopG family putative transcriptional regulator/formate dehydrogenase formation protein
Reads supporting (aligned to +/- strand):  ref base T (63/74);  new base C (5/5);  total (68/79)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.51e-01

Only 100 of 148 total aligned reads displayed.

CAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTGAGGCCAAACCGGC  >  NC_000913/4079894‑4080192
                                                                                                                                                    |                                                                                                                                                      
cAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTa                                                                                                                                                       <  1:1541630/150‑1 (MQ=255)
    gtTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                                   >  2:133607/1‑150 (MQ=255)
     ttCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                                  <  2:763167/150‑1 (MQ=39)
     ttCCTGGCAGCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                                  <  2:436466/150‑1 (MQ=255)
      tCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                                 <  1:1520399/150‑1 (MQ=255)
      tCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                                 <  1:1278300/150‑1 (MQ=255)
       ccTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                                <  1:856985/150‑1 (MQ=255)
       ccTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATAGCAGCGTCAGCTGCGCGAGGAATGGTACTGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                                <  1:1124112/150‑1 (MQ=255)
        cTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGcat                                                                                                                                               <  1:1391596/150‑3 (MQ=38)
        cTGGCATCTGGCAAGGGGGTGAGGAAGATGGTGTCGAATATCACCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCt                                                                                                                                               <  2:911907/150‑1 (MQ=255)
         tGGCATCTGGCAAGGGTGTGAGGAAGATGTTGTCGAATATCAGTGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCTGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGcatg                                                                                                                                              <  2:193048/150‑4 (MQ=38)
         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)
                 ggCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCg                                                                                                                                      >  2:1688019/1‑150 (MQ=255)
                    aggggTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGtat                                                                                                                                   <  2:531769/148‑1 (MQ=255)
                    aaGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGtat                                                                                                                                   <  1:1482985/150‑1 (MQ=255)
                     aGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGtata                                                                                                                                  <  1:1252619/150‑1 (MQ=255)
                      gggTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATaa                                                                                                                                 >  1:774128/1‑150 (MQ=255)
                        gtgtGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTAGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCCTATAAAt                                                                                                                               >  2:449614/1‑150 (MQ=255)
                          gtgAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTc                                                                                                                             >  1:902317/1‑150 (MQ=255)
                              ggAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATg                                                                                                                         >  2:1275750/1‑150 (MQ=255)
                              ggAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATg                                                                                                                         >  2:149058/1‑150 (MQ=255)
                                  gATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTAtt                                                                                                                     <  1:1106648/150‑1 (MQ=255)
                                  gATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTAtt                                                                                                                     <  2:505462/150‑1 (MQ=255)
                                    tggtggAGAATATCAGCGTACGCTGCGCGAGGCATGGTAATGGGAAAAAATGATTTCAATCAGATTGCTGACACCGTGCGCGTTGTTCATGCCGGAGGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGc                                                                                                                   <  2:89832/150‑1 (MQ=25)
                                    tGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGc                                                                                                                   <  2:41416/150‑1 (MQ=255)
                                    tGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGc                                                                                                                   <  1:6204/150‑1 (MQ=255)
                                      gtgtCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAg                                                                                                                 >  2:1114079/1‑150 (MQ=255)
                                           gAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTAt                                                                                                            >  2:1716726/1‑150 (MQ=255)
                                             atatCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTg                                                                                                          >  2:249930/1‑150 (MQ=255)
                                                tCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTa                                                                                                       >  2:1075283/1‑150 (MQ=255)
                                                tCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTa                                                                                                       <  1:1713429/150‑1 (MQ=255)
                                                  aGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAgg                                                                                                     >  1:1522723/1‑150 (MQ=255)
                                                    cGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACACCGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGcc                                                                                                   >  2:1217199/1‑150 (MQ=255)
                                                    cGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGcc                                                                                                   <  1:65800/150‑1 (MQ=255)
                                                     gTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGCCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCt                                                                                                  <  1:472288/150‑1 (MQ=255)
                                                         tcTGCGCGAGGAATGTTAATGGGTATAAATTATTTCAATCAGATTGCTGAGAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGACTACAAAGTAGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATa                                                                                              <  1:302746/149‑1 (MQ=25)
                                                          cTGCGCGAGGAATGGTAATGGGAATAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATaa                                                                                             <  2:859894/150‑1 (MQ=255)
                                                            gcgcGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAagc                                                                                           <  1:409442/150‑1 (MQ=255)
                                                             cgcgAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAagcg                                                                                          >  1:1343389/1‑150 (MQ=255)
                                                              gcgAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAagcgc                                                                                         >  1:1012328/1‑150 (MQ=255)
                                                              gcgAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGACTACAAAGTCGAATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAaacgc                                                                                         >  1:1106329/1‑150 (MQ=255)
                                                                  ggtattgTATTGGGAAAAAATGTTGTCAATCAGATTGCTTACAATGTGCGCGTGGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAGAAAGTCGTTAAATTTCAATGTATTGCTGTTTATTTGGAGGCCTGATAagcgcagcg                                                                                     <  2:1141393/144‑1 (MQ=21)
                                                                  ggAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAagcgcagcg                                                                                     >  1:1500054/1‑150 (MQ=255)
                                                                   gAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAagcgcagcgc                                                                                    <  1:1041648/150‑1 (MQ=255)
                                                                    aaTGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAgcgcagcgca                                                                                   <  1:891659/150‑1 (MQ=255)
                                                                    aaTGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAgcgcagcgca                                                                                   >  1:686650/1‑150 (MQ=255)
                                                                       ggTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCa                                                                                <  2:74241/150‑1 (MQ=255)
                                                                       ggTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCa                                                                                >  1:1718727/1‑150 (MQ=255)
                                                                          aaTGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGc                                                                             <  1:1019717/150‑1 (MQ=255)
                                                                          aaTGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGc                                                                             >  2:1068132/1‑150 (MQ=255)
                                                                                  aaaaTGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTAc                                                                     <  2:557738/150‑1 (MQ=255)
                                                                                  aaaaTGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTAc                                                                     >  1:1452383/1‑150 (MQ=255)
                                                                                         tgTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGtc                                                              >  2:1006812/1‑150 (MQ=255)
                                                                                          gTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGtca                                                             >  2:316621/1‑150 (MQ=255)
                                                                                             aaTCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGtcatca                                                          >  1:601338/1‑150 (MQ=255)
                                                                                              aTCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAg                                                         <  1:267403/150‑1 (MQ=255)
                                                                                                cAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGtc                                                       >  2:1498296/1‑150 (MQ=255)
                                                                                                cAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCAGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGtc                                                       >  1:629271/1‑150 (MQ=255)
                                                                                                    ttGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCAt                                                   <  1:1102921/150‑1 (MQ=255)
                                                                                                     tGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATg                                                  <  1:446356/150‑1 (MQ=255)
                                                                                                     tGCTGAAAACGTGCGCGGTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATg                                                  <  2:635976/150‑1 (MQ=255)
                                                                                                       cttaCAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGcc                                                <  2:1548213/147‑1 (MQ=255)
                                                                                                       cTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGcc                                                >  2:396601/1‑150 (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:1415043/1‑150 (MQ=255)
                                                                                                              cGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAAcc                                         >  2:637172/1‑150 (MQ=255)
                                                                                                              cGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAAcc                                         <  2:1324129/150‑1 (MQ=255)
                                                                                                                 gcgcGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTa                                      >  1:920782/1‑150 (MQ=255)
                                                                                                                  cgcgTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTAt                                     <  2:586273/150‑1 (MQ=255)
                                                                                                                  cgagttgttCATTCCGGATGCTGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTAt                                     <  1:356304/147‑1 (MQ=255)
                                                                                                                    cgTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTTTTGGTAGGCCTGTTAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTg                                   <  2:1371756/150‑1 (MQ=255)
                                                                                                                    cgTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTg                                   >  1:513914/1‑150 (MQ=255)
                                                                                                                       tgttCATGCCGGAGGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTGTTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGAtg                                <  2:714184/150‑1 (MQ=255)
                                                                                                                           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAAGTCAATGTATTGCAGTTTAGTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATttggt                            >  1:1240922/4‑150 (MQ=37)
                                                                                                                           cATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTAtt                                    >  1:360535/1‑142 (MQ=255)
                                                                                                                           cATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTAtt                                    <  2:360535/142‑1 (MQ=255)
                                                                                                                           cATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATgtggt                            <  1:1608088/150‑1 (MQ=255)
                                                                                                                           cATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATgtggt                            <  1:1601800/150‑1 (MQ=255)
                                                                                                                           aaTGCCGGATTCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATgtggt                            <  1:131568/149‑1 (MQ=255)
                                                                                                                            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATgtggtg                           >  2:1163331/3‑150 (MQ=38)
                                                                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTTCGCTGGTCATCAGTCTCATGCCCATAACCTTAGTGATgtggtgg                          >  2:121995/2‑150 (MQ=37)
                                                                                                                               ccGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATgtggtggtg                        <  2:1619837/150‑1 (MQ=255)
                                                                                                                                 ggATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCTCATAACCTTATTGATGTGGTGGTggg                      <  1:29567/150‑1 (MQ=255)
                                                                                                                                        gcgTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTg               <  2:686650/150‑1 (MQ=255)
                                                                                                                                         cgTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTGa              >  2:935247/1‑150 (MQ=255)
                                                                                                                                         cgTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGAAATCTTACGCTTGTCATCAGTCTCATGCCCATACCCTTAGTGATGTGGTGGGGGGTAAATTGa              >  2:1727333/1‑150 (MQ=255)
                                                                                                                                           tGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTGAgg            >  2:1546200/1‑150 (MQ=255)
                                                                                                                                            gAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTGAGGc           >  2:30360/1‑150 (MQ=255)
                                                                                                                                            gAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTGAGGc           <  1:1275750/150‑1 (MQ=255)
                                                                                                                                             aaCGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGACGTGGTGGTGGGTAAATTGAGGcc          <  2:451734/150‑1 (MQ=255)
                                                                                                                                                gCCTTATCCGGCCTACAAAGTCGTATAAATTCAAGGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTGAGGCCaaa       >  1:1292726/1‑150 (MQ=255)
                                                                                                                                                  cTTCTCCGGCCTAGAACGTAGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCGTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAAACTTATTGATGTGGTGGTGGGGAAATTGAGGACAAAcc     >  2:136770/1‑150 (MQ=38)
                                                                                                                                                    tATCCGGCCTACAAAGTAGTATAATTTCAATTTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTGAGGCCAAACCGGc  <  2:1443069/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTGAGGCCAAACCGGC  >  NC_000913/4079894‑4080192

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: