Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 3,440,086 C→A G315V (GGC→GTC)  rpoA ← RNA polymerase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi3,440,0860CA98.5% 230.8 / NA 68G315V (GGC→GTC) rpoARNA polymerase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (37/30);  total (37/31)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGC  >  NC_000913_3_bme_pgi/3439938‑3440231
                                                                                                                                                    |                                                                                                                                                 
aaCATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGacc                                                                                                                                                 <  1:286419/151‑1 (MQ=255)
 aCATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCa                                                                                                                                                <  1:26533/151‑1 (MQ=255)
          ggCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGt                                                                                                                                        >  3:416582/1‑150 (MQ=255)
               tGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACgg                                                                                                                                  >  2:457606/1‑151 (MQ=255)
                   tgctgTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGaa                                                                                                                               <  3:323562/150‑1 (MQ=255)
                       gTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAg                                                                                                                          >  2:4159/1‑151 (MQ=255)
                              ttCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGATAGTCCACGGGAAGCCAGCACGTc                                                                                                                    >  4:267989/1‑150 (MQ=255)
                                  gTTGACGACCACTCTTACGATGGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTa                                                                                                                >  3:427737/1‑150 (MQ=255)
                                    tGACGACCACTCTTACGATGGCGCATGACCTTATCCTTTTCAGTAAAACCTTAACCTGTGCTCCGGTTACTCGTCAGCGATGCTTGCCGGCGGCCAGTTTTCCAGGCGCATGACCAGTGACAGCCCACGGGACGCCAGCACGTCTTTAAt                                                                                                              <  4:18620/150‑1 (MQ=255)
                                     gacgacCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAAtct                                                                                                            >  1:311580/1‑151 (MQ=255)
                                      acgacCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAAtct                                                                                                            <  4:416582/150‑1 (MQ=255)
                                       cgacCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCa                                                                                                          >  1:209678/1‑151 (MQ=255)
                                             ctctTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGAc                                                                                                                                           >  1:35815/1‑112 (MQ=255)
                                             ctctTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGAc                                                                                                                                           <  2:35815/112‑1 (MQ=255)
                                              tctTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAgaga                                                                                                                                            <  4:242243/110‑1 (MQ=255)
                                              tctTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAgaga                                                                                                                                            >  3:242243/1‑110 (MQ=255)
                                              tctTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAagag                                                                                                   <  1:234621/151‑1 (MQ=255)
                                                ttACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAgaga                                                                                                  >  3:330935/1‑150 (MQ=255)
                                                   cGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGAtttt                                                                                              >  2:638427/1‑151 (MQ=255)
                                                     aTGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGAtttttt                                                                                            >  1:70471/1‑151 (MQ=255)
                                                      tGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGAtttttt                                                                                            >  3:452120/1‑150 (MQ=255)
                                                       ggCGCCTGACCGTATCCTTCTCAGTAAAACCTTAGCCTGTGATCCGGTTACTCGTCAGTGATGCTTGCTGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTAc                                                                                          <  2:77364/151‑1 (MQ=255)
                                                          gcATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCaa                                                                                       >  2:953006/1‑151 (MQ=255)
                                                                 cTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGAc                                                                                                                                           <  2:345164/92‑1 (MQ=255)
                                                                 cTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGAc                                                                                                                                           >  1:345164/1‑92 (MQ=255)
                                                                     tCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAagag                                                                                                   <  3:96012/128‑1 (MQ=255)
                                                                     tCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAagag                                                                                                   >  4:96012/1‑128 (MQ=255)
                                                                       cTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGttt                                                                           >  4:114165/1‑150 (MQ=255)
                                                                       cTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACAGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGttt                                                                           <  4:383235/150‑1 (MQ=255)
                                                                        ttCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTa                                                                         <  1:457606/151‑1 (MQ=255)
                                                                           tcAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAg                                                                       <  4:237506/150‑1 (MQ=255)
                                                                             aGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAg                                                                    >  1:802159/1‑151 (MQ=255)
                                                                              gTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGacc                                                                                                                                                 >  1:869409/1‑73 (MQ=255)
                                                                              gTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGacc                                                                                                                                                 <  2:869409/73‑1 (MQ=255)
                                                                                 aaaCCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGCGCTCa                                                                >  2:647748/1‑151 (MQ=255)
                                                                                          ccTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTa                                                       <  2:575537/151‑1 (MQ=255)
                                                                                           cTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTa                                                       >  3:272989/1‑150 (MQ=255)
                                                                                             gtgATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGCGAAGCCAGCACGTCTGTAATCTCAGTAGGAGATTTTTTACCACGGTTAAGCGTTTTAACGAGCTCAATGTCAGTACg                                                     <  4:225760/150‑1 (MQ=255)
                                                                                              ttaTCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTGTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAACTCAACCTCGGTACGCt                                                   <  2:499494/149‑1 (MQ=255)
                                                                                              tgATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCt                                                   >  2:854900/1‑151 (MQ=255)
                                                                                                aTCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTg                                                  >  3:180886/1‑150 (MQ=255)
                                                                                                 tCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTa                                                <  1:946324/151‑1 (MQ=255)
                                                                                                  ccGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTa                                                >  4:282219/1‑150 (MQ=255)
                                                                                                    ggTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCa                                             <  1:638427/151‑1 (MQ=255)
                                                                                                        aCTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGc                                                                                                                             >  3:231467/1‑67 (MQ=255)
                                                                                                        aCTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGc                                                                                                                             <  4:231467/67‑1 (MQ=255)
                                                                                                         cTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCa                                        <  1:191202/151‑1 (MQ=255)
                                                                                                         cTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCa                                        >  1:95427/1‑151 (MQ=255)
                                                                                                           cGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAc                                       >  3:566820/1‑150 (MQ=255)
                                                                                                              cAGCGATGCTTGCCGGTAGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGat                                   >  1:836948/1‑151 (MQ=255)
                                                                                                                   aTGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGt                              <  2:633441/151‑1 (MQ=255)
                                                                                                                    tGCTTTCCGGTGGCCAGTTTTGCCGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTg                             <  1:904312/151‑1 (MQ=255)
                                                                                                                    tGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGt                              <  4:452120/150‑1 (MQ=255)
                                                                                                                      cTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAg                                                                       >  1:605709/1‑107 (MQ=255)
                                                                                                                      cTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAg                                                                       <  2:605709/107‑1 (MQ=255)
                                                                                                                       ttGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGAt                          <  1:647748/151‑1 (MQ=255)
                                                                                                                          ccGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGc                       <  1:842937/151‑1 (MQ=255)
                                                                                                                                gCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGc                 >  2:820917/1‑151 (MQ=255)
                                                                                                                                 ccAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCt                <  2:628162/151‑1 (MQ=255)
                                                                                                                                  cAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCt                >  4:831/1‑150 (MQ=255)
                                                                                                                                   aGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCttt              >  1:103979/1‑151 (MQ=255)
                                                                                                                                    gTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGAAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTa             >  2:188724/1‑151 (MQ=255)
                                                                                                                                      tttCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAg           <  1:4159/151‑1 (MQ=255)
                                                                                                                                       ttCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAgg          <  2:311580/151‑1 (MQ=255)
                                                                                                                                         ccAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGc         >  4:559564/1‑150 (MQ=255)
                                                                                                                                          cAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATATCTTCTGCTTTAAGGCAg       <  2:488365/151‑1 (MQ=255)
                                                                                                                                          cAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAg       >  1:205299/1‑151 (MQ=255)
                                                                                                                                           aGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCg                                                                              <  2:113406/79‑1 (MQ=255)
                                                                                                                                           aGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCg                                                                              >  1:113406/1‑79 (MQ=255)
                                                                                                                                            ggCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGtt     >  1:651070/1‑151 (MQ=255)
                                                                                                                                               gcATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGc  >  1:320061/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGC  >  NC_000913_3_bme_pgi/3439938‑3440231

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: