Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913501,5720AC46.3% 95.9 / 11.0 67V556G (GTG→GGG) ybaLinner membrane putative NAD(P)‑binding transporter
Reads supporting (aligned to +/- strand):  ref base A (32/4);  new base C (1/30);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 1.44e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCG  >  NC_000913/501423‑501719
                                                                                                                                                       |                                                                                                                                                   
tCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATCTTTTTCATATTTTACATCCGGCAA‑CCATCGTTT‑ACCCCGTca                                                                                                                                                     >  2:135620/1‑150 (MQ=255)
tCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTca                                                                                                                                                     >  1:289505/1‑150 (MQ=255)
tCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTca                                                                                                                                                     >  2:444371/1‑150 (MQ=255)
 cGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTcacc                                                                                                                                                   <  2:1086672/151‑1 (MQ=255)
     aaGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTcaccac                                                                                                                                                 >  1:797016/1‑149 (MQ=255)
       gTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCATCATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCTTCCCCACCTCa                                                                                                                                             >  2:820865/1‑151 (MQ=255)
                 gttggtggcccAGATAGGGCGTCCCCGCCGCATCTGCCTTTTGGCTCTCGATGCCTGAGGCGACGCTGGCGCGTCTTATCATGCCTCCATATTTTTCATATTTTCCATCCGGCAA‑CCCCCTTTT‑CCCCCGTCCCCCCCTCCCCCGCCGGt                                                                                                                                    <  2:803778/141‑1 (MQ=255)
                 gTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTg                                                                                                                                   >  2:117440/1‑151 (MQ=255)
                   tCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGACACCACCTCACCCGCCGGTgg                                                                                                                                  >  2:452343/1‑150 (MQ=255)
                             gATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCccgca                                                                                                                       >  2:940437/1‑151 (MQ=255)
                             gATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagca                                                                                                                       >  2:561995/1‑151 (MQ=255)
                             gATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagca                                                                                                                       >  1:448337/1‑151 (MQ=255)
                                  ggCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTc                                                                                                                   >  1:676203/1‑150 (MQ=255)
                                     gtccccgccgcTTCGGGTTTTTGGCTCTCGATGCTTGATGCGCCGCGGGCGCGTTTTATCATCCCTACATTTTTTTCATTTTTTACATCCGGCAC‑CCCCCGTTT‑CCCCCGTCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAg                                                                                                                <  1:478710/145‑1 (MQ=255)
                                          cgccgcATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGt                                                                                                          >  2:667466/1‑151 (MQ=255)
                                          cgccgcATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGt                                                                                                          >  1:188777/1‑151 (MQ=255)
                                              gcATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACgg                                                                                                      >  2:160726/1‑151 (MQ=255)
                                                    ggCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAtt                                                                                                >  1:260275/1‑151 (MQ=255)
                                                    ggCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAtt                                                                                                >  1:758573/1‑151 (MQ=255)
                                                      cttttGCCTCTCGATGCCTGATGCGACGCTGGCCCTTCTTAACAGCCCTACATTTTTTTCATATTTTACATCCGGCAA‑CCCCCTTTT‑CCCCCCCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTc                                                                                              <  1:172683/149‑1 (MQ=255)
                                                       ttttGCCTCTCGTTGCCTGATGCGCCGCGGGCGCGTCTTATCATCCCTACATTTTTTTAATATTTTACATCCGGCAC‑CCACCTTTT‑CCCCCGTCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTc                                                                                              <  2:648759/149‑1 (MQ=255)
                                                          tgtctctGGATGCCGGTGGGGCCGCGGGCGCGTCTTATCTTCCCTACATTTTTTTCATTTTTTACTCCCGGAAC‑CCCCCTTTT‑CCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGc                                                                                          <  2:835375/148‑1 (MQ=255)
                                                           ggCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACg                                                                                           >  1:1041234/1‑149 (MQ=255)
                                                            gCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCt                                                                                         >  2:25226/1‑150 (MQ=255)
                                                               ctcGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc                                                                                      >  1:1134267/1‑150 (MQ=255)
                                                               ctcGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc                                                                                      >  1:842770/1‑150 (MQ=255)
                                                                     gCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCa                                                                                >  2:242764/1‑150 (MQ=255)
                                                                      cccgaggcgcCGCTGGCGCGTCTTCTCATCCCCCCTTTTTTTTAATATTTCCCCCCCGGCAC‑CCCCCTTTT‑CCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACt                                                                              <  1:824348/145‑1 (MQ=255)
                                                                      ccTGATGCGCCGCGGGCGCTTTTTTTCATCCCTACATTTTTTTCATTTTTTCCACCCGGCAA‑CCCCCGTTT‑CCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACt                                                                              <  1:117464/151‑1 (MQ=255)
                                                                        tggtgcgccgggggCGCTTCTTCTCTTGCCTACATTTTTTCCATATTTTCCACCCGGCAA‑CCCCCGTTT‑CCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAc                                                                            <  1:205877/139‑1 (MQ=255)
                                                                          atgccgccctcgcgcTTCTTCTCATCCCTATATATTTTTCATATTTTCCATCCGGCAC‑CCCCCTTTC‑CCCCCCCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTGcc                                                                           <  2:212457/139‑1 (MQ=255)
                                                                           tGCGCCGCGGGCGCGTCTTATCATGCCTAAATTTTTTTCATATTTTACACCCGGCAC‑CCCCCTTTT‑CCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTg                                                                         <  1:138460/151‑1 (MQ=255)
                                                                             cgccgcgGGCGCGTTTTATCATGCCTACATTTTTTTCATATTTTACACCCGGCAA‑CCCCCGTTT‑CCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAt                                                                       <  1:566720/148‑1 (MQ=255)
                                                                               ccGCTGGCGCGTCTTATCATGCCTACATTTTTTTCATATTTTACATCCGGCAC‑CCCCCGTTT‑CCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAtt                                                                      <  1:744571/149‑1 (MQ=255)
                                                                                         tttttATCACGCCTCCTTTTTTTTCATATTTTCCACCCGCAAC‑CCACCTTTC‑CCCCCGTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt                                                           <  2:512877/148‑1 (MQ=255)
                                                                                         ttCTTATCATGCCTACATATTTTTCATTTTTTACACCCGGCAA‑CCACCTTTT‑CCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTc                                                                     <  2:328610/140‑1 (MQ=255)
                                                                                         gTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTc                                                                     >  1:328686/1‑141 (MQ=255)
                                                                                         gTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt                                                           >  2:229466/1‑151 (MQ=255)
                                                                                         gTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGt                                                            >  2:76557/1‑150 (MQ=255)
                                                                                             tATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGt                                                       >  2:133014/1‑151 (MQ=255)
                                                                                               tCATCCCTAAATTTTTTTCATTTTTACCACCCGGCAC‑CCCCCTTTT‑CCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGa                                                     <  1:954023/151‑1 (MQ=255)
                                                                                                  tccctcattttttttAATTTTTTCCACCCGGCAC‑CCCCCTTTT‑CCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGatat                                                  <  1:935397/141‑1 (MQ=255)
                                                                                                  tccccaCATTTTTTTCATTTTTCCCCCCCGGCCA‑CCCCCTTTT‑CCCCCCCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGatat                                                  <  1:433987/146‑1 (MQ=255)
                                                                                                  tGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGatat                                                  >  1:704576/1‑151 (MQ=255)
                                                                                                   gCCTACATATTTTTAATATTTTACATCCTCCAA‑CCACCGTTT‑ACACCGTCACCACCTCACCAGCCGGTGGCGTTTCCAGCAGTTACAGCATGGTACCGGCGATTTCACGCTCGTCCATCACTACCTGATTCGCACCACGATCGGTGAtata                                                 >  2:1162543/1‑151 (MQ=255)
                                                                                                    ccTACATTTTTTTCATTTTTTCCATCCGGCAA‑CCCCCGTTT‑CCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                                                <  1:797248/151‑1 (MQ=255)
                                                                                                    ccTACATTTTTTTCATATTTTACATCCGGCAC‑CCCCCTTTT‑CCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                                                <  1:541071/151‑1 (MQ=255)
                                                                                                    ccTACATTTTTTTCATATTTTACACCCGCCCC‑CCCCCGTTC‑CCCCCCTCCCCCCCTCTCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                                                <  2:842436/151‑1 (MQ=255)
                                                                                                    ccTACATATTTTTCATATTTTACATCCGGCAA‑CCACCTTTT‑CCCCCGTCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                                                <  2:482663/151‑1 (MQ=255)
                                                                                                    ccTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGAtata                                                 >  2:648545/1‑150 (MQ=255)
                                                                                                     ctccATATTTTTCATATTTTACACCCGGCAC‑CCCCCGTTT‑CCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACg                                               <  1:775573/148‑1 (MQ=255)
                                                                                                     cTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACg                                               >  1:479246/1‑151 (MQ=255)
                                                                                                      cactttttttcCATTTTTTCCCTCGGGAAC‑CCACCTTTT‑ACCCCCCCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGc                                              <  1:964200/145‑1 (MQ=255)
                                                                                                      cactttttttacttttttcACAGCCGGCAA‑CCCCCTTTT‑CCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGc                                              <  2:1007387/136‑1 (MQ=255)
                                                                                                       aCATTTTTTTCATTTTTTACACCCGGCAC‑CCCCCGTTT‑CCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGcc                                             <  1:977999/151‑1 (MQ=255)
                                                                                                        cATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCa                                            >  2:449704/1‑151 (MQ=255)
                                                                                                              tttCATTTTTTCCCCCCGGCAA‑CCCCCTTTT‑CCCCCCCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCAt                                      <  1:532024/151‑1 (MQ=255)
                                                                                                                taattttttccccccGGCAA‑CCCCCGTTT‑CCCCCGCCCCCACCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCg                                    <  2:1200857/138‑1 (MQ=255)
                                                                                                                   tatTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCa                                 >  1:135666/1‑151 (MQ=255)
                                                                                                                        tacccccGGCAC‑CCCCCTTTT‑CCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTCCGGTGATATACGCCACTTCATCGTCATAATg                            <  1:76575/146‑1 (MQ=255)
                                                                                                                        tACATCCGGCAC‑CCCCCTTTT‑CCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATg                            <  1:449850/151‑1 (MQ=255)
                                                                                                                         cacaccGGCAACCCCCCGTTT‑CCCCCGCCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATg                            <  2:164608/147‑1 (MQ=255)
                                                                                                                          cATCCGGCAC‑CCCCCTTTT‑CCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATggg                          <  1:260942/151‑1 (MQ=255)
                                                                                                                             cccgCAA‑CCCCCCTTTCACCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcg                        <  1:783925/148‑1 (MQ=255)
                                                                                                                                        ccGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCCGCATGGTACGGGCGATTTCACGCTCGCCCGTCCCTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaatcat             >  2:98833/1‑151 (MQ=255)
                                                                                                                                        ccGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaataat             >  1:675901/1‑151 (MQ=255)
                                                                                                                                        ccGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaataat             >  2:395425/1‑151 (MQ=255)
                                                                                                                                        ccGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaataat             >  2:213745/1‑151 (MQ=255)
                                                                                                                                         ctttt‑CCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAtc            <  2:307364/149‑1 (MQ=255)
                                                                                                                                          tttt‑CCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAtct           <  1:285956/150‑1 (MQ=255)
                                                                                                                                                ccccGTCACCACCTCACCCGCCGGCGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCCTCATAATGGGCGCGGGCAATAATCTCAata      >  1:58340/1‑151 (MQ=255)
                                                                                                                                                    gTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCg  >  2:970643/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                   
TCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCG  >  NC_000913/501423‑501719

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