Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,614,808 T→C 7.0% intergenic (+26/+12) focB → / ← yfgO putative formate transporter/putative transporter YfgO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,614,8080TC 7.0% 209.0 / 5.5 72intergenic (+26/+12)focB/yfgOputative formate transporter/putative transporter YfgO
Reads supporting (aligned to +/- strand):  ref base T (35/31);  new base C (3/2);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

CTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCG  >  NC_000913/2614663‑2614956
                                                                                                                                                 |                                                                                                                                                    
cTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCggg                                                                                                                                                  >  2:1470129/1‑150 (MQ=255)
   aCGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGACGAACTGTTTCGCCGGGCTg                                                                                                                                               <  2:1162636/150‑1 (MQ=255)
    cGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGa                                                                                                                                              <  2:335328/150‑1 (MQ=255)
     ggTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGaa                                                                                                                                             <  2:1280102/150‑1 (MQ=255)
        cAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGACGAACTGTTTCGCCGGGCTGatttt                                                                                                                                          <  2:12902/150‑5 (MQ=255)
          gCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGatta                                                                                                                                        >  2:122885/1‑150 (MQ=255)
            cATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGattatt                                                                                                                                      >  1:669367/1‑150 (MQ=255)
             aTTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTc                                                                                                                                     <  2:887078/150‑1 (MQ=255)
                ttattaCCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTg                                                                                                                                  <  1:571241/150‑1 (MQ=255)
                            aTCTGCTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGccc                                                                                                                      <  2:1434611/150‑1 (MQ=255)
                              cTGCTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCAt                                                                                                                    >  2:650882/1‑150 (MQ=255)
                                 cTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCgg                                                                                                                 >  2:402361/1‑150 (MQ=255)
                                 cTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCgg                                                                                                                 >  1:357022/1‑150 (MQ=255)
                                 cTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCgg                                                                                                                 <  1:121242/150‑1 (MQ=255)
                                 cTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCgg                                                                                                                 >  1:362599/1‑150 (MQ=255)
                                   cccGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGc                                                                                                               >  2:448207/1‑150 (MQ=255)
                                    ccGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGcc                                                                                                              <  2:178508/150‑1 (MQ=255)
                                           tGCTGGGTAATATTATCGGCGGTGCGGTGCTGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATg                                                                                                       <  2:728196/150‑1 (MQ=255)
                                           tGCTGGGTAATATTATCGGCGGTGCGGTGCTGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATg                                                                                                       >  1:879498/1‑150 (MQ=255)
                                               gggTAATATTATCGGCGGTGCGGTGCTGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATc                                                                                                   >  1:842258/1‑150 (MQ=255)
                                                  tAATATTATCGGCGGTGCGGTGCTGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACg                                                                                                <  1:448207/150‑1 (MQ=255)
                                                   aaTATTATCGGCGGTGCGGTGCGGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCACGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACgg                                                                                               <  2:774644/150‑1 (MQ=255)
                                                      tttatCGGCGGTGCGGTTCTTGTGATTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCTTGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCtt                                                                                            <  2:552535/149‑1 (MQ=255)
                                                        tatCGGCGGTGCGGTGCTGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTg                                                                                          <  1:1010597/150‑1 (MQ=255)
                                                         atCGGCGGTGCGGTGCTGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGa                                                                                         <  1:1407345/150‑1 (MQ=255)
                                                         atCGGCGGTGCGGTGCTGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGa                                                                                         <  1:122885/150‑1 (MQ=255)
                                                           cggcggTGCGGTGCTGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATc                                                                                       >  2:1469644/1‑150 (MQ=255)
                                                           cggcggTGCGGTGCTGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATc                                                                                       <  1:101899/150‑1 (MQ=255)
                                                           cggcggTGCGGTGCTGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATc                                                                                       >  2:989413/1‑150 (MQ=255)
                                                            ggcggTGCGGTGCTGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCa                                                                                      <  2:1316943/150‑1 (MQ=255)
                                                                 ggcggtgcTGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTc                                                                                 <  2:987161/149‑1 (MQ=255)
                                                                       gcTGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGACTCACGGCTTTGATCAGCGTCGCCAAc                                                                           <  1:650882/150‑1 (MQ=255)
                                                                         tGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACgg                                                                         <  2:411670/150‑1 (MQ=255)
                                                                         tGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACgg                                                                         >  1:272513/1‑150 (MQ=255)
                                                                           gTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGaa                                                                       >  1:1108268/1‑150 (MQ=255)
                                                                             gAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATg                                                                     >  2:778274/1‑150 (MQ=255)
                                                                               gTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGc                                                                   <  1:1450554/150‑1 (MQ=255)
                                                                                   gtgtTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCgaag                                                               >  1:1350655/1‑150 (MQ=255)
                                                                                   gtgtTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCgaag                                                               >  2:1461496/1‑150 (MQ=255)
                                                                                    tgtTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCgaaga                                                              >  2:523193/1‑150 (MQ=255)
                                                                                    tgtTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCgaaga                                                              >  1:1442387/1‑150 (MQ=255)
                                                                                    tgtTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCgaaga                                                              >  2:1307770/1‑150 (MQ=255)
                                                                                     gtTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCgaagaa                                                             >  1:34618/1‑150 (MQ=255)
                                                                                      ttATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAAt                                                            <  1:1461496/150‑1 (MQ=255)
                                                                                       tATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATa                                                           <  2:1442387/150‑1 (MQ=255)
                                                                                        aTCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATAc                                                          >  1:1039230/1‑150 (MQ=255)
                                                                                           gggCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGcc                                                       <  1:1469644/150‑1 (MQ=255)
                                                                                           gggCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGcc                                                       >  1:1140124/1‑150 (MQ=255)
                                                                                           gggCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGcc                                                       >  2:522577/1‑150 (MQ=255)
                                                                                            ggCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGccc                                                      <  2:1039230/150‑1 (MQ=255)
                                                                                                atttatttACGTCAGGAACCCTGAAAAATCAGCCCGGCGATACAGTTCGGCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAAGGGAATGGCGAAGGATTGGCCCGAgc                                                  >  2:803332/1‑149 (MQ=255)
                                                                                                 tttatttaCGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGaa                                                 <  2:414529/150‑1 (MQ=255)
                                                                                                    atttaCGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATcc                                              <  1:1451526/150‑1 (MQ=255)
                                                                                                    atttaCGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATcc                                              >  2:787225/1‑150 (MQ=255)
                                                                                                      ttaCGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATcccc                                            <  1:1231114/150‑1 (MQ=255)
                                                                                                       taCGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCca                                           <  2:879498/150‑1 (MQ=255)
                                                                                                         cGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCcaca                                         >  1:1241690/1‑150 (MQ=255)
                                                                                                              ggAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAacca                                    >  2:43622/1‑150 (MQ=255)
                                                                                                                    cTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGGATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAg                              >  2:709731/1‑150 (MQ=255)
                                                                                                                    cTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAg                              >  1:669432/1‑150 (MQ=255)
                                                                                                                    cTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAg                              >  2:1454641/1‑150 (MQ=255)
                                                                                                                    cTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAg                              <  1:1470129/150‑1 (MQ=255)
                                                                                                                    cTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAg                              >  2:143965/1‑150 (MQ=255)
                                                                                                                        atcttCAGCCCGACGAACTGTTTCGCCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCCCGGCTTTGATCAGCGTCGCCCACGGAATGGCGATGAATACGCCCCAGAATCCCCACCAACCACCGAAGATcc                          >  2:494702/5‑149 (MQ=255)
                                                                                                                         tcttCAGCCCGACGAACTGTTTCGCCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATcac                         >  1:133190/4‑150 (MQ=255)
                                                                                                                           ttCAGCCCGACGAACTGTTTCGCCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATcaccc                       >  1:1501853/2‑149 (MQ=255)
                                                                                                                           aTCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATcacca                       <  2:1108268/150‑1 (MQ=255)
                                                                                                                             cAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCaccacc                     >  1:288475/1‑150 (MQ=255)
                                                                                                                                 ccGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATa                 >  1:731271/1‑150 (MQ=255)
                                                                                                                                     cGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAaat             >  2:98929/1‑150 (MQ=255)
                                                                                                                                      gAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAaata            <  1:710294/150‑1 (MQ=255)
                                                                                                                                       aaaCAGTTACTCCGGCTTAAGATTATTCTACCTCGATATGCACATCGGGCCATGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAaataa           >  1:1423620/1‑150 (MQ=255)
                                                                                                                                                gtcgGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCg  >  1:62367/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTGGTGAGTATGTGTTATCGGGCTATTTATTTACGTCAGGAACCCTGAAAAATCAGCCCGGCGAAACAGTTCGTCGGGCTGAAGATTATTCTTGCGCGATTTGCCCATCGGGCCAGGCATGAATCACGGCTTTGATCAGCGTCGCCAACGGAATGGCGAAGAATACGCCCCAGAATCCCCACAAACCACCGAAGATCACCACCGATAAAATAATCACCAGCG  >  NC_000913/2614663‑2614956

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: