Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,559,217 A→C D317A (GAT→GCT)  rtcR → sigma 54‑dependent transcriptional regulator of rtcBA expression

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,559,2170AC96.7% 199.1 / ‑0.9 61D317A (GAT→GCT) rtcRsigma 54‑dependent transcriptional regulator of rtcBA expression
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (33/26);  total (34/27)
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.73e-01

CCGGGGCGCGGGAATCTCGTGAAGGTTTATTACGCAGCGCCAACGGCGGAATGTTGTTTCTTGATGAGATTGGCGAACTGGGCGCAGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGATCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGAAGATCTGTACGCGCGGATCAATCTCTGGACCTTCACCCTGCCGGGTCTACGCCAGCGCCAG  >  NC_000913/3559070‑3559362
                                                                                                                                                   |                                                                                                                                                 
ccGGGGCGCGGGAATCTCGTGAAGGTTTATTACGCAGCGCCAACGGCGGAATGTTGTTTCTTGATGAGATTGGCGAACTGGGCGCAGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCg                                                                                                                                                <  1:418760/151‑1 (MQ=255)
  ggggCGCGGGAATCTCGTGAAGGTTTATTACGCAGCGCCAACGGCGGAATGTTGTTTCTTGATGAGATTGGCGAACTGGGCGCAGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGcc                                                                                                                                              <  1:930701/151‑1 (MQ=255)
          ggAATCTCGTGAAGGTTTATTACGCAGCGCCAACGGCGGAATGTTGTTTCTTGATGAGATTGGCGAACTGGGCGCAGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGc                                                                                                                                      <  2:364461/151‑1 (MQ=255)
           gAATCTCGTGAAGGTTTATTACGCAGCGCCAACGGCGGAATGTTGTTTCTTGATGAGATTGGCGAACTGGGCGCAGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCa                                                                                                                                     <  2:388470/151‑1 (MQ=255)
             aTCTCGTGAAGGTTTATTACGCAGCGCCAACGGCGGAATGTTGTTTCTTGATGAGATTGGCGAACTGGGCGCAGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGt                                                                                                                                   >  2:801831/1‑151 (MQ=255)
               ctcGTGAAGGTTTATTACGCAGCGCCAACGGCGGAATGTTGTTTCTTGATGAGATTGGCGAACTGGGCGCAGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTg                                                                                                                                  >  2:591014/1‑150 (MQ=255)
                  gTGAAGGTTTATTACGCAGCGCCAACGGCGGAATGTTGTTTCTTGATGAGATTGGCGAACTGGGCGCAGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGAttt                                                                                                                              >  2:121980/1‑151 (MQ=255)
                                   aGCGCCAACGGCGGAATGTTGTTTCTTGATGAGATTGGCGAACTGGGCGCAGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGa                                                                                                              >  1:637663/1‑150 (MQ=255)
                                      gcCAACGGCGGAATGTTGTTTCTTGATGAGATTGGCGAACTGGGCGCAGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACgg                                                                                                          <  2:624625/151‑1 (MQ=255)
                                        cAACGGCGGAATGTTGTTTCTTGATGAGATTGGCGAACTGGGCGCAGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTg                                                                                                        >  2:14925/1‑151 (MQ=255)
                                         aaCGGCGGAATGTTGTTTCTTGATGAGATTGGCGAACTGGGCGCAGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTgc                                                                                                       <  2:4740/151‑1 (MQ=255)
                                         aaCGGCGGAATGTTGTTTCTTGATGAGATTGGCGAACTGGGCGCAGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTg                                                                                                        >  2:116983/1‑150 (MQ=255)
                                           cggcggAATGTTGTTTCTTGATGAGATTGGCGAACTGGGCGCAGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTgcgc                                                                                                     >  1:523188/1‑151 (MQ=255)
                                                gAATGTTGTTTCTTGATGAGATTGGCGAACTGGGCGCAGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGAttt                                                                                                >  1:751681/1‑151 (MQ=255)
                                                        tttCTTGATGAGATTGGCGAACTGGGCGCAGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGc                                                                                        >  2:164450/1‑151 (MQ=255)
                                                         ttCTTGATGAGATTGGCGAACTGGGCGCAGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCt                                                                                       >  1:575367/1‑151 (MQ=255)
                                                          tCTTGATGAGATTGGCGAACTGGGCGCAGACGAACAGTCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTg                                                                                      >  2:323268/1‑151 (MQ=255)
                                                            ttGATGAGATTGGCGAACTGGGCGCAGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGt                                                                                    >  1:285706/1‑151 (MQ=255)
                                                                 gagaTTGGCGAACTGGGCGCAGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCg                                                                               <  1:557304/151‑1 (MQ=255)
                                                                  agaTTGGCGAACTGGGCGCAGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCg                                                                               >  2:740373/1‑150 (MQ=255)
                                                                    aTTGGCGAACTGGGCGCAGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAg                                                                            <  2:644280/151‑1 (MQ=255)
                                                                      tGGCGAACTGGGCGCAGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAgg                                                                           <  1:383343/150‑1 (MQ=255)
                                                                       ggCGAACTGGGCGCAGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGc                                                                          >  1:863637/1‑150 (MQ=255)
                                                                        gCGAACTGGGCGCAGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCaa                                                                        <  2:318795/151‑1 (MQ=255)
                                                                          gAACTGGGCGCAGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAAt                                                                      >  2:252664/1‑151 (MQ=255)
                                                                          gAACTGGGCGCAGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAAt                                                                      >  2:242124/1‑151 (MQ=255)
                                                                                 gcgcAGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGa                                                               <  1:395242/151‑1 (MQ=255)
                                                                                 gcgcAGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGa                                                               <  1:144880/151‑1 (MQ=255)
                                                                                  cgcAGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGaa                                                              <  2:538298/151‑1 (MQ=255)
                                                                                  cgcAGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGa                                                               >  1:158288/1‑150 (MQ=255)
                                                                                  cgcAGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGa                                                               >  2:381572/1‑150 (MQ=255)
                                                                                    cAGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGa                                                               >  1:496684/1‑148 (MQ=255)
                                                                                     aGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGAAGAt                                                           <  1:526993/151‑1 (MQ=255)
                                                                                         gAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGATCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGAAGATCTGt                                                       <  1:247471/151‑1 (MQ=255)
                                                                                          aaCAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGAAGATCTGTa                                                      <  2:448042/151‑1 (MQ=255)
                                                                                           aCAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGAAGATCTGTAc                                                     >  2:570435/1‑151 (MQ=255)
                                                                                             aGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGAAGATCTGTAcgc                                                   <  2:855935/151‑1 (MQ=255)
                                                                                                             aGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGAAGATCTGTACGCGCGGATCAATCTCTgg                                   <  2:882925/151‑1 (MQ=255)
                                                                                                               ccATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGAAGATCTGTACGCGCGGATCAATCTCTGGAc                                 <  1:481440/151‑1 (MQ=255)
                                                                                                                  ttGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGAAGATCTGTACGCGCGGATCAATCTCTGGACCtt                              >  2:235293/1‑151 (MQ=255)
                                                                                                                   tGAAGAGAAAACCTTTTACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGAAGATCTGTACGCGCGGATCAATCTCTGGACCTTc                             >  2:287492/1‑151 (MQ=255)
                                                                                                                                 tttACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGAAGATCTGTACGCGCGGATCAATCTCTGGACCTTCACCCTGCCGGGTCt               >  2:149101/1‑151 (MQ=255)
                                                                                                                                  ttACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGAAGATCTGTACGCGCGGATCAATCTCTGGACCTTCACCCTGCCGGGTc                >  2:425239/1‑149 (MQ=255)
                                                                                                                                  ttACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGAAGATCTGTACGCGCGGATCAATCTCTGGACCTTCACCCTGCCGGGTCt               <  2:343191/150‑1 (MQ=255)
                                                                                                                                  ttACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGAAGATCTGTACGCGCGGATCAATCTCTGGACCTTCACCCTGCCGGGTCt               <  1:753495/150‑1 (MQ=255)
                                                                                                                                   tACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGAAGATCTGTGCGCGCGGATCAATCTCTGGACCTTCGCCCTGCCGGGTCt               <  1:77088/149‑1 (MQ=255)
                                                                                                                                   tACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGAAGATCTGTACGCGCGGATCAATCTCTGGACCTTCACCCTGCCGGGTCTa              <  2:67312/150‑1 (MQ=255)
                                                                                                                                   tACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGAAGATCTGTACGCGCGGATCAATCTCTGGACCTTCACCCTGCCGGGTCTAc             >  1:641255/1‑151 (MQ=255)
                                                                                                                                   tACCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGAAGATCTGTACGCGCGGATCAATCTCTGGACCTTCACCCTGCCGGGTCTAc             >  2:459163/1‑151 (MQ=255)
                                                                                                                                    aCCCGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGAAGATCTGTACGCGCGGATCAATCTCTGGACCTTCACCCTGCCGGGTCTAcg            >  1:915583/1‑151 (MQ=255)
                                                                                                                                     cccGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGAAGATCTGTACGCGCGGATCAATCTCTGGACCTTCACCCTGCCGGGTCTAcgc           >  2:250775/1‑151 (MQ=255)
                                                                                                                                      ccGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGAAGATCTGTACGCGCGGATCAATCTCTGGACCTTCACCCTGCCGGGTCTAcgcc          <  2:98850/151‑1 (MQ=255)
                                                                                                                                      ccGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGAAGATCTGTACGCGCGGATCAATCTCTGGACCTTCACCCTGCCGGGTCTAcgcc          >  1:175451/1‑151 (MQ=255)
                                                                                                                                      ccGTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGAAGATCTGTACGCGCGGATCAATCTCTGGACCTTCACCCTGCCGGGTCTAcgcc          >  1:84173/1‑151 (MQ=255)
                                                                                                                                        gTTTGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGAAGATCTGTACGCGCGGATCAATCTCTGGACCTTCACCCTGCCGGGTCTAcgcc          >  1:655960/1‑149 (MQ=255)
                                                                                                                                        gTTTGGCAGCGATCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGAAGATCTGTACGCGCGGATCAATCTCTGGACCTTCACCCTGCCGGGTCTAcgcc          >  1:449944/1‑149 (MQ=255)
                                                                                                                                         tttGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGAAGATCTGTACGCGCGGATCAATCTCTGGACCTTCACCCTGCCGGGTCTAcgccagc       >  1:626231/1‑151 (MQ=255)
                                                                                                                                         tttGGCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGAAGATCTGTACGCGCGGATCAATCTCTGGACCTTCACCCTGCCGGGTCTAcgccagc       >  2:368340/1‑151 (MQ=255)
                                                                                                                                            ggCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGAAGATCTGTACGCGCGGATCAATCTCTGGACCTTCACCCTGCCGGGTCTAcgccagcgcc    <  2:401094/151‑1 (MQ=255)
                                                                                                                                             gCAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGAAGATCTGTACGCGCGGATCAATCTCTGGACCTTCACCCTGCCGGGTCTAcgccagcgcca   <  2:9838/151‑1 (MQ=255)
                                                                                                                                              cAGCGCTCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGAAGATCTGTACGCGCGGATCAATCTCTGGACCTTCACCCTGCCGGGTCTAcgccagcgccag  <  1:557316/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCGGGGCGCGGGAATCTCGTGAAGGTTTATTACGCAGCGCCAACGGCGGAATGTTGTTTCTTGATGAGATTGGCGAACTGGGCGCAGACGAACAGGCAATGCTGCTGAAAGCCATTGAAGAGAAAACCTTTTACCCGTTTGGCAGCGATCGCCAGGTGAGCAGTGATTTTCAGCTTATCGCCGGAACGGTGCGCGATTTGCGCCAGCTGGTTGCCGAAGGCAAATTTCGCGAAGATCTGTACGCGCGGATCAATCTCTGGACCTTCACCCTGCCGGGTCTACGCCAGCGCCAG  >  NC_000913/3559070‑3559362

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