Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,080,035 A→C 7.3% intergenic (+66/+264) 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,0350AC 7.3% 411.6 / 16.5 153intergenic (+66/+264)yiiF/fdhEputative thymol sensitivity protein, CopG family putative transcriptional regulator/formate dehydrogenase formation protein
Reads supporting (aligned to +/- strand):  ref base A (64/77);  new base C (6/6);  total (70/83)
Fisher's exact test for biased strand distribution p-value = 7.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.34e-01

Only 100 of 160 total aligned reads displayed.

GCAAGAACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTGAGGCC  >  NC_000913/4079886‑4080184
                                                                                                                                                     |                                                                                                                                                     
 cAAGAACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGaa                                                                                                                                                      >  2:1608088/1‑150 (MQ=255)
 cAAGAACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGaa                                                                                                                                                      >  2:1132655/1‑150 (MQ=255)
  aaGAACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAAc                                                                                                                                                     >  1:1324129/1‑150 (MQ=255)
  aaGAACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCACGCCGCATCCGGCATGAAc                                                                                                                                                     <  1:1732403/150‑1 (MQ=255)
    gAACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGc                                                                                                                                                   >  2:1344011/1‑150 (MQ=255)
     aaCCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGcc                                                                                                                                                  <  1:1410191/150‑1 (MQ=255)
     aaCCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGAGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGAAGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGcc                                                                                                                                                  <  1:142776/150‑1 (MQ=255)
     aaCCAGTGATCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGcc                                                                                                                                                  <  1:673828/150‑1 (MQ=255)
      aCCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCt                                                                                                                                                 >  2:735315/1‑150 (MQ=255)
      aCCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCt                                                                                                                                                 >  1:1158666/1‑150 (MQ=255)
        cAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTa                                                                                                                                               <  1:1541630/150‑1 (MQ=255)
        cAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATAAGGCGTTCACGccgca                                                                                                                                               <  1:1198563/150‑4 (MQ=39)
            gtTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                           >  2:133607/1‑150 (MQ=255)
             ttCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                          >  1:557738/1‑150 (MQ=255)
             ttCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                          <  1:1104051/150‑1 (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)
               gcTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTGGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGca                                                                                                                                        <  2:472153/149‑2 (MQ=38)
               ccTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                        <  1:856985/150‑1 (MQ=255)
               ccTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATAGCAGCGTCAGCTGCGCGAGGAATGGTACTGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                        <  1:1124112/150‑1 (MQ=255)
                cTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGcat                                                                                                                                       <  1:1391596/150‑3 (MQ=38)
                 tGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                                                                                                      <  1:452228/150‑1 (MQ=255)
                       cTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGt                                                                                                                                <  1:1274557/150‑1 (MQ=255)
                       cTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGt                                                                                                                                >  2:19900/1‑150 (MQ=255)
                       cTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGt                                                                                                                                >  2:267403/1‑150 (MQ=255)
                       cTGGCAAGGGTGTGAGGAAGATGGTGGCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGt                                                                                                                                <  1:1281367/150‑1 (MQ=255)
                         ggCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCg                                                                                                                              >  2:1688019/1‑150 (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:149058/1‑150 (MQ=255)
                                      ggAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATg                                                                                                                 >  2:1275750/1‑150 (MQ=255)
                                          gATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTAtt                                                                                                             <  2:505462/150‑1 (MQ=255)
                                          gATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTAtt                                                                                                             <  1:1106648/150‑1 (MQ=255)
                                            tggtggAGAATATCAGCGTACGCTGCGCGAGGCATGGTAATGGGAAAAAATGATTTCAATCAGATTGCTGACACCGTGCGCGTTGTTCATGCCGGAGGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGc                                                                                                           <  2:89832/150‑1 (MQ=25)
                                            tGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGc                                                                                                           <  2:41416/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:1716726/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)
                                                          aGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAgg                                                                                             >  1:1522723/1‑150 (MQ=255)
                                                           gCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGCCAACGTGCGAGTTGTTCATGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGGTTATTGGTAGGc                                                                                            >  2:1073696/1‑150 (MQ=255)
                                                           gCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGc                                                                                            <  1:1609132/150‑1 (MQ=255)
                                                            cGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACACCGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGcc                                                                                           >  2:1217199/1‑150 (MQ=255)
                                                             gTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCt                                                                                          <  2:1312058/150‑1 (MQ=255)
                                                              tAAGATGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTg                                                                                         <  2:589524/150‑1 (MQ=255)
                                                                    gcgcGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAagc                                                                                   <  1:409442/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)
                                                                         aGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAagcgcagc                                                                              >  1:445576/1‑150 (MQ=255)
                                                                          ggtattgTATTGGGAAAAAATGTTGTCAATCAGATTGCTTACAATGTGCGCGTGGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAGAAAGTCGTTAAATTTCAATGTATTGCTGTTTATTTGGAGGCCTGATAagcgcagcg                                                                             <  2:1141393/144‑1 (MQ=21)
                                                                          ggAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAagcgcagcg                                                                             >  1:1500054/1‑150 (MQ=255)
                                                                            aaTGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAgcgcagcgca                                                                           <  1:891659/150‑1 (MQ=255)
                                                                                  aaTGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGc                                                                     <  1:1019717/150‑1 (MQ=255)
                                                                                       gAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCt                                                                >  1:470298/1‑150 (MQ=255)
                                                                                          aaaaTGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTAc                                                             <  2:557738/150‑1 (MQ=255)
                                                                                            aaTGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGc                                                           <  2:902317/150‑1 (MQ=255)
                                                                                                 tgTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGtc                                                      >  2:1006812/1‑150 (MQ=255)
                                                                                                  gTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGtca                                                     >  2:316621/1‑150 (MQ=255)
                                                                                                     aaTCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGtcatca                                                  >  1:601338/1‑150 (MQ=255)
                                                                                                      aTCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAg                                                 <  2:1216204/150‑1 (MQ=255)
                                                                                                      aTCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAg                                                 <  1:267403/150‑1 (MQ=255)
                                                                                                        cAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGtc                                               >  2:1498296/1‑150 (MQ=255)
                                                                                                        cAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCAGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGtc                                               >  1:629271/1‑150 (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)
                                                                                                                  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)
                                                                                                                          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)
                                                                                                                                   cATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTAtt                            <  2:360535/142‑1 (MQ=255)
                                                                                                                                   cATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATgtggt                    <  1:1601800/150‑1 (MQ=255)
                                                                                                                                    aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATgtggtg                   >  2:1163331/3‑150 (MQ=38)
                                                                                                                                    aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTAGTGCAGGTTAGTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCTTCAGTCTCATGCCCATAACCTTATTGCTgtggtg                   >  1:1057357/1‑150 (MQ=25)
                                                                                                                                       ccGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCCGTTTATTGGTAGGCCTTATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATgtggtggtg                <  2:554394/150‑1 (MQ=255)
                                                                                                                                       ccGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATgtggtggtg                <  2:1619837/150‑1 (MQ=255)
                                                                                                                                         ggATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCTCATAACCTTATTGATGTGGTGGTggg              <  1:29567/150‑1 (MQ=255)
                                                                                                                                          gATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGt             <  1:1114079/150‑1 (MQ=255)
                                                                                                                                              cggcgTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAAt         >  1:38701/1‑150 (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:193644/1‑150 (MQ=255)
                                                                                                                                                   tGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTGAgg    <  2:794224/150‑1 (MQ=255)
                                                                                                                                                   tGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTGAgg    >  2:1546200/1‑150 (MQ=255)
                                                                                                                                                     aaCGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATGGATGTGGTGGTGGGTAAATTGAGGcc  >  2:1273163/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCAAGAACCAGTGTTCCTGGCATCTGGCAAGGGTGTGAGGAAGATGGTGTCGAATATCAGCGTAAGCTGCGCGAGGAATGGTAATGGGAAAAAATGATGTCAATCAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCAATCTTACGCTTGTCATCAGTCTCATGCCCATAACCTTATTGATGTGGTGGTGGGTAAATTGAGGCC  >  NC_000913/4079886‑4080184

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: