Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,440,212 T→C E273G (GAA→GGA)  rpoA ← RNA polymerase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,440,2120TC98.4% 251.3 / ‑4.2 62E273G (GAA→GGA) rpoARNA polymerase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (44/17);  minor base A (1/0);  total (45/17)
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

GGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCC  >  NC_000913/3440065‑3440358
                                                                                                                                                   |                                                                                                                                                  
ggCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTg                                                                                                                                                 >  2:399283/1‑151 (MQ=255)
     gTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCttt                                                                                                                                             >  1:851078/1‑150 (MQ=255)
     gTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTa                                                                                                                                            <  2:608573/151‑1 (MQ=255)
     gTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTa                                                                                                                                            <  1:839571/151‑1 (MQ=255)
        ttCCAGGCGCATGCCCAGAGACAGTCCACGCGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCCAGGTTAGGCGTTTTAAGGAGCTCAACATCGGTACGCTGTACCAGATCACTGATATAGTGGATAGCTCCTGCTTTAAgg                                                                                                                                         >  1:195675/1‑151 (MQ=255)
         tCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAgg                                                                                                                                         >  2:672389/1‑150 (MQ=255)
          ccAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGc                                                                                                                                        >  1:753960/1‑150 (MQ=255)
               cgcATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGtt                                                                                                                                    >  2:738486/1‑149 (MQ=255)
                gcATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGc                                                                                                                                 >  1:137939/1‑151 (MQ=255)
                  aTGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGc                                                                                                                                 >  2:229202/1‑149 (MQ=255)
                      ccAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGc                                                                                                                            >  1:348800/1‑150 (MQ=255)
                        agagACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCgg                                                                                                                          >  1:22713/1‑150 (MQ=255)
                         gagaCAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGAc                                                                                                                        <  2:2664/151‑1 (MQ=255)
                          agaCAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACa                                                                                                                       <  2:396107/151‑1 (MQ=255)
                              aGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTc                                                                                                                    >  1:842189/1‑150 (MQ=255)
                                 ccACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAAtt                                                                                                                >  1:647131/1‑151 (MQ=255)
                                    cGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTcc                                                                                                              >  1:818106/1‑150 (MQ=255)
                                        aaGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTcc                                                                                                              >  2:243553/1‑146 (MQ=255)
                                         aGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCg                                                                                                        <  2:156042/151‑1 (MQ=255)
                                          gCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGt                                                                                                       <  2:820163/151‑1 (MQ=255)
                                                aCGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAAc                                                                                                   >  1:416300/1‑149 (MQ=255)
                                                 cGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAgg                                                                                                >  2:752652/1‑151 (MQ=255)
                                                 cGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAgg                                                                                                <  1:856806/151‑1 (MQ=255)
                                                   tCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGgc                                                                                              >  2:169797/1‑151 (MQ=255)
                                                      ttAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGgc                                                                                              <  2:535107/148‑1 (MQ=255)
                                                               gTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAggatcg                                                                                  >  1:271434/1‑151 (MQ=255)
                                                                    agaTTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGgatcggatcg                                                                             >  1:735316/1‑151 (MQ=255)
                                                                      aTTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGa                                                                            <  1:415880/150‑1 (MQ=255)
                                                                      aTTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGa                                                                            >  1:957726/1‑150 (MQ=255)
                                                                          tttACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTATGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAActct                                                                       >  2:918218/1‑151 (MQ=255)
                                                                              ccAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGtt                                                                   >  2:249383/1‑151 (MQ=255)
                                                                                 aGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTtct                                                                >  1:59952/1‑151 (MQ=255)
                                                                                      aGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCtcttct                                                           >  2:964144/1‑151 (MQ=255)
                                                                                        gCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCttt                                                         >  2:84138/1‑151 (MQ=255)
                                                                                        gCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCttt                                                         >  2:404228/1‑151 (MQ=255)
                                                                                            tttAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACt                                                     >  2:766139/1‑151 (MQ=255)
                                                                                             ttAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACtt                                                    >  2:310342/1‑151 (MQ=255)
                                                                                              tAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACtt                                                    <  1:460434/150‑1 (MQ=255)
                                                                                               ccttcgCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTc                                                   >  2:453936/6‑150 (MQ=255)
                                                                                               aaGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTc                                                   >  1:18503/1‑150 (MQ=255)
                                                                                                aGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGGGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTc                                                   >  1:274807/1‑149 (MQ=255)
                                                                                                  gAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGc                                               <  2:551206/151‑1 (MQ=255)
                                                                                                    gCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTg                                             <  1:809008/151‑1 (MQ=255)
                                                                                                       cAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTATGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACg                                          <  1:918207/151‑1 (MQ=255)
                                                                                                       cAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACg                                          >  2:295621/1‑151 (MQ=255)
                                                                                                         aCCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGt                                         >  2:715108/1‑150 (MQ=255)
                                                                                                                tACGCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGCTTCTCTTCTTTCACTTCAGGCTGACGTACATCACg                                 >  2:168930/1‑151 (MQ=255)
                                                                                                                   gCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTaa                              <  1:656126/151‑1 (MQ=255)
                                                                                                                   gCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTa                               >  2:468682/1‑150 (MQ=255)
                                                                                                                   gCTGTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTa                               >  1:76816/1‑150 (MQ=255)
                                                                                                                      gTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTc                           <  2:704669/151‑1 (MQ=255)
                                                                                                                      gTACCAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTc                           >  2:588443/1‑151 (MQ=255)
                                                                                                                         ccAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAAc                        <  1:672248/151‑1 (MQ=255)
                                                                                                                          cAGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACg                       >  2:75023/1‑151 (MQ=255)
                                                                                                                           aGATCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACg                       >  1:331374/1‑150 (MQ=255)
                                                                                                                              tCACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGaa                     >  1:134674/1‑149 (MQ=255)
                                                                                                                               cACCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGc                  <  1:479360/151‑1 (MQ=255)
                                                                                                                                aCCGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCt                 >  2:453044/1‑151 (MQ=255)
                                                                                                                                  cGATATAGTGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTc               >  2:718526/1‑151 (MQ=255)
                                                                                                                                          tGGATAGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAgttgtt       >  2:750571/1‑151 (MQ=255)
                                                                                                                                               aGCTCCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGcc  >  2:947870/1‑151 (MQ=255)
                                                                                                                                               aGCTACTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGcc  >  2:828245/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACATCACGTAAGTCAACGAAAGCTTCCAGTTGTTCAGCC  >  NC_000913/3440065‑3440358

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: