Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,780,450 A→C S208S (TCA→TCC)  lldD → L‑lactate dehydrogenase, FMN‑linked

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,780,4500AC98.4% 224.1 / NA 62S208S (TCA→TCC) lldDL‑lactate dehydrogenase, FMN‑linked
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (31/30);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CACCGGGCGCACGCTACCGTGATGCGCATTCAGGTATGAGCGGCCCGAACGCGGCAATGCGCCGCTACTTGCAAGCGGTGACACATCCGCAATGGGCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCAGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAATGGATCCGCGATTTCTGGGATGGCCCGATGGTGATCAAAGGGATCCTCGATCCGGA  >  NC_000913/3780305‑3780596
                                                                                                                                                 |                                                                                                                                                  
cacCGGGCGCACGCTACCGTGATGCGCATTCAGGTATGAGCGGCCCGAACGCGGCAATGCGCCGCTACTTGCAAGCGGTGACACATCCGCAATGGGCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCtt                                                                                                                                                >  2:268437/1‑150 (MQ=255)
    gggCGCACGCTACCGTGATGCGCATTCAGGTATGAGCGGCCCGAACGCGGCAATGCGCCGCTACTTGCAAGCGGTGACACATCCGCAATGGGCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTAtctc                                                                                                                                           <  1:135777/151‑1 (MQ=255)
          aCGCTACCGTGATGCGCATTCAGGTATGAGCGGCCCGAACGCGGCAATGCGCCGCTACTTGCAAGCGGTGACACATCCGCAATGGGCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCaa                                                                                                                                      >  2:634104/1‑150 (MQ=255)
            gCTACCGTGATGCGCATTCAGGTATGAGCGGCCCGAACGCGGCAATGCGCCGCTACTTGCAAGCGGTGACACATCCGCAATGGGCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAacc                                                                                                                                   <  2:406843/151‑1 (MQ=255)
              tACCGTGATGCGCATTCAGGTATGAGCGGCCCGAACGCGGCAATGCGCCGCTACTTGCAAGCGGTGACACATCCGCAATGGGCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAaccg                                                                                                                                  >  1:742956/1‑150 (MQ=255)
               aCCGTGATGCGCATTCAGGTATGAGCGGCCCGAACGCGGCAATGCGCCGCTACTTGCAAGCGGTGACACATCCGCAATGGGCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAaccg                                                                                                                                  <  2:536965/149‑1 (MQ=255)
                  gTGATGCGCATTCAGGTATGAGCGGCCCGAACGCGGCAATGCGCCGCTACTTGCAAGCGGTGACACATCCGCAATGGGCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCgg                                                                                                                             <  2:555061/151‑1 (MQ=255)
                   tGATGCGCATTCAGGTATGAGCGGCCCGAACGCGGCAATGCGCCGCTACTTGCAAGCGGTGACACATCCGCAATGGGCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCgg                                                                                                                             >  1:584715/1‑150 (MQ=255)
                           aTTCAGGTATGAGCGGCCCGAACGCGGCAATGCGCCGCTACTTGCAAGCGGTGACACATCCGCAATGGGCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGa                                                                                                                    <  1:841459/151‑1 (MQ=255)
                               aGGTATGAGCGGCCCGAACGCGGCAATGCGCCGCTACTTGCAAGCGGTGACACATCCGCAATGGGCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTAc                                                                                                                <  1:661177/151‑1 (MQ=255)
                                   aTGAGCGGCCCGAACGCGGCAATGCGCCGCTACTTGCAAGCGGTGACACATCCGCAATGGGCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCg                                                                                                            >  1:636981/1‑151 (MQ=255)
                                      aGCGGCCCGAACGCGGCAATGCGCCGCTACTTGCAAGCGGTGACACATCCGCAATGGGCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCggct                                                                                                         <  2:676165/151‑1 (MQ=255)
                                                 cgcgGCAATGCGCCGCTACTTGCAAGCGGTGACACATCCGCAATGGGCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAAt                                                                                              >  1:530424/1‑151 (MQ=255)
                                                   cgGCAATGCGCCGCTACTTGCAAGCGGTGACACATCCGCAATGGGCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATa                                                                                             >  1:711538/1‑150 (MQ=255)
                                                      cAATGCGCCGCTACTTGCAAGCGGTGACACATCCGCAATGGGCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACtt                                                                                         <  1:434617/151‑1 (MQ=255)
                                                       aaTGCGCCGCTACTTGCAAGCGGTGACACATCCGCAATGGGCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTc                                                                                        <  2:633668/151‑1 (MQ=255)
                                                        aTGCGCCGCTACTTGCAAGCGGTGACACATCCGCAATGGGCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCg                                                                                       >  1:591438/1‑151 (MQ=255)
                                                         tGCGCCGCTACTTGCAAGCGGGGAGACATCCGCAATGGGCGGGGGATGTGGGCCTGAACGGTCGTCAACATGATTTAGGCAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGACGATTACATCGGCTGGCTGGGGAATAACtt                                                                                         >  2:360067/1‑148 (MQ=255)
                                                             ccgcTACTTGCAAGCGGTGACACATCCGCAATGGGCGTGGGATGGGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTGATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATc                                                                                    >  2:369515/1‑149 (MQ=255)
                                                               gcTACTTGCAAGCGGTGACACATCCGCAATGGGCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTc                                                                                <  1:9434/151‑1 (MQ=255)
                                                               gcTACTTGCAAGCGGTGACACATCCGCAATGGGCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTc                                                                                <  2:334726/151‑1 (MQ=255)
                                                                 tACTTGCAAGCGGTGACACATCCGCAATGGGCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTcc                                                                               >  1:624343/1‑150 (MQ=255)
                                                                        aaGCGGTGACACATCCGCAATGGGCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATg                                                                       <  2:712066/151‑1 (MQ=255)
                                                                             gTGACACATCCGCAATGGGCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAg                                                                  <  2:683611/151‑1 (MQ=255)
                                                                              tGACACATCCGCAATGGGCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGaa                                                                    >  2:555130/1‑148 (MQ=255)
                                                                                  acaTCCGCAATGGGCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCtt                                                             <  1:104247/151‑1 (MQ=255)
                                                                                   cagccGCAATGGGCGTGAGATGTGGGCCTGAGCGGTCGTCCACATGATTTAGGTAATATCTCAGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTACGATCCGTCCATCTCATGGTAAGACCTTg                                                            <  2:489888/148‑1 (MQ=255)
                                                                                     tCCGCAATGGGCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGa                                                           >  2:613955/1‑150 (MQ=255)
                                                                                     tCCGCAATGGGCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGa                                                           >  2:69863/1‑150 (MQ=255)
                                                                                      ccGCAATGGGCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAAt                                                         <  2:533871/151‑1 (MQ=255)
                                                                                      ccGCAATGGGCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAAt                                                         >  2:151006/1‑151 (MQ=255)
                                                                                        gCAATGGGCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAATgg                                                       >  1:696029/1‑151 (MQ=255)
                                                                                             gggCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAATGGATcc                                                   >  2:857293/1‑150 (MQ=255)
                                                                                             gggCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAATGGATCcg                                                  >  1:481932/1‑151 (MQ=255)
                                                                                               gCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAATGGATCcgcg                                                <  1:421757/151‑1 (MQ=255)
                                                                                                   gggATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAATGGATCCGCGAttt                                            <  2:501439/151‑1 (MQ=255)
                                                                                                       tgtgGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAATGGATCCGCGATTTCTgg                                        >  2:687425/1‑151 (MQ=255)
                                                                                                        gtgGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAATGGATCCGCGATTTCTggg                                       <  1:640685/151‑1 (MQ=255)
                                                                                                        gtgGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAATGGATCCGCGATTTCTgg                                        <  1:308652/150‑1 (MQ=255)
                                                                                                         tgGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAATGGATCCGCGATTTCTggg                                       >  2:676864/1‑150 (MQ=255)
                                                                                                          gggCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAATGGATCCGCGATTTCTGGGa                                      <  2:610640/150‑1 (MQ=255)
                                                                                                          gggCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAATGGATCCGCGATTTCTGGGAt                                     <  1:309822/151‑1 (MQ=255)
                                                                                                           ggCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAATGGATCCGCGATTTCTGGGATg                                    <  1:351871/151‑1 (MQ=255)
                                                                                                               tGAACGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAATGGATCCGCGATTTCTGGGATGGccc                                >  2:106196/1‑151 (MQ=255)
                                                                                                                 aaCGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAATGGATCCGCGATTTCTGGGATGGCCCg                               >  2:724247/1‑150 (MQ=255)
                                                                                                                 aaCGGTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAATGGATCCGCGATTTCTGGGATGGCCCGa                              <  2:163416/151‑1 (MQ=255)
                                                                                                                    ggTCGTCCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAATGGATCCGCGATTTCTGGGATGGCCCGATgg                           >  1:762633/1‑151 (MQ=255)
                                                                                                                     gtcgtcCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAATGGATCCGCGATTTCTGGGATGGCCCGATGGt                          >  1:539078/1‑151 (MQ=255)
                                                                                                                     gtcgtcCACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAATGGATCCGCGATTTCTGGGATGGCCCGATGGt                          <  2:412640/151‑1 (MQ=255)
                                                                                                                          ccACATGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAATGGATCCGCGATTTCTGGGATGGCCCGATGGTGATc                      >  2:155893/1‑150 (MQ=255)
                                                                                                                             caTGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAATGGATCCGCGATTTCTGGGATGGCCCGATGGTGATCAAAg                  <  2:200031/151‑1 (MQ=255)
                                                                                                                              aTGATTTAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAATGGATCCGCGATTTCTGGGATGGCCCGATGGTGATCAAAgg                 <  1:561566/151‑1 (MQ=255)
                                                                                                                                   ttAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAATGGATCCGCGATTTCTGGGATGGCCCGATGGTGATCAAAGGGATc             <  1:395101/150‑1 (MQ=255)
                                                                                                                                   ttAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAATGGATCCGCGATTTCTGGGATGGCCCGATGGTGATCAAAGGGATc             <  2:665094/150‑1 (MQ=255)
                                                                                                                                    tAGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAATGGATCCGCGATTTCTGGGATGGCCCGATGGTGATCAAAGGGATCCt           >  2:437362/1‑151 (MQ=255)
                                                                                                                                     aGGTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAATGGATCCGCGATTTCTGGGATGGCCCGATGGTGATCAAAGGGATCCt           <  2:614318/150‑1 (MQ=255)
                                                                                                                                       gTAATATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAATGGATCCGCGATTTCTGGGATGGCCCGATGGTGATCAAAGGGATCCTCg         >  2:276577/1‑150 (MQ=255)
                                                                                                                                         aaTATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAATGGATCCGCGATTTCTGGGATGGCCCGATGGTGATCAAAGGGATCCTCGAt       <  1:501365/150‑1 (MQ=255)
                                                                                                                                         aaTATCTCCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAATGGATCCGCGATTTCTGGGATGGCCCGATGGTGATCAAAGGGATCCTCGATc      >  2:174323/1‑151 (MQ=255)
                                                                                                                                             tctcCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAATGGATCCGCGATTTCTGGGATGGCCCGATGGTGATCAAAGGGATCCTCGATCCgg   >  1:874774/1‑150 (MQ=255)
                                                                                                                                              ctcCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAATGGATCCGCGATTTCTGGGATGGCCCGATGGTGATCAAAGGGATCCTCGATCCGGa  >  2:395510/1‑150 (MQ=255)
                                                                                                                                              ctcCGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAATGGATCCGCGATTTCTGGGATGGCCCGATGGTGATCAAAGGGATCCTCGATCCGGa  >  1:13356/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CACCGGGCGCACGCTACCGTGATGCGCATTCAGGTATGAGCGGCCCGAACGCGGCAATGCGCCGCTACTTGCAAGCGGTGACACATCCGCAATGGGCGTGGGATGTGGGCCTGAACGGTCGTCCACATGATTTAGGTAATATCTCAGCTTATCTCGGCAAACCGACCGGACTGGAAGATTACATCGGCTGGCTGGGGAATAACTTCGATCCGTCCATCTCATGGAAAGACCTTGAATGGATCCGCGATTTCTGGGATGGCCCGATGGTGATCAAAGGGATCCTCGATCCGGA  >  NC_000913/3780305‑3780596

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑