Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,483,781 G→C 5.6% G194R (GGC→CGC) ‡ yheT → putative hydrolase YheT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,483,7810GC 5.6% 304.1 / 6.8 89G194R (GGC→CGC) ‡yheTputative hydrolase YheT
Reads supporting (aligned to +/- strand):  ref base G (48/36);  new base C (3/2);  total (51/38)
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.67e-01

ATGCGCCAACGGCTGCCGTCGGCTATTCGCTCGGCGGTAATATGCTGGCCTGTTTGCTGGCAAAGGAAGGCAATGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTACGTCGCATCCGTGAATTTGACGATCTCATCACC  >  NC_000913/3483632‑3483930
                                                                                                                                                     |                                                                                                                                                     
aTGCGCCAACGGCTGCCGTCGGCTATTCGCTCGGCGGTAATATGCTGGCCTGTTTGCTGGCAAAGGAAGGCAATGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAgg                                                                                                                                                       >  1:957328/1‑150 (MQ=255)
aTGCGCCAACGGCTGCCGTCGGCTATTCGCTAGGCGGTAATATGCTGGCATGTTTGATGGCAAAGGAAGGCAATGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTACGCCCTTTATGCTGGAAGCCTGTAGATATCATATGTAAAAgg                                                                                                                                                       >  2:655666/1‑150 (MQ=255)
 tGCGCCAACGGCTGCCGTCGGCTATTCGCTCGGCGGTAATATGCTGGCCTGTTTGCTGGCAAAGGAAGGCAATGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAggg                                                                                                                                                      >  2:733813/1‑150 (MQ=255)
  gcgcCAACGGCTGCCGTCGGCTATTCGCTCGGCGGTAATATGCTGGCCTGTTTGCTGGCAAAGGAAGGCAATGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGc                                                                                                                                                     <  1:686372/150‑1 (MQ=255)
     ccAACGGCTGCCGTCGGCTATTCGCTCGGCGGTAATATGCTGGCCTGTTTGCTGGCAAAGGAAGGCAATGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCttt                                                                                                                                                  >  2:1152397/1‑150 (MQ=255)
            cTGCCGTCGGCTATTCGCTCGGCGGTAATATGCTGGCCTGTTTGCTGGCAAAGGAAGGCAATGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCcgcg                                                                                                                                           <  1:463351/150‑1 (MQ=255)
             tGCCGTCGGCTATTCGCTCGGCGGTAATATGCTGGCCTGTTTGCTGGCAAAGGAAGGCAATGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGCCCTTTTccatatg                                                                                                                                          <  2:1472945/150‑6 (MQ=255)
              gCCGTCGGCTATTCGCTCGGCGGTAATATGCTGGCCTGTTTGCTGGCAAAGGAAGGCAATGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGtt                                                                                                                                         <  1:981497/150‑1 (MQ=255)
               ccGTCGGCTATTCGCTCGGCGGTAATATGCTGGCCTGTTTGCTGGCAAAGGAAGGCAATGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGCCCTTTTccatatgat                                                                                                                                        <  2:1218033/150‑8 (MQ=255)
                   cGGCTATTCGCTCGGCGGTAATATGCTGGCCTGTTTGCTGGCAAAGGAAGGCAATGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCa                                                                                                                                    <  2:331792/150‑1 (MQ=255)
                   cGGCTATTCGCTCGGCGGTAATATGCTGGCCTGTTTGCTGGCAAAGGAAGGCAATGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCa                                                                                                                                    >  1:1348277/1‑150 (MQ=255)
                      gtATTCGCTCGGCGGTAATATGCTGTCCTGTTTGCTTGCAAAGGAATTCAATGATCTCCCTGTTGATTCGGCTGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAATGGCTTTTCCCGCGTTTATCAGCg                                                                                                                                 <  2:917535/149‑1 (MQ=255)
                       tATTCGCTCGGCGGTAATATGCTGGCCTGTTTGCTGGCAAAGGAAGGCAATGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGt                                                                                                                                >  1:821093/1‑150 (MQ=255)
                           cGCTCGGCGGTAATATGCTGGCCTGTTTGCTGGCAAAGGAAGGCAATGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACt                                                                                                                            <  1:1268824/150‑1 (MQ=255)
                           cGCTCGGCGGTAATATGCTGGCCTGTTTGCTGGCAAAGGAAGGCAATGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCAGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACt                                                                                                                            >  1:352264/1‑150 (MQ=255)
                               cggcggTAATATGCTGGCCTGTTTGCTGGCAAAGGAAGGCAATGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCt                                                                                                                        <  2:1490518/150‑1 (MQ=255)
                                   ggTAATATGCTGGCCTGTTTGCTGGCAAAGGAAGGCAATGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAAc                                                                                                                    <  1:1107220/150‑1 (MQ=255)
                                   ggTAATATGCCGGCCTGTTTGCTGGCAAAGGAAGGCAATGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAAc                                                                                                                    <  2:826212/150‑1 (MQ=255)
                                    gTAATATGCTGGCCTGTTTGCTGGCAAAGGAAGGCAATGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAAcc                                                                                                                   <  1:1017751/150‑1 (MQ=255)
                                      aaTATGCTGGCCTGTTTGCTGGCAAAGGAAGGCAATGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTg                                                                                                                 <  2:1180069/150‑1 (MQ=255)
                                      aaTATGCTGGCCTGTTTGCTGGCAAAGGAAGGCAATGATCTCCCGGTTGAAGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTg                                                                                                                 >  1:231266/1‑150 (MQ=255)
                                        tatGCTGGCCTGTTTGCTGGCAAAGGAAGGCAATGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGtt                                                                                                               <  2:1348277/150‑1 (MQ=255)
                                         atGCTGGCCTGTTTGCTGGCAAAGGAAGGCAATGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTa                                                                                                              <  2:1313827/150‑1 (MQ=255)
                                              ggCCTGTTTGCTGGCAAAGGAAGGCAATGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGc                                                                                                         >  1:717872/1‑150 (MQ=255)
                                                ccTGTTTGCTGGCAAAGGAAGGCAATGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCa                                                                                                       >  1:807303/1‑150 (MQ=255)
                                                   gTTTGCTGGCAAAGGAAGGCAATGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATg                                                                                                    >  1:1277550/1‑150 (MQ=255)
                                                    tttGCTGGCAAAGGAAGGCAATGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGc                                                                                                   <  1:1016092/150‑1 (MQ=255)
                                                          ggCAAAGGAAGGCAATGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCcgcgcg                                                                                             >  1:415323/1‑150 (MQ=255)
                                                          ggCAAAGGAAGGCAATGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCcgcgcg                                                                                             >  1:420561/1‑150 (MQ=255)
                                                           gCAAAGGAAGGCAATGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGGTTATGCTGGAAGCCGGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCgcgcgc                                                                                            <  1:1365070/150‑1 (MQ=255)
                                                                  aaggCAATGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTgg                                                                                     >  1:1369018/1‑150 (MQ=255)
                                                                    ggCAATGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCa                                                                                   >  1:89381/1‑150 (MQ=255)
                                                                     gCAATGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTTGCTATGATATGGAAAATGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAg                                                                                  <  1:655666/150‑1 (MQ=255)
                                                                      cAATGATCTCCCGGTTCATGCCGCGGTGATTGTCTCTGCGCCGATTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTGGCTGATCCTCTTAAAAGCCAATGCCGCGCGCAAGCTGGCacc                                                                                 >  1:304360/1‑148 (MQ=255)
                                                                         tGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTa                                                                              <  2:98959/150‑1 (MQ=255)
                                                                            tctcCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTAccc                                                                           <  1:70058/150‑1 (MQ=255)
                                                                            tctcCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTAccc                                                                           >  2:1400208/1‑150 (MQ=255)
                                                                             ctcCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTAcccc                                                                          >  2:1399771/1‑150 (MQ=255)
                                                                                  ggTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAAc                                                                     >  1:784681/1‑150 (MQ=255)
                                                                                     tGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGc                                                                                    <  1:1287132/132‑1 (MQ=255)
                                                                                     tGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGc                                                                                    >  2:1287132/1‑132 (MQ=255)
                                                                                      gATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTg                                                                 >  2:1494534/1‑150 (MQ=255)
                                                                                       aTGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGc                                                                >  1:1007221/1‑150 (MQ=255)
                                                                                         gcggcgGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCg                                                              >  2:295037/1‑150 (MQ=255)
                                                                                         gcggcgGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCg                                                              <  2:1163714/150‑1 (MQ=255)
                                                                                            gcggTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGAtt                                                           >  1:953640/1‑150 (MQ=255)
                                                                                                 gATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAAtct                                                      >  1:1127049/1‑150 (MQ=255)
                                                                                                  aTTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAAtctc                                                     <  2:957328/150‑1 (MQ=255)
                                                                                                   ttGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTcg                                                    <  2:405283/150‑1 (MQ=255)
                                                                                                   ttGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTcg                                                    >  2:276511/1‑150 (MQ=255)
                                                                                                       ctctGCGCTGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCa                                                >  2:1147965/1‑150 (MQ=255)
                                                                                                       ctctGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCa                                                >  2:1120905/1‑150 (MQ=255)
                                                                                                       ctctGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCa                                                >  2:106993/1‑150 (MQ=255)
                                                                                                           gcgcCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTa                                            >  2:945491/1‑150 (MQ=255)
                                                                                                              ccGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTaaaa                                         <  2:1434551/150‑1 (MQ=255)
                                                                                                               cGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAAt                                        >  2:1142141/1‑150 (MQ=255)
                                                                                                                gTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATc                                       <  1:1152397/150‑1 (MQ=255)
                                                                                                                gTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATc                                       <  2:65308/150‑1 (MQ=255)
                                                                                                                  ttATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCgg                                     >  2:1486361/1‑150 (MQ=255)
                                                                                                                   tATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGt                                    >  2:599184/1‑150 (MQ=255)
                                                                                                                   tATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCAGAGAGCAAGCTCGCAGCCTAAACCGGAACGATGCCGATTAATCTCGCGCAGTTAAAATCgtt                                    >  1:1018605/1‑148 (MQ=255)
                                                                                                                     tGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTAc                                  >  1:371148/1‑150 (MQ=255)
                                                                                                                     tGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTAc                                  >  2:938468/1‑150 (MQ=255)
                                                                                                                      gCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTACg                                 >  1:276050/1‑150 (MQ=255)
                                                                                                                       cTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTACGt                                >  1:1512792/1‑150 (MQ=255)
                                                                                                                         ggAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTCAAAATCGGTACGTCg                              <  1:1235525/150‑1 (MQ=255)
                                                                                                                             gCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTACGTCGCATc                          >  1:182417/1‑150 (MQ=255)
                                                                                                                              ccTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTACGTCGCATcc                         >  2:318833/1‑150 (MQ=255)
                                                                                                                               cTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTACGTCGCATCCg                        >  1:246875/1‑150 (MQ=255)
                                                                                                                                tGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTACGTCGCATCCGt                       <  2:1277550/150‑1 (MQ=255)
                                                                                                                                 gTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTACGTCGCATCCGTg                      <  2:89381/150‑1 (MQ=255)
                                                                                                                                    gCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTACGTCGCATCCGTGAAt                   >  2:1309050/1‑150 (MQ=255)
                                                                                                                                     cTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTACGTCGCATCCGTGAAtt                  >  1:860961/1‑150 (MQ=255)
                                                                                                                                      tATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTACGTCGCATCCGTGAAttt                 <  2:415323/150‑1 (MQ=255)
                                                                                                                                      tATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTACGTCGCATCCGTGAAttt                 <  2:420561/150‑1 (MQ=255)
                                                                                                                                      aaTCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTACGTCGCATCCGTGAAttt                 >  1:1343625/2‑150 (MQ=255)
                                                                                                                                       aTCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTACGTCGCATCCGTGAATTTg                <  2:717872/150‑1 (MQ=255)
                                                                                                                                         acgcgggAAAAGCCCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTACGTCGCATCCGTGAATTTGAc              >  1:823478/6‑150 (MQ=255)
                                                                                                                                          atatGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTACGTCGCATCCGTGAATTTGACg             >  1:528656/1‑150 (MQ=255)
                                                                                                                                            cgggAAAAGCCCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTACGTCGCATCCGTGAATTTGACGAt           >  2:291896/3‑150 (MQ=255)
                                                                                                                                            cgggAAAAGCCCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTACGTCGCATCCGTGAATTTGACGAt           >  1:692700/3‑150 (MQ=255)
                                                                                                                                            atGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTACGTCGCATCCGTGAATTTGACGAt           >  2:372213/1‑150 (MQ=255)
                                                                                                                                             tGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTACGTCGCATCCGTGAATTTGACGAtc          <  2:231266/150‑1 (MQ=255)
                                                                                                                                             tGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTACGTCGCATCCGTGAATTTGACGAtc          >  2:1502642/1‑150 (MQ=255)
                                                                                                                                               gAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTACGTCGCATCCGTGAATTTGACGAtctc        >  2:347633/1‑150 (MQ=255)
                                                                                                                                               gAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTACGTCGCATCCGTGAATTTGACGAtctc        <  2:1007221/150‑1 (MQ=255)
                                                                                                                                               gAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTACGTCGCATCCGTGAATTTGACGAtctc        >  2:1369048/1‑150 (MQ=255)
                                                                                                                                                aaaaGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCTCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTACGTCGCATCCGTGAATTTGACGATCtca       <  1:1120905/150‑1 (MQ=255)
                                                                                                                                                aaaaGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTACGTCGCATCCGTGAATTTGAGGATCtca       <  2:1343625/150‑1 (MQ=255)
                                                                                                                                                 aaaGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTACGTCGCATCCGTGAATTTGACGATCtcat      >  2:377356/1‑150 (MQ=255)
                                                                                                                                                 aaaGGGCATTTCCCGAGTTTATCAGCGTTACTTGCTCAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCAGATTAATCTCGCGCAGTTAAAATCGGAACGTCGCATCCGTGAATTTGACGATCtcaa      >  2:547684/1‑149 (MQ=255)
                                                                                                                                                  aaGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTACGTCGCATCCGTGAATTTGACGATCtcatc     <  2:1127049/150‑1 (MQ=255)
                                                                                                                                                    gggCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTACGTCGCATCCGTGAATTTGACGATCTCATCAc   >  1:1201528/1‑150 (MQ=255)
                                                                                                                                                     ggCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTACGTCGCATCCGTGAATTTGACGATCTCATCAcc  >  2:73212/1‑150 (MQ=255)
                                                                                                                                                     ggCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTACGTCGCATCCGTGAATTTGACGATCTCATCAcc  >  2:1055936/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATGCGCCAACGGCTGCCGTCGGCTATTCGCTCGGCGGTAATATGCTGGCCTGTTTGCTGGCAAAGGAAGGCAATGATCTCCCGGTTGATGCGGCGGTGATTGTCTCTGCGCCGTTTATGCTGGAAGCCTGTAGCTATCATATGGAAAAGGGCTTTTCCCGCGTTTATCAGCGTTACTTGCTGAACCTGTTAAAAGCCAATGCCGCGCGCAAGCTGGCAGCCTACCCCGGAACGCTGCCGATTAATCTCGCGCAGTTAAAATCGGTACGTCGCATCCGTGAATTTGACGATCTCATCACC  >  NC_000913/3483632‑3483930

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: