Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,252,599 A→G 6.2% S512G (AGT→GGT)  dhaR → DNA‑binding transcriptional dual regulator DhaR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,252,5990AG 6.2% 192.1 / 2.9 64S512G (AGT→GGT) dhaRDNA‑binding transcriptional dual regulator DhaR
Reads supporting (aligned to +/- strand):  ref base A (31/29);  new base G (2/2);  total (33/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 = 7.77e-01

AACAACGACCGCGGACCTCGCAATGCTGGTGGAACAAAATCGTTTTAGTCGCCAGCTGTATTACGCGCTGCATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCG  >  NC_000913/1252451‑1252746
                                                                                                                                                    |                                                                                                                                                   
aacaacGACCGCGGACCTCGCAATGCTGGTGGAACAAAATCGTTTTAGTCGCCAGCTGTATTACGCGCTGCATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAg                                                                                                                                                    <  2:1226933/150‑1 (MQ=255)
  caacGACCGCGGACCTCGCAATGCTGGTGGAACAAAATCGTTTTAGTCGCCAGCTGTATTACGCGCTGCATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTc                                                                                                                                                  <  1:109452/150‑1 (MQ=255)
  caacGACCGCGGACCTCGCAATGCTGGTGGAACAAAATCGTTTTAGTCGCCAGCTGTATTACGCGCTGCATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTc                                                                                                                                                  >  1:740055/1‑150 (MQ=255)
   aacGACCGCGGACCTCGCAATGCTGGTGGAACAAAATCGTTTTAGTCGCCAGCTGTATTACGCGCTGCATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCt                                                                                                                                                 >  1:604132/1‑150 (MQ=255)
     cGACCGCGGACCTCGCAATGCTGGTGGAACAAAATCGTTTTAGTCGCCAGCTGTATTACGCGCTGCATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTg                                                                                                                                               <  2:1445905/150‑1 (MQ=255)
        ccGCGGACCTCGCAATGCTGGTGGAACAAAATCGTTTTAGTCGCCAGCTGTATTACGCGCTGCATGCATCTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGaaa                                                                                                                                            >  2:385546/1‑150 (MQ=255)
                 tCGCAATGCTGGTGGAACAAAATCGTTTTAGTCGCCAGCTGTATTACGCGCTGGATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTtct                                                                                                                                   <  2:929069/150‑1 (MQ=255)
                   gCAATGCTGGTGGAACAAAATCGTTTTAGTCGCCAGCTGTATTACGCGCTGCATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTctct                                                                                                                                 <  1:1481654/150‑1 (MQ=255)
                    cAATGCTGGTGGAACAAAATCGTTTTAGTCGCCAGCTGTATTACGCGCTGCATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTa                                                                                                                                >  2:321851/1‑150 (MQ=255)
                    cAATGCTGGTGGAACAAAATCGTTTTAGTCGCCAGCTGTATTACGCGCTGCATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTa                                                                                                                                >  2:125054/1‑150 (MQ=255)
                           ggtggAACAAAATCGTTTTAGTCGCCAGCTGTATTACGCGCTGCATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCt                                                                                                                         <  1:591293/150‑1 (MQ=255)
                           ggtggAACAAAATCGTTTTAGTCGCCAGCTGTATTACGCGCTGCATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCt                                                                                                                         >  2:1221442/1‑150 (MQ=255)
                           ggtgcAACAAAATCGTTTTAGTCGCCAGCTCTATTACGCGCTCCATCCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCACCAATCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACCCTACTCTACGCGCCt                                                                                                                         >  2:1304053/1‑150 (MQ=255)
                             tggAACAAAATCGTTTTAGTCGCCAGCTGTATTACGCGCTGCATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGa                                                                                                                       >  1:443149/1‑150 (MQ=255)
                              ggAACAAAATCGTTTTAGTCGCCAGCTGTATTACGCGCTGCATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGaa                                                                                                                      >  1:873224/1‑150 (MQ=255)
                                 aCAAAATCGTTTTAGTCGCCAGCTGTATTACGCGCTGCATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAAt                                                                                                                   >  1:1133583/1‑150 (MQ=255)
                                     aaTCGTTTTAGTCGCCAGCTGTATTACGCGCTGCATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATtgat                                                                                                               <  1:1345740/150‑1 (MQ=255)
                                         gTTTTAGTCGCCAGCTGTATTACGCGCTGCATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACg                                                                                                           <  2:873224/150‑1 (MQ=255)
                                         gTTTTAGTCGCCAGCTGTATTACGCGCTGCATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCAGCTGAAAATTGATGACg                                                                                                           >  1:552944/1‑150 (MQ=255)
                                          ttttAGTCGCCAGCTGTATTACGCGCTGCATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGa                                                                                                          <  1:321851/150‑1 (MQ=255)
                                            ttAGTCGCCAGCTGTATTACGCGCTGCATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATg                                                                                                        >  2:558203/1‑150 (MQ=255)
                                            ttAGTCGCCAGCTGTATTACGCGCTGCATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATg                                                                                                        <  1:856371/150‑1 (MQ=255)
                                                 cGCCAGCTGTATTACGCGCTGCGTGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCctc                                                                                                   >  1:1304765/1‑150 (MQ=255)
                                                 cGCCAGCTGTATTACGCGCTGCATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCGGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATtg                                                                                                                 >  1:837158/1‑136 (MQ=255)
                                                 cGCCAGCTGTATTACGCGCTGCATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCGGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATtg                                                                                                                 >  1:842603/1‑136 (MQ=255)
                                                 cGCCAGCTGTATTACGCGCTGCATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCGGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATtg                                                                                                                 <  2:837158/136‑1 (MQ=255)
                                                 cGCCAGCTGTATTACGCGCTGCATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCGGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATtg                                                                                                                 <  2:842603/136‑1 (MQ=255)
                                                         gTATTACGCGCTGCATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCt                                                                                           >  1:161324/1‑150 (MQ=255)
                                                         gTATTACGCGCTGCATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCt                                                                                           >  2:608007/1‑150 (MQ=255)
                                                          tATTACGCGCTGCATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTg                                                                                          <  2:740055/150‑1 (MQ=255)
                                                          tATTACGCGCTGCATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTg                                                                                          <  1:1343458/150‑1 (MQ=255)
                                                           aTTACGCGCTGCATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTgg                                                                                         <  2:604132/150‑1 (MQ=255)
                                                                 cgcTGCATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCtt                                                                                   >  1:240651/1‑150 (MQ=255)
                                                                  gcTGCATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTACGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTtg                                                                                  <  1:411578/150‑1 (MQ=255)
                                                                   cTGCATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTtgt                                                                                 >  2:1378033/1‑150 (MQ=255)
                                                                       atgcatTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCAt                                                                             >  2:860333/1‑150 (MQ=255)
                                                                          catTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGc                                                                          >  1:171261/1‑150 (MQ=255)
                                                                            tttGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTAGGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCa                                                                        <  1:1304053/150‑1 (MQ=255)
                                                                              tGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAgg                                                                      >  2:896533/1‑150 (MQ=255)
                                                                               gAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGc                                                                     >  1:199328/1‑150 (MQ=255)
                                                                                     aCCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGAt                                                               >  2:27507/1‑150 (MQ=255)
                                                                                      ccATCCCGCCTCTGCGTATGCGGCTTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGAtt                                                              <  1:1490281/150‑1 (MQ=255)
                                                                                         tCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTg                                                           <  1:1406054/150‑1 (MQ=255)
                                                                                          cccGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGa                                                          <  1:1261946/150‑1 (MQ=255)
                                                                                          cccGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGa                                                          <  1:922295/150‑1 (MQ=255)
                                                                                            cGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTGTTGAAACACGCTTCTCTCCGCGGGTGAATATTGATGCCGATGCCCTCGCTCGCCTGGGTTCTTGTGCATGTCCAGGCACCGATTTTGAAc                                                        >  1:1239191/1‑150 (MQ=255)
                                                                                                       aTGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTc                                             >  1:229191/1‑150 (MQ=255)
                                                                                                          cggcgTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATc                                          >  2:161859/1‑150 (MQ=255)
                                                                                                           ggcgTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCg                                         <  2:443149/150‑1 (MQ=255)
                                                                                                                  cagcaTTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCt                                  <  1:162241/150‑1 (MQ=255)
                                                                                                                       ttCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGctc                             >  1:674218/1‑150 (MQ=255)
                                                                                                                           ggCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGag                         >  2:1085447/1‑150 (MQ=255)
                                                                                                                              gcTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGagtag                      <  1:385546/150‑1 (MQ=255)
                                                                                                                              gcTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGagtag                      <  1:1078788/150‑1 (MQ=255)
                                                                                                                                 ggTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGGATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGa                   <  2:253942/150‑1 (MQ=255)
                                                                                                                                  gTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGAt                  <  2:161324/150‑1 (MQ=255)
                                                                                                                                   tGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATa                 >  1:1252397/1‑150 (MQ=255)
                                                                                                                                      aTAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACg              <  1:608007/150‑1 (MQ=255)
                                                                                                                                        aaCAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGCGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACggg            <  2:983551/150‑1 (MQ=255)
                                                                                                                                        aaCAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCAGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACggg            <  2:552944/150‑1 (MQ=255)
                                                                                                                                         aCAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGgc           >  2:1086998/1‑150 (MQ=255)
                                                                                                                                         aCAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGgc           >  1:780166/1‑150 (MQ=255)
                                                                                                                                            aaTTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCa        >  1:487085/1‑150 (MQ=255)
                                                                                                                                                  gCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTcgcg  <  1:661409/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AACAACGACCGCGGACCTCGCAATGCTGGTGGAACAAAATCGTTTTAGTCGCCAGCTGTATTACGCGCTGCATGCATTTGAAATTACCATCCCGCCTCTGCGTATGCGGCGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCG  >  NC_000913/1252451‑1252746

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