Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 774,633 G→T 6.5% intergenic (+31/‑119) ybgE → / → ybgC cyd operon protein YbgE/esterase/thioesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913774,6330GT 6.5% 326.4 / 6.1 99intergenic (+31/‑119)ybgE/ybgCcyd operon protein YbgE/esterase/thioesterase
Reads supporting (aligned to +/- strand):  ref base G (51/40);  new base T (3/5);  total (54/45)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTC  >  NC_000913/774489‑774782
                                                                                                                                                |                                                                                                                                                     
ggTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCa                                                                                                                                                  >  1:143092/1‑150 (MQ=255)
 gtgtGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGATCCAg                                                                                                                                                 >  2:1405009/1‑150 (MQ=255)
  tgtgATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAgg                                                                                                                                                <  2:1018344/150‑1 (MQ=255)
            ggCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAAtt                                                                                                                                      >  1:1774075/1‑150 (MQ=255)
            ggCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAAtt                                                                                                                                      <  1:114248/150‑1 (MQ=255)
               gTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAAttttt                                                                                                                                   <  2:320299/150‑1 (MQ=255)
               gTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAAttttt                                                                                                                                   <  2:1684203/150‑1 (MQ=255)
                tGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAATTTATGGGCCTGGAGCGGCCCATAATTTTTTa                                                                                                                                  <  2:73975/150‑1 (MQ=38)
                 gggCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAATTTATGGGCCTGGAGCGGCCCATAATTTTTTaa                                                                                                                                 <  1:831150/150‑2 (MQ=38)
                  ggCTTTCGTCCGCAGAAGGTTCTTTGGCACGGGATTTTTTGCCCATTGCTTGCCGATATTCTTCTCATTGTCGGGCTGATTTTTTTCTTCTTTTTATTCTGATTTCTCTTTAAAATATTATTGTCCGCTCTATTCCCATAAATTTTTatg                                                                                                                                <  2:1689798/150‑3 (MQ=255)
                   gCTTTCGTCCGCAGAAGGTGCTTTGGCAACGCATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAATTTATGGGCCTGGAGCGGCCCATAATTTTTTaaag                                                                                                                               <  2:1403933/150‑4 (MQ=37)
                    cTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTGTTTTAAATCAGAATTCTCTTTAAAAATTTATGGGCCTGGAGCGGCCCATAATTTTTTaaaga                                                                                                                              <  1:407375/150‑5 (MQ=37)
                    cTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAATTTATGGGCCTGGAGCGGCCCATAATTTTTTaaaga                                                                                                                              >  1:1450324/1‑146 (MQ=38)
                    cTTTCGTCCGCAGAAGGTGCTTTGGCAAGCGATTTTTTGCCCATTGCTTCCCGATATTGTTCTCATTGTCCGGCGGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAAATATGGGCCGCTCCAGGCCCATAAATTTTTACGCt                                                                                                                              <  1:1842997/150‑1 (MQ=255)
                      ttCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTcc                                                                                                                            <  1:1122360/150‑1 (MQ=255)
                         gTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCtt                                                                                                                         >  2:312858/1‑150 (MQ=255)
                          tCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTa                                                                                                                        >  1:781538/1‑150 (MQ=255)
                            cGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAAc                                                                                                                      >  1:1567194/1‑150 (MQ=255)
                            cGCAGAAGGCGCTTTGGCAAGGCATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAAc                                                                                                                      <  2:1002935/150‑1 (MQ=255)
                             gCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACt                                                                                                                     >  2:1759155/1‑150 (MQ=255)
                             gCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACt                                                                                                                     >  2:1771255/1‑150 (MQ=255)
                                    gTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTc                                                                                                              >  2:10767/1‑150 (MQ=255)
                                    gTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTc                                                                                                              >  1:1590577/1‑150 (MQ=255)
                                       gtttGGCAATTGATTTTTTGCCCATTGCTTGCCGATTTTGTTCTCATTTTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCAtt                                                                                                           <  2:1594335/149‑1 (MQ=255)
                                       cTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCAtt                                                                                                           >  1:1095931/1‑150 (MQ=255)
                                         ttGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTcc                                                                                                         >  2:1608167/1‑150 (MQ=255)
                                             cAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCaaa                                                                                                     >  1:949348/1‑150 (MQ=255)
                                              aaGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAAc                                                                                                    <  1:363381/150‑1 (MQ=255)
                                               aGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAAcc                                                                                                   <  1:431824/150‑1 (MQ=255)
                                               aGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAAcc                                                                                                   <  1:720205/150‑1 (MQ=255)
                                               aGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAAcc                                                                                                   >  2:137036/1‑150 (MQ=255)
                                                gggATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTc                                                                                                          <  2:425538/142‑1 (MQ=255)
                                                gggATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTc                                                                                                          >  1:425538/1‑142 (MQ=255)
                                                gggATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCt                                                                                                  >  2:128717/1‑150 (MQ=255)
                                                gggATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTCCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCt                                                                                                  >  1:1418553/1‑150 (MQ=255)
                                                 ggATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTc                                                                                                 <  1:1784317/150‑1 (MQ=255)
                                                  gATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCa                                                                                                >  1:1578925/1‑150 (MQ=255)
                                                      ttttGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATcg                                                                                            >  1:224949/1‑150 (MQ=255)
                                                      ttttGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTATTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATcg                                                                                            <  1:444497/150‑1 (MQ=255)
                                                         tGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATcgcgc                                                                                         >  1:568647/1‑150 (MQ=255)
                                                         tGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATcgcgc                                                                                         <  1:1405009/150‑1 (MQ=255)
                                                           cccATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCGTAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATcgcgcgc                                                                                       >  2:749586/1‑150 (MQ=255)
                                                            ccATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATcgcgcgcg                                                                                      <  1:312858/150‑1 (MQ=255)
                                                              aTTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGta                                                                                    <  1:159153/150‑1 (MQ=255)
                                                                  cttgcCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTAtagt                                                                                <  1:24664/150‑1 (MQ=255)
                                                                  cttgcCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTAtagt                                                                                >  2:1227473/1‑150 (MQ=255)
                                                                    tgcCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAAt                                                                                              >  1:1473689/1‑134 (MQ=255)
                                                                    tgcCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAAt                                                                                              <  2:1473689/134‑1 (MQ=255)
                                                                     gcCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCa                                                                                                >  1:1100690/1‑131 (MQ=255)
                                                                     gcCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCa                                                                                                <  2:1100690/131‑1 (MQ=255)
                                                                             tGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTa                                                                     <  2:224949/150‑1 (MQ=255)
                                                                                    aTTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGctc                                                              <  1:1859162/150‑1 (MQ=255)
                                                                                     ttGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGctct                                                             >  1:1530391/1‑150 (MQ=255)
                                                                                       gTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTa                                                                     >  1:919374/1‑140 (MQ=255)
                                                                                       gTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTa                                                                     <  2:919374/140‑1 (MQ=255)
                                                                                           ggCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACttt                                                       >  1:946076/1‑150 (MQ=255)
                                                                                           ggCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGACCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACttt                                                       >  1:1375625/1‑150 (MQ=255)
                                                                                            gCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCTTTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACtttt                                                      <  2:1134047/150‑1 (MQ=255)
                                                                                             cTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTaaaa                                                                  >  2:353394/1‑137 (MQ=255)
                                                                                             cTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTaaaa                                                                  <  1:353394/137‑1 (MQ=255)
                                                                                             cTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTACCAGCGTTTAAAAGCTCTAACTTttg                                                     >  2:1670296/1‑150 (MQ=255)
                                                                                              tGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTttgt                                                    <  1:1518372/150‑1 (MQ=255)
                                                                                               gATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAAttttt                                                                                                                                   >  2:665051/1‑70 (MQ=255)
                                                                                               gATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAAttttt                                                                                                                                   <  1:665051/70‑1 (MQ=255)
                                                                                               gATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTttgtt                                                   <  1:137036/150‑1 (MQ=255)
                                                                                                aTTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTttgttg                                                  <  1:1771255/150‑1 (MQ=255)
                                                                                                 ttttCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGc                                                 <  2:1774075/150‑1 (MQ=255)
                                                                                                    tcttcttctTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTTTAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCAtt                                              <  2:1567194/150‑1 (MQ=255)
                                                                                                         ttcttTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCgg                                         >  1:25352/1‑150 (MQ=255)
                                                                                                         ttcttTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCgg                                         <  2:1578925/150‑1 (MQ=255)
                                                                                                          tcttTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCggg                                        <  2:1648745/150‑1 (MQ=255)
                                                                                                           cttTTAAATCAGAATTCTCTTTAAAAAATTTTGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGa                                       >  2:1016805/1‑150 (MQ=255)
                                                                                                           cttTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGa                                       >  2:1531958/1‑150 (MQ=255)
                                                                                                           cttTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGa                                       >  1:592749/1‑150 (MQ=255)
                                                                                                                aaaTCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTaa                                  >  1:1579554/1‑150 (MQ=255)
                                                                                                                aaaTCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTCCGGGATGTaaa                                 >  1:1946410/1‑150 (MQ=255)
                                                                                                                 aaTCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTaaa                                 >  1:462750/1‑150 (MQ=255)
                                                                                                                  aTCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAg                                >  2:212819/1‑150 (MQ=255)
                                                                                                                  aTCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAg                                >  2:1980912/1‑150 (MQ=255)
                                                                                                                  aTCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAg                                >  1:1169941/1‑150 (MQ=255)
                                                                                                                   acAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGt                               <  2:556747/149‑1 (MQ=255)
                                                                                                                    cAGAATTCTCTTTAAAAAATTATGGGCCTCTCCAGGCCCATAAATTTTTACTCTCCCTTAACTTTTTATAATTACAAAACCTCAATAGCGAGAGTATATTATAAGCTTTTAAAATCTCTAAATTTTGTTGCATTAACTTGATGTAAATTg                              <  2:1925557/150‑1 (MQ=255)
                                                                                                                     aGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGa                             >  2:1592449/1‑150 (MQ=255)
                                                                                                                     aGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGa                             >  2:43776/1‑150 (MQ=255)
                                                                                                                        aTTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATa                          >  1:1620700/1‑150 (MQ=255)
                                                                                                                             ctTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACg                     >  1:170945/1‑150 (MQ=255)
                                                                                                                              tttAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGc                    <  2:1095931/150‑1 (MQ=255)
                                                                                                                               ttAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCGCTTAACTTGCCCTCATTCCCAAACCTCAAGCGCGCGCGTAGAGTAGCAGCGGTTAAAAGCTCTAACTTTTGTTGCATTCCCGGGATGTATAGTGAATACAACGCt                   >  1:1908120/1‑150 (MQ=255)
                                                                                                                               gtAAAAATTTATGGGCCTGGAGCGGCCCATAAATTTTTACGCTCCCTTAACTTGACCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGAAGCGTTTAAAAGATATAAATTTTGTTAAATTACCGGGAAGTAAAATGAATaaaaataa                   >  1:851817/2‑146 (MQ=25)
                                                                                                                                tAAAAATTTATGGGCCTGGAGCGGCCCATAATTTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTg                  >  1:1487194/1‑150 (MQ=38)
                                                                                                                                  aaaaaTTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGtt                >  2:1371656/1‑150 (MQ=255)
                                                                                                                                  aaaaaTTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGtt                >  1:115393/1‑150 (MQ=255)
                                                                                                                                  aaaaaTTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGtt                >  1:1345327/1‑150 (MQ=255)
                                                                                                                                   aaaaTTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGttt               >  2:855823/1‑150 (MQ=255)
                                                                                                                                    aaaTTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAACCCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTc              >  1:1658219/1‑150 (MQ=255)
                                                                                                                                         aTGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATgg         <  2:568647/150‑1 (MQ=255)
                                                                                                                                         aTGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATgg         <  1:128717/150‑1 (MQ=255)
                                                                                                                                           gggCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGcc       <  1:212819/150‑1 (MQ=255)
                                                                                                                                             gCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCgg     <  1:10767/150‑1 (MQ=255)
                                                                                                                                                gCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTc  >  2:131383/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTC  >  NC_000913/774489‑774782

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