Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913501,5750AC44.4% 124.3 / 14.8 72V555G (GTG→GGG) ybaLinner membrane putative NAD(P)‑binding transporter
Reads supporting (aligned to +/- strand):  ref base A (35/5);  new base C (0/32);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 1.55e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTT  >  NC_000913/501428‑501724
                                                                                                                                                    |                                                                                                                                                     
aaGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCt                                                                                                                                                    <  1:450591/151‑1 (MQ=255)
         tcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACccgcc                                                                                                                                            <  1:1117573/150‑1 (MQ=255)
           ggTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCgg                                                                                                                                          <  2:1237923/150‑1 (MQ=255)
                   ggCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTc                                                                                                                                 >  1:421622/1‑151 (MQ=255)
                     ccAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTcc                                                                                                                                >  1:884091/1‑150 (MQ=255)
                      cAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTcc                                                                                                                                >  2:90871/1‑149 (MQ=255)
                       aGATAGGGCGTTCCCGCCGCATTTGGTTTTTGGCTCTCGATCCCTGATGCGCCGCTGGCGCGTCTTATCATCCCTACATTTTTTTCATATTTTCCACCCGGAACCCCCCTTTTC‑CCCCGTCCCCCCCCCACCCGCCGGTGGCGTTTcc                                                                                                                                <  1:1119601/148‑1 (MQ=255)
                       aGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagc                                                                                                                             >  1:1094251/1‑151 (MQ=255)
                        gATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagc                                                                                                                             >  2:1105609/1‑150 (MQ=255)
                          tAAGGCTTCCCCCCCTCACCTGTCATTTGGCTCTCGATGCCTGATGCCCCGCGGGCGCGTCTTATCATGCCTACATTTTTTTCTTATTTTCCACCCGGCACCCCCCTTTTC‑CCCCGCCCCCCCCACCCCCGCCGGTGGCGTTTCca                                                                                                                               <  1:683391/146‑1 (MQ=255)
                          tAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGt                                                                                                                          >  2:603859/1‑151 (MQ=255)
                          tAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGt                                                                                                                          >  1:70047/1‑151 (MQ=255)
                           aaGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGtt                                                                                                                         >  1:700243/1‑151 (MQ=255)
                            aGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTc                                                                                                                        >  2:587561/1‑151 (MQ=255)
                                     cgccgcTTCGGGCTTTTGCCTCTCGATGCCTGATGCGACGCTGGCGCTTTTTATCATCCCTACATTTTTTTCATTTTTTACCCCCGCCACCCCCCGTTTC‑CCCCGCCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCAt                                                                                                                  <  1:296140/148‑1 (MQ=255)
                                     cgccccATCTGGATTTTGGCTCTCGATCCCTGATGCGCCGCGGGCGCTTTTTTTCATCCCTACATATTTTTCATTTTTTACACCCGCCACCCCCCTTTTC‑CCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGt                                                                                                               <  1:178236/151‑1 (MQ=255)
                                       ccgcATCTGGCTTTTGGCTCTCGATGCCTGAGGCGACGCTGGCGCTTTTTATCATGCCTAATTTTTTTTCATATTTTACACCCGCCACCCCCCTTTTC‑CCCCTCCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATgg                                                                                                                <  2:69116/148‑1 (MQ=255)
                                       ccgcATCTGGCTTTTGCCTCTCGATGCCTGTGGCAACGCGGGCGCTTTTTATCATCCCTACATTTTTTTCATTTTTTACCTCCGGCACCCCCTGTTTC‑CCCCTTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGt                                                                                                               <  1:382616/149‑1 (MQ=255)
                                       ccgcATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTAc                                                                                                             >  2:464999/1‑151 (MQ=255)
                                        cgcATCTGGTTTTTGGCTCTCGATCCCTGATGCGCCGCGGGCGCTTTTTATCATCCCTAATTTTTTTTCATTTTTTCCACCGGGCACCCCCCGTTTC‑CCCCTCCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTa                                                                                                              <  1:782838/149‑1 (MQ=255)
                                        cccttctgctttttGGCTCTCGATGCCTGATGCGCCGCGGGCGCGTTTTTTCATCCCTACATATTTTTCATATTTTACCCCCGGCACCCCCCTTTTC‑CCCCTCCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGc                                                                                                                    <  1:477340/132‑1 (MQ=255)
                                          cATCTGGCTTTTGGCTCTCGATGCCTGAGGCGACGCTGGCGCTTCTTATCATCCCTACATATTTTTCATATTTTACATCCGGCAACCCCCGTTTC‑CCCCGTCCCCACCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACgg                                                                                                           <  2:924184/150‑1 (MQ=255)
                                             cTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCg                                                                                                        >  1:668936/1‑150 (MQ=255)
                                                   tttGGCTCTCGATGCCTGATGCGACGCGGGGGCTTCTTTTCATGCTTACTTTTTTTTCATTTTTTACACCCGGCACCCACCTTTTC‑CCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAttt                                                                                                    <  1:981847/148‑1 (MQ=255)
                                                     tGGCTCTCGATGCCTGATGCGCCGCTGGCGCGTTTTATCATGCCTACATTTTTTTCATATTTTACACCCGGCAACCACCGTTTC‑CCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACg                                                                                                <  1:424366/150‑1 (MQ=255)
                                                        ctcccgagcccTGTTGCGCCGCGGGCGCTTTTTATCATCCCTACATTTTTTTCATTTTTTACACCCGCCACCCCCCTTTTC‑CCCCGTCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCt                                                                                              <  1:520281/140‑1 (MQ=255)
                                                            cGATCCCTGAGGCGACGCGGGCGCTTTTTATCATCCCTACATTTTTTTCATTTTTTACACCCGGCACCCCCCGTTTC‑CCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc                                                                                           <  1:1050816/148‑1 (MQ=255)
                                                                  cTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACt                                                                                   >  2:770464/1‑150 (MQ=255)
                                                                    gATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAcc                                                                                >  2:69399/1‑151 (MQ=255)
                                                                     aTGCGCCGCTGGCGCTTCTTATCATCCCTACATATTTTCCATTTTTTCCACCCGGCACCCCCCTTTTC‑CCCCTCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACt                                                                                   <  2:1042315/147‑1 (MQ=255)
                                                                         gccGCTGGCGCTTTTTATCTTCCCTACATTTTTTTCTTATTTTCCATCCGGCAACCCCCGTTTC‑CCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAt                                                                            <  2:1142080/148‑1 (MQ=255)
                                                                          aCGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTc                                                                          >  1:661883/1‑151 (MQ=255)
                                                                          aCGCTGGCGCGTCTTATCATCCCTACATATTTTTCATTTTTTACACCCGGCAACCCCCGTTTC‑CCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAtt                                                                           <  1:699131/150‑1 (MQ=255)
                                                                              tGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcac                                                                      >  1:1026880/1‑151 (MQ=255)
                                                                               ggCGCGTTTTATCATGCCTACATTTTTTTCATTTTTTCCCCCCGGCACCCCCCTTTTC‑CCCCGCCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGca                                                                       <  2:228275/149‑1 (MQ=255)
                                                                                 cgcgTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCTCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcaccac                                                                   >  2:833707/1‑151 (MQ=255)
                                                                                  gcgTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACg                                                                  >  1:871390/1‑151 (MQ=255)
                                                                                   cgTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGt                                                                 >  1:425747/1‑151 (MQ=255)
                                                                                    gTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt                                                                >  1:1249296/1‑151 (MQ=255)
                                                                                      cTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCg                                                              >  2:37797/1‑151 (MQ=255)
                                                                                            aTGCCTACATTTTTTTCATATTTTACACCCGGCAACCCCCTTTTC‑CCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGat                                                         <  2:120253/150‑1 (MQ=255)
                                                                                            aTGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGat                                                         >  2:573737/1‑150 (MQ=255)
                                                                                              gCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGatat                                                       >  2:561460/1‑150 (MQ=255)
                                                                                               ccTACATTTTTTCCATATTTTACACCCGCAAACCCCCTTTTC‑CCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGa                                                          <  1:1170347/146‑1 (MQ=255)
                                                                                                  acttttttttCATATTTTCCACCCGCCAACCCCCTTTTC‑CCCCGCCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCgg                                                             <  1:1025645/135‑1 (MQ=255)
                                                                                                     tttttttCATTTTTTACACCCGGCACCCACCTTTTC‑CCCCGCCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCAc                                                <  2:255181/148‑1 (MQ=255)
                                                                                                     tatTTTTCATATTTTACCCCCGCCACCCCCCGTTTC‑CCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCAc                                                <  1:768635/150‑1 (MQ=255)
                                                                                                      atTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACtt                                              >  1:721295/1‑151 (MQ=255)
                                                                                                       tttttCATATTTTACACCCGGCAACCACCGTTTC‑CCCCCCAACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACt                                               <  2:87692/149‑1 (MQ=255)
                                                                                                       tttttCATATTTTACACCCGGCAACCACCGTTTA‑CCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACt                                               <  2:23877/149‑1 (MQ=255)
                                                                                                         ttttATTTTTTCCATCCGCCACCCCCCGTTCC‑CCCCCCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCa                                            <  1:750223/150‑1 (MQ=255)
                                                                                                             attttttACACCCGGCACCCCCCGTTTC‑CCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATc                                          <  1:748104/145‑1 (MQ=255)
                                                                                                               atTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCAt                                     >  2:1091537/1‑151 (MQ=255)
                                                                                                                ttttACACCCGCCACCCCCCGTTTC‑CCCCTCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCAt                                     <  1:1177645/150‑1 (MQ=255)
                                                                                                                 tttACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATa                                    >  2:689990/1‑150 (MQ=255)
                                                                                                                 tttACACCCGCCACCCCCCTTTTC‑CCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATa                                    <  1:630519/150‑1 (MQ=255)
                                                                                                                     cATCCGGCACCCCCCTTTTC‑CCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAAt                                  <  1:409689/148‑1 (MQ=255)
                                                                                                                     cATCCGGCACCCCCCGTTTC‑CCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCACACCACGTTCGGTGATATACGCCACTTCATCGTCATAATgg                                <  2:884104/150‑1 (MQ=255)
                                                                                                                     cATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATggg                               >  1:908784/1‑151 (MQ=255)
                                                                                                                      cacccgCAACCCCCCTTTCCCCCCGTCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAAt                                  <  2:224459/143‑1 (MQ=255)
                                                                                                                      aTCCGGCACCCCCCGTTTC‑CCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAAt                                  <  2:37781/147‑1 (MQ=255)
                                                                                                                         cGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGcgcg                           >  2:605062/1‑151 (MQ=255)
                                                                                                                          ggCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCgg                          >  2:1246497/1‑151 (MQ=255)
                                                                                                                               cccccGTTTC‑CCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGAGCaa                      <  1:11844/147‑1 (MQ=255)
                                                                                                                                caccGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaat                     >  2:41539/1‑150 (MQ=255)
                                                                                                                                caccGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaat                     >  1:1128005/1‑150 (MQ=255)
                                                                                                                                  ccGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaataat                  >  2:1091566/1‑151 (MQ=255)
                                                                                                                                    gTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAtct                >  2:503229/1‑151 (MQ=255)
                                                                                                                                          ccccGTCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCa              <  2:465134/148‑1 (MQ=255)
                                                                                                                                          ccccGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAat            >  1:1138622/1‑150 (MQ=255)
                                                                                                                                          ccccGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAat            >  1:647006/1‑150 (MQ=255)
                                                                                                                                            ccGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATc         >  1:1018650/1‑151 (MQ=255)
                                                                                                                                                 accaccTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGAt    >  2:85335/1‑151 (MQ=255)
                                                                                                                                                   caccTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGAttt  <  1:246957/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTT  >  NC_000913/501428‑501724

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