Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,409,489 A→G 5.2% intergenic (+28/+31) 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,4890AG 5.2% 298.1 / 4.8 96intergenic (+28/+31)yfbR/yfbSdCMP phosphohydrolase/putative transporter YfbS
Reads supporting (aligned to +/- strand):  ref base A (56/35);  new base G (3/2);  total (59/37)
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.72e-01

CGCGACTGGAAGCGACGCTTGAAGCGCGTCGCAGCCAGGAGATGGACTACTTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTG  >  NC_000913/2409343‑2409638
                                                                                                                                                  |                                                                                                                                                     
cgcgACTGGAAGCGACGCTTGAAGCGCGTCGCAGCCAGGAGATGGACTACTTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATa                                                                                                                                                    <  1:60382/150‑1 (MQ=255)
  cgACTGGAAGCGACGCTTGAAGCGCGTCGCAGCCAGGAGATGGACTACTTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATacg                                                                                                                                                  >  1:1000261/1‑148 (MQ=255)
  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                                                                                                                                            >  1:1715330/1‑150 (MQ=255)
        gAAGCGACGCTTGAAGCGCGTCGCAGCCAGGAGATGGACTACTTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCgg                                                                                                                                            <  2:664056/150‑1 (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                                                                                                                      >  2:1745470/1‑150 (MQ=255)
                              gCAGCCAGGAGATGGACTACTTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCa                                                                                                                      >  1:1187610/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:1613000/150‑1 (MQ=255)
                                                              tttGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACcacac                                                                                      <  2:1187610/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:387007/1‑150 (MQ=255)
                                                                                              gaTTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTAcc                                                      >  1:283205/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                                                                      >  1:1985704/1‑130 (MQ=255)
                                                                                                  aGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGAt                                                                      <  2:1985704/130‑1 (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:1178634/1‑150 (MQ=255)
                                                                                                            cACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACaaa                                        >  1:970788/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       >  1:1364432/1‑150 (MQ=255)
                                                                                                                                             cAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAAc       <  1:1719890/150‑1 (MQ=255)
                                                                                                                                             cAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAAc       <  2:1878809/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  <  1:1957428/150‑1 (MQ=255)
                                                                                                                                                  aaTAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAgtg  <  2:1254424/150‑1 (MQ=255)
                                                                                                                                                  aaTAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAgtg  <  2:90528/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CGCGACTGGAAGCGACGCTTGAAGCGCGTCGCAGCCAGGAGATGGACTACTTCATGGAAATATTTGTTCCCAGCTTCCATCTTTCGCTCGATGAGATTAGCCAGGATTCACCGCTGTAAGCAGCCGGAGTCTGCGTCGCATCAGGCAATAAGCGCCGGATGCGACATCAGGCTCTTGTCAAAACGGAAACAGCATCGGGATCATCACCACACAAACCGCCATCACGATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTG  >  NC_000913/2409343‑2409638

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: