Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,299,267:1 +A 5.2% intergenic (+22/‑1915) ychE → / → oppA UPF0056 family inner membrane protein/oligopeptide ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,299,2671.A 5.2% 459.1 / 19.4 115intergenic (+22/‑1915)ychE/oppAUPF0056 family inner membrane protein/oligopeptide ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base . (52/57);  new base A (3/3);  total (55/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

Only 100 of 116 total aligned reads displayed.

TGGTACGGGTTTTACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTT  >  NC_000913/1299118‑1299413
                                                                                                                                                      |                                                                                                                                                  
tGGTACGGGTTTTACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGc                                                                                                                                                     >  2:1013768/1‑150 (MQ=255)
     cGGGTTTTACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCa                                                                                                                                               >  2:414916/1‑150 (MQ=255)
        gTTTTACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGc                                                                                                                                            >  1:1239638/1‑150 (MQ=255)
         ttttACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCt                                                                                                                                           <  1:628653/150‑1 (MQ=255)
          tttACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTc                                                                                                                                          >  2:301567/1‑150 (MQ=255)
          tttACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTc                                                                                                                                          >  2:1833900/1‑150 (MQ=255)
           ttACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGTTGGCATTGGGGTTTGTATTTATCGTTAATTGTATTTAGTGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACTGAGC‑TGCCATGCTcc                                                                                                                                         <  1:664384/150‑1 (MQ=255)
           ttACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTcc                                                                                                                                         >  2:1494096/1‑150 (MQ=255)
              cGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCAATGCTCCGtt                                                                                                                                      >  1:1631263/1‑150 (MQ=255)
               gCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGttt                                                                                                                                     <  2:1175397/150‑1 (MQ=255)
                ccAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTTCTGATGGCATTGGGGATTGATTTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTa                                                                                                                                    >  2:1759611/1‑150 (MQ=255)
                 cAGTCCGGCATCAACGTGATTACGCTTATTATGGCGCTATAGCTGATGGCATTGTGGATTGAATTTATCGTTACTGTTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCATGGCA‑GCTCCGtttca                                                                                                                                   <  2:1223518/150‑3 (MQ=255)
                   gACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCATGGCA‑GCTCCGtttcatt                                                                                                                                 <  2:775521/150‑5 (MQ=255)
                     ccGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCg                                                                                                                               >  1:1030207/1‑150 (MQ=255)
                      cGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGTGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGt                                                                                                                              <  1:1124450/150‑1 (MQ=255)
                        gCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCa                                                                                                                            <  1:1198403/150‑1 (MQ=255)
                         cATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATACCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCat                                                                                                                           >  2:837506/1‑150 (MQ=255)
                          aTCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCatt                                                                                                                          >  2:290146/1‑150 (MQ=255)
                          aTCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCAACTGACTGCTTAATTAATTCCTTTCAAATGAAAAGGAGC‑TGCAATGCTCCGTATACTTCGTCatt                                                                                                                          >  1:171826/1‑150 (MQ=255)
                            cAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCattat                                                                                                                        >  2:1653123/1‑150 (MQ=255)
                             aaCGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCattatt                                                                                                                       >  1:1905024/1‑150 (MQ=255)
                                gTGCTTACGCGTTATATGGGGCTATTGCTGAGGGCCTTGGGGAATGAATTTATCGTTACCCTTATAAAGGGTATTTTCCCCGGCCTGCTTAATTAACTCCTTTCAAATTAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATAAttttt                                                                                                                    <  2:367980/150‑1 (MQ=255)
                                 tGATTACGCGTATTATGGGGATATTGCTGATGGAATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTa                                                                                                                   <  2:813871/150‑1 (MQ=255)
                                  gATTACGCGTATTATGGGGCTATTTCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTAc                                                                                                                  <  1:8036/150‑1 (MQ=255)
                                      aCGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTg                                                                                                              <  2:1600483/150‑1 (MQ=255)
                                       cgcgTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGt                                                                                                             <  1:637635/150‑1 (MQ=255)
                                              tatGGGGCTATTGCTGATGGCATTGGGGTTTGAATTTATGGTTACTTGTATTAAGGGGATTTTCCGCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCcgcg                                                                                                      <  2:1098644/150‑1 (MQ=255)
                                               atGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCgcgc                                                                                                     <  2:232999/150‑1 (MQ=255)
                                                      tATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATc                                                                                              <  2:1957130/150‑1 (MQ=255)
                                                      tATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATc                                                                                              <  2:70373/150‑1 (MQ=255)
                                                        ttGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATcaa                                                                                            >  2:1431391/1‑150 (MQ=255)
                                                           cTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATcaaaag                                                                                         <  2:1003972/150‑1 (MQ=255)
                                                             gATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATcaaaagca                                                                                       <  1:626928/150‑1 (MQ=255)
                                                              aTGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATcaaaagcaa                                                                                      <  2:1792101/150‑1 (MQ=255)
                                                              aTGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATcaaaagcaa                                                                                      >  1:1074648/1‑150 (MQ=255)
                                                              aTGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAAGTAATTCCTTTCAAATGAAACTGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATcaaaagcaa                                                                                      >  1:86405/1‑150 (MQ=255)
                                                                  cATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAgg                                                                                  <  1:590833/150‑1 (MQ=255)
                                                                    ttGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGaa                                                                                >  1:1917186/1‑150 (MQ=255)
                                                                      ggggATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATa                                                                              <  2:380098/150‑1 (MQ=255)
                                                                      ggggATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATa                                                                              <  1:161176/150‑1 (MQ=255)
                                                                      ggggATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTACTTTCAAAAGAAACGCAGC‑AGCCATGCTCCGATTACTTCGTCATTATTTTTACTATGTTCCCGCGCAGTTATCAAAACCAAAAGTAATa                                                                              >  2:779430/1‑150 (MQ=255)
                                                                        ggATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCGTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAgg                                                                            >  1:1068673/1‑150 (MQ=255)
                                                                        ggATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAgg                                                                            >  1:245702/1‑150 (MQ=255)
                                                                          aTTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTTTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTa                                                                          >  2:580649/1‑150 (MQ=255)
                                                                                  tATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTTTCAAAAGCAAAAGGAATAGGTAAAAATAtt                                                                  >  1:71586/1‑150 (MQ=255)
                                                                                  tATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATAtt                                                                  <  2:138010/150‑1 (MQ=255)
                                                                                   aTCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATAttc                                                                 <  1:440779/150‑1 (MQ=255)
                                                                                   aTCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATAttc                                                                 <  1:414916/150‑1 (MQ=255)
                                                                                      gTTACTGGTATTAAGGGGATTTTCCCCGACCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATAttcttc                                                              >  2:1052856/1‑150 (MQ=255)
                                                                                      gTTACTGGTATTAAGGGGATTTTCACCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATAttcttc                                                              <  1:1290294/150‑1 (MQ=255)
                                                                                         aCTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAATGGAATAGTTAAAAATATTCTTCTCa                                                           >  1:1235384/1‑150 (MQ=255)
                                                                                          cTGGTATTAAGGGTATTTTCCCCGG‑CTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCaaa                                                         <  2:873511/150‑1 (MQ=255)
                                                                                          cTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCaa                                                          >  2:1286677/1‑150 (MQ=255)
                                                                                            ggTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCAATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGAAGTTATCAAAAGCAAAAGGAATAGGTAATAATATTCTTCTCAAAt                                                        >  1:1265201/1‑150 (MQ=255)
                                                                                             gTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAAtt                                                       <  1:1494096/150‑1 (MQ=255)
                                                                                                 tAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACag                                                   <  2:1971374/150‑1 (MQ=255)
                                                                                                      ggATTTTCCCCGGCCTGCTTTATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTTCTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACagttagt                                              <  1:1770296/150‑1 (MQ=255)
                                                                                                      ggATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTTAAAATATTCTTCTCAAATTACagttagt                                              >  1:1663034/1‑150 (MQ=255)
                                                                                                        aTTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAgttagtta                                            <  2:541211/150‑1 (MQ=255)
                                                                                                         ttttCCCCGGCCTGCTTAATTTATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTtat                                           >  1:323270/1‑150 (MQ=255)
                                                                                                             ccccGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAgg                                       >  2:290336/1‑150 (MQ=255)
                                                                                                              cccGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGa                                      >  2:469739/1‑150 (MQ=255)
                                                                                                              cccGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCAGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGa                                      >  2:1295355/1‑150 (MQ=255)
                                                                                                                 ggCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGAttt                                   <  1:1973116/150‑1 (MQ=255)
                                                                                                                  gCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTc                                  >  1:67987/1‑150 (MQ=255)
                                                                                                                  gCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTc                                  <  1:574035/150‑1 (MQ=255)
                                                                                                                  gCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTc                                  >  1:468905/1‑150 (MQ=255)
                                                                                                                   ccTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTcc                                 <  1:1689795/150‑1 (MQ=255)
                                                                                                                     tGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCtt                               <  2:602683/150‑1 (MQ=255)
                                                                                                                     tGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCtt                               <  1:180282/150‑1 (MQ=255)
                                                                                                                      tcTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAACAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTa                              <  1:301567/149‑1 (MQ=255)
                                                                                                                      gCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTa                              >  1:1546005/1‑150 (MQ=255)
                                                                                                                      gCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTa                              <  1:1808291/150‑1 (MQ=255)
                                                                                                                         taattaatTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACt                           <  1:1204690/150‑1 (MQ=255)
                                                                                                                         taattaatTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACt                           <  2:1271960/150‑1 (MQ=255)
                                                                                                                          aattaattCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTg                          <  2:1239638/150‑1 (MQ=255)
                                                                                                                            ttcattCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCt                        <  2:680374/147‑1 (MQ=255)
                                                                                                                              aattCCTTTCAAATGAAACGGAGC‑TGCAATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTtc                      <  1:1013768/150‑1 (MQ=255)
                                                                                                                               attCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTtct                     <  1:1801423/150‑1 (MQ=255)
                                                                                                                               attCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTtct                     <  1:1689136/150‑1 (MQ=255)
                                                                                                                               attCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTtct                     >  2:600741/1‑150 (MQ=255)
                                                                                                                               attCCTATCAAATGAAACGGAGC‑TGCCAAGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTtct                     <  1:192331/150‑1 (MQ=255)
                                                                                                                                   cTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTctcctc                 <  2:1905024/150‑1 (MQ=255)
                                                                                                                                   cTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTctcctc                 <  1:1653123/150‑1 (MQ=255)
                                                                                                                                   cTTTCAAATGAAACGGAGC‑TGCAATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAGTTACAGTTTGTTATATGGTTTTCCTTAACTGCTTctcctc                 >  1:269258/1‑150 (MQ=255)
                                                                                                                                    tttCAAATGACACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCa                <  1:1085829/150‑1 (MQ=255)
                                                                                                                                    tttCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCa                <  2:1817101/150‑1 (MQ=255)
                                                                                                                                         aaTGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACcatc           >  2:483784/1‑150 (MQ=255)
                                                                                                                                          agtaaaCGGAGCATGGCA‑GCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCATCACcaaca          >  1:128177/4‑150 (MQ=255)
                                                                                                                                          aTGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACcatca          <  1:332277/150‑1 (MQ=255)
                                                                                                                                          aTGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACcatca          >  2:590492/1‑150 (MQ=255)
                                                                                                                                          aTGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACcatca          <  1:578246/150‑1 (MQ=255)
                                                                                                                                            taaaCGGAGCATGGCA‑GCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATg        >  1:1111561/2‑150 (MQ=255)
                                                                                                                                            gAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATg        >  2:1963697/1‑150 (MQ=255)
                                                                                                                                             aaaCGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGt       >  1:1381672/1‑150 (MQ=255)
                                                                                                                                               aCGGAGC‑TGCCCTGCTCCGTTTACTTCGGCATTATTTTTACTTTGTTCCCGCGCAGGTATCAAAAGCAAAAGGAAGTGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTa     >  2:1515107/1‑150 (MQ=255)
                                                                                                                                               aCGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTa     <  2:1631263/150‑1 (MQ=255)
                                                                                                                                               aCGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAAAATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTa     <  2:1622882/150‑1 (MQ=255)
                                                                                                                                                 ggAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATTTTAtt   >  2:771026/1‑150 (MQ=255)
                                                                                                                                                 ggAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTAATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTCCAGTTAGTTTTAAGGATTTACTTAACTGATTCTCCTCACCATCATGTTAtt   >  1:1233676/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
TGGTACGGGTTTTACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTT  >  NC_000913/1299118‑1299413

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