Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,409,494 G→T 5.2% intergenic (+33/+26) yfbR → / ← yfbS dCMP phosphohydrolase/putative transporter YfbS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,409,4940GT 5.2% 326.6 / 5.6 96intergenic (+33/+26)yfbR/yfbSdCMP phosphohydrolase/putative transporter YfbS
Reads supporting (aligned to +/- strand):  ref base G (54/37);  new base T (3/2);  total (57/39)
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.46e-01

CGACTGGAAGCGACGCTTGAAGCGCGTCGCAGCCAGGAGATGGACTACTTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTA  >  NC_000913/2409345‑2409643
                                                                                                                                                     |                                                                                                                                                     
cgACTGGAAGCGACGCTTGAAGCGCGTCGCAGCCAGGAGATGGACTACTTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAg                                                                                                                                                       <  1:1768891/150‑1 (MQ=255)
cgACTGGAAGCGACGCTTGAAGCGCGTCGCAGCCAGGAGATGGACTACTTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGAGTCTGCGTCGCATCAGGCAATAAgc                                                                                                                                                      >  2:227546/1‑150 (MQ=255)
    tGGAAGCGACGCTTGAAGCGCGTCGCAGCCAGGAGATGGACTACTTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGcc                                                                                                                                                   <  1:15688/150‑1 (MQ=255)
     ggAAGCGACGCTTGAAGCGCGTCGCAGCCAGGAGATGGACTACTTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCg                                                                                                                                                  >  2:1624622/1‑150 (MQ=255)
      gAAGCGACGCTTGAAGCGCGTCGCAGCCAGGAGATGGACTACTTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCgg                                                                                                                                                 <  2:664056/150‑1 (MQ=255)
      gAAGCGACGCTTGAAGCGCGTCGCAGCCAGGAGATGGACTACTTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCgg                                                                                                                                                 >  1:1715330/1‑150 (MQ=255)
          cGACGCTTGAAGCGCGTCGCAGCCAGGAGATGGACTACTTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGc                                                                                                                                             >  1:1281219/1‑150 (MQ=255)
            aCGCTTGAAGCGCGTCGCAGCCAGGAGATGGACTACTTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGa                                                                                                                                           <  1:1624622/150‑1 (MQ=255)
              gCTTGAAGCGCGTCGCAGCCAGGAGATGGACTACTTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCCGGCGCTTATTGCCTGATGCGAcg                                                                                                                                         <  1:1740893/150‑2 (MQ=255)
              gCTTGAAGCGCGTCGCAGCCAGGAGATGGACTACTTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACa                                                                                                                                         >  1:1878809/1‑150 (MQ=255)
               cTTGAAGCGCGTCGCAGCCAGGAGATGGACTACTTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCCGGCGCTTATTGCCTGATGCGAcgc                                                                                                                                        <  1:1264037/150‑3 (MQ=255)
                 tGAAGCGCGTCGCAGCCAGGAGATGGACTACTTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCa                                                                                                                                      <  1:1532350/150‑1 (MQ=255)
                     gcgcGTCGCAGCCAGGAGATGGACTACTTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGct                                                                                                                                  <  2:1281219/150‑1 (MQ=255)
                         gTCGCAGCCAGGAGATGGACTACTTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTg                                                                                                                              <  2:1740536/150‑1 (MQ=255)
                           cGCAGCCAGGAGATGGACTACTTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTc                                                                                                                            >  1:1701227/1‑150 (MQ=255)
                            gCAGCCAGGAGATGGACTACTTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCa                                                                                                                           >  1:1187610/1‑150 (MQ=255)
                            gCAGCCAGGAGATGGACTACTTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCa                                                                                                                           >  2:1745470/1‑150 (MQ=255)
                            gCAGCCAGGAGATGGACAACTTCATGGAAATATTTGTTCCCAGCATCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCCGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCa                                                                                                                           >  2:1161035/1‑150 (MQ=255)
                             cAGCCAGGAGATGGACTACTTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCaa                                                                                                                          >  2:401369/1‑150 (MQ=255)
                                    gagaTGGACTACTTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGaa                                                                                                                   >  2:609731/1‑150 (MQ=255)
                                      gaTGGACTACTTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAAc                                                                                                                 >  1:1641737/1‑150 (MQ=255)
                                      gaTGGACTACTTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAAc                                                                                                                 <  1:1660316/150‑1 (MQ=255)
                                        tGGACTACTTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAg                                                                                                               >  1:1254424/1‑150 (MQ=255)
                                        tGGACTACTTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAg                                                                                                               >  1:1738023/1‑150 (MQ=255)
                                        tGGACTACTTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGAAGAGATTAGCCAGGAATCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAg                                                                                                               >  2:1596441/1‑150 (MQ=255)
                                          gACTACTTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCa                                                                                                             >  1:662687/1‑150 (MQ=255)
                                              actTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCgg                                                                                                         >  1:286788/1‑150 (MQ=255)
                                               ctTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCggg                                                                                                        <  1:401369/150‑1 (MQ=255)
                                               ctTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCggg                                                                                                        <  2:1701227/150‑1 (MQ=255)
                                                 tCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGat                                                                                                      >  1:1185823/1‑150 (MQ=255)
                                                  cATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGatc                                                                                                     >  1:72757/1‑150 (MQ=255)
                                                    tGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCCGGatcat                                                                                                   >  2:392056/1‑150 (MQ=255)
                                                     ggAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGatcatc                                                                                                  <  1:609731/150‑1 (MQ=255)
                                                      gAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGAtcatca                                                                                                 >  1:1056079/1‑150 (MQ=255)
                                                            tttGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACcacac                                                                                           <  2:1187610/150‑1 (MQ=255)
                                                            tttGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACcacac                                                                                           <  2:1613000/150‑1 (MQ=255)
                                                             ttGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACcacaca                                                                                          <  2:1049535/150‑1 (MQ=255)
                                                                 tCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCCCCACACAAAcc                                                                                      >  1:908846/1‑150 (MQ=255)
                                                                 tCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAAcc                                                                                      >  2:1419868/1‑150 (MQ=255)
                                                                    cAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGcc                                                                                   >  1:589479/1‑150 (MQ=255)
                                                                       cTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATc                                                                                >  2:1957428/1‑150 (MQ=255)
                                                                       cTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATc                                                                                >  1:187957/1‑150 (MQ=255)
                                                                        ttCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCa                                                                               >  2:1719890/1‑150 (MQ=255)
                                                                         tCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCAc                                                                              >  2:785082/1‑150 (MQ=255)
                                                                             tCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACgg                                                                                                                     <  2:1921993/107‑1 (MQ=255)
                                                                             tCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACgg                                                                                                                     >  1:1921993/1‑107 (MQ=255)
                                                                                ttCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATg                                                                       >  1:1574331/1‑150 (MQ=255)
                                                                                  cGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGt                                                                     >  1:794062/1‑150 (MQ=255)
                                                                                   gCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTg                                                                    >  1:412611/1‑150 (MQ=255)
                                                                                      cGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAAc                                                                 <  1:482379/150‑1 (MQ=255)
                                                                                            gaTTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTAcc                                                           >  1:283205/1‑150 (MQ=255)
                                                                                            gaTTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTAcc                                                           >  1:387007/1‑150 (MQ=255)
                                                                                              ttAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCCCACAAACCGCCATCACGATAATGGTGAACGGTAcccc                                                         >  1:1650069/1‑150 (MQ=255)
                                                                                               tAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTAccccc                                                        <  2:1715330/150‑1 (MQ=255)
                                                                                               tAGCCAGGATTCACCGCTGTAAGCAGCCGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCa                                                       <  1:227546/150‑1 (MQ=255)
                                                                                                aGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGAt                                                                           <  2:1985704/130‑1 (MQ=255)
                                                                                                aGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGAt                                                                           >  1:1985704/1‑130 (MQ=255)
                                                                                                  ccAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAAc                                                     <  1:392056/150‑1 (MQ=255)
                                                                                                   cAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACt                                                    >  1:1222500/1‑150 (MQ=255)
                                                                                                   cAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGAAACAGGCAATAAGCGCCGGATGCAACATCAGGCTCTTCTCAAAACGGAAAAAGCATAGGGAACATCACCACACAAACCACAATCACGATAATGATGAACGGAAACCCCAAca                                                    >  2:850650/1‑149 (MQ=255)
                                                                                                     ggATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTc                                                  <  2:286788/150‑1 (MQ=255)
                                                                                                        ttCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTcaca                                               >  1:1919027/1‑150 (MQ=255)
                                                                                                          cACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCCCACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACaaa                                             >  1:1011290/1‑150 (MQ=255)
                                                                                                          cACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACaaa                                             >  1:970788/1‑150 (MQ=255)
                                                                                                          cACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACaaa                                             >  1:1178634/1‑150 (MQ=255)
                                                                                                          cACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACaaa                                             >  1:619950/1‑150 (MQ=255)
                                                                                                          cACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATAGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACaaa                                             >  1:457009/1‑150 (MQ=255)
                                                                                                           aCCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAg                                            >  1:1999116/1‑150 (MQ=255)
                                                                                                                 gTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTg                                      <  1:1506049/150‑1 (MQ=255)
                                                                                                                     gcagcCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGc                                  >  1:542825/1‑150 (MQ=255)
                                                                                                                       agcCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCTCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTg                                <  2:313653/150‑1 (MQ=255)
                                                                                                                       agcCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTg                                <  1:1853974/150‑1 (MQ=255)
                                                                                                                         ccGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGGCATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTa                              >  2:1921480/1‑150 (MQ=255)
                                                                                                                          cGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTaa                             <  1:956254/150‑1 (MQ=255)
                                                                                                                              gTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTc                         <  2:662687/150‑1 (MQ=255)
                                                                                                                                cTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTccc                       >  1:1486211/1‑150 (MQ=255)
                                                                                                                                 tGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCg                      <  2:72757/150‑1 (MQ=255)
                                                                                                                                 gatgTCGCATCCGGCGCTTATTGCCTGATGCGACATCCGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCg                      >  2:1652645/4‑150 (MQ=255)
                                                                                                                                 gatgTCGCATCCGGCGCTTATTGCCTGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCg                      >  2:759849/4‑150 (MQ=255)
                                                                                                                                  gCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCgg                     >  2:1780825/1‑150 (MQ=255)
                                                                                                                                   tgTCGCATCCGGCGCTTATTGCCTGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGa                    >  1:190924/2‑150 (MQ=255)
                                                                                                                                     tCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGAcc                  <  1:1026400/150‑1 (MQ=255)
                                                                                                                                     tCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGAcc                  <  2:1185823/150‑1 (MQ=255)
                                                                                                                                      cGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCt                 <  2:1641737/150‑1 (MQ=255)
                                                                                                                                       gCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTa                >  1:278790/1‑150 (MQ=255)
                                                                                                                                       gCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTa                >  1:90528/1‑150 (MQ=255)
                                                                                                                                         aTCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACAGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTaaa              >  2:1279915/1‑150 (MQ=255)
                                                                                                                                           cAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAAc            <  2:1878809/150‑1 (MQ=255)
                                                                                                                                           cAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAAc            >  1:1364432/1‑150 (MQ=255)
                                                                                                                                           cAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAAc            <  1:1719890/150‑1 (MQ=255)
                                                                                                                                            aGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAAcc           <  2:970788/150‑1 (MQ=255)
                                                                                                                                             ggCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCa          >  1:486805/1‑150 (MQ=255)
                                                                                                                                              gCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAg         <  2:1919027/150‑1 (MQ=255)
                                                                                                                                               cAATAAGAGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAgt        >  2:1273023/1‑150 (MQ=255)
                                                                                                                                                aaTAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAgtg       <  2:1254424/150‑1 (MQ=255)
                                                                                                                                                aaTAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAgtg       <  1:1957428/150‑1 (MQ=255)
                                                                                                                                                aaTAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAgtg       <  2:90528/150‑1 (MQ=255)
                                                                                                                                                    aGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGCATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGTACCTAAAACAAGTGTGtt   <  2:1602611/150‑1 (MQ=255)
                                                                                                                                                     gcgcCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTa  <  2:387007/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGACTGGAAGCGACGCTTGAAGCGCGTCGCAGCCAGGAGATGGACTACTTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTA  >  NC_000913/2409345‑2409643

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