Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913501,5750AC47.0% 106.1 / 11.2 66V555G (GTG→GGG) ybaLinner membrane putative NAD(P)‑binding transporter
Reads supporting (aligned to +/- strand):  ref base A (33/2);  new base C (0/31);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.65e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCG‑TCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGG  >  NC_000913/501426‑501720
                                                                                                                                                      |                                                                                                                                                 
ttAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCG‑Tcaccac                                                                                                                                                  >  2:553310/1‑151 (MQ=255)
         cgtcgGTTCGTAGGCCAGATAAGGCGTTCGCGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCG‑TCACCACCTCACccgc                                                                                                                                         <  2:688871/151‑1 (MQ=255)
                tCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCG‑TCACCACCTCACCCGCCGGTgg                                                                                                                                   >  1:706235/1‑150 (MQ=255)
                  gTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCG‑TCACCACCTCACCCGCCGGTGGCGt                                                                                                                                >  2:1045770/1‑151 (MQ=255)
                    aGGCCAGATAAGGCGTTCGCGCCGCTTCTGGCATTTGGCTCTCGATGCCTGATGCGGCGCTGGCGCGTTTTATCTTGCCACCATATTTTTCATTTTTCACGTCCGGCACCCTCCTTTTACCCCG‑TCACCCCCTCACCCGCCGGTggt                                                                                                                                  <  2:50932/147‑2 (MQ=255)
                                     cACGCCGCATCTGGCTTTTGGCTCTCGTTGCCGGATGCGCCGCTGGCGCGTCTTATCATGCCTACATTTTTTTCATTTTTTCCATCCGGCACCCCCCGTTTCCCCCG‑TCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCa                                                                                                               <  1:777744/149‑1 (MQ=255)
                                      aCGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCG‑TCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATgg                                                                                                            >  1:873875/1‑151 (MQ=255)
                                       cgccgcATCTGGCTTTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCCCCGTTTCCCCCG‑TCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGt                                                                                                           <  1:346083/151‑1 (MQ=255)
                                          cgcATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCG‑TCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACg                                                                                                        >  2:808034/1‑151 (MQ=255)
                                           gcATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCG‑TCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACgg                                                                                                       >  2:149125/1‑151 (MQ=255)
                                                tGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCG‑TCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAt                                                                                                  >  2:763099/1‑151 (MQ=255)
                                                           tctcGATGCGTGAGGCGGCGCTGGCGCGTCTTATCATCCCAACATATTTTTCATATTTTCCATCCGGCAACCACCGTTTCCCCCG‑TCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACg                                                                                            <  1:811102/146‑1 (MQ=255)
                                                           tctcGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCG‑TCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc                                                                                       >  2:243956/1‑151 (MQ=255)
                                                             tcGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCG‑TCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGccc                                                                                     >  2:68223/1‑151 (MQ=255)
                                                              cGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACACCCGGCACCCCCCGTTTCCCCCG‑TCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCa                                                                                    <  1:553466/151‑1 (MQ=38)
                                                               gATGCCGGATGCGCCGCTGGCGCGTTTTATCATGCCTACATATTTTTCATTTTTTCCATCCGGCAACCCCCGTTTCCCCCG‑TCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCAt                                                                                   <  1:1167571/151‑1 (MQ=255)
                                                                   ccTGATGCGACGCTGGCGCGTTTTATCTTGCCTATATTTTTTAATTATTTTACACCCGGCACCCACCGTTCCCCCCC‑TCCCCACCACCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACt                                                                               <  1:641606/151‑1 (MQ=255)
                                                                   ccTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCG‑TCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACt                                                                               >  2:842310/1‑151 (MQ=255)
                                                                   ccTGAGGCGACGCGGGCGCGTCTTATCATGCCTACATTTTTTCAATATTTTCCACCCGGCACCCCCCTTTCCCCCCG‑TCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATc                                                                                  <  2:1064282/148‑1 (MQ=255)
                                                                      gATGCGACGCGGGCGCTTTTAATCATCCCTACATTTTTTTCTTTTTTTACACCCGGCACCCCCCTTTTCCCCCG‑CCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAcc                                                                            <  2:344564/151‑1 (MQ=255)
                                                                        tGCGCCGCTGGCGCGTTTTATCATCCCTACATATTTTTCATATTTTACACCCGCCACCCCCCGTTTCCCCCG‑TCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCt                                                                           <  2:27890/150‑1 (MQ=255)
                                                                          cgccgcgGGCGCGTTTTTTCATCCCTACATTTTTTTCATTTTTTACACCCGCCTCCCCCCTTTTCCCCCG‑TCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGa                                                                         <  2:543162/147‑1 (MQ=255)
                                                                           gACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCG‑TCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAtt                                                                       >  2:325837/1‑151 (MQ=255)
                                                                           gACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCG‑TCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAtt                                                                       >  1:173980/1‑151 (MQ=255)
                                                                            ccgcgGGCGCGTTTTACCATGCCTACATTTTTTTCATATTTTACACCCGCCACCCCCCGTTTCCCCCG‑CCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAcc                                                                            <  2:321517/144‑1 (MQ=255)
                                                                             cgcgGGCGCTTTTTTTCTTCCCTACATTTTTTTAATTTTTTACACCCGGCACCCCCCTTTTCCCCCG‑TCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTc                                                                      <  1:253188/150‑1 (MQ=255)
                                                                             cGCTGGCGCTTTTTATCATGCCTAAATTTTTTTCATATTTTCCCCCGGGCACCCCCCGTTTCCCCCT‑TCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTc                                                                      <  2:486676/150‑1 (MQ=255)
                                                                              gcgggCGCGTCTTATCATGCCTACATTTTTTTCATATTTTCCACCCGGCAACCCCCTTTTCCCCCG‑TCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTc                                                                      <  1:874293/146‑1 (MQ=255)
                                                                              gCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCG‑TCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGc                                                                    >  1:157094/1‑151 (MQ=255)
                                                                                  gcgcGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCG‑TCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcacc                                                                 >  1:1170706/1‑150 (MQ=255)
                                                                                     cttCTTATCATCCCTACATATTTTTCATATTTTACACCCGGCACCCCCCGTTTCCCCCG‑CCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcaccac                                                               <  1:149157/147‑1 (MQ=255)
                                                                                       ttttttcttcccTACATTTTTTTCATTTTTTACATCCGGCACCCCCCTTTTCCCCCT‑CCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGt                                                             <  2:1056425/139‑1 (MQ=255)
                                                                                        cTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCG‑TCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCg                                                          >  2:338100/1‑151 (MQ=255)
                                                                                          tATCATCCCTACATTTTTTTCATATTTAACTCCCAGAAACCCCCTTTTCCCCCT‑TTCCCCCCTTTCCCGCCGGTGGCGTTTCCAGCACTTCCAGCATGGTACGGGCGATTTCACGCTCCCCCTTCCCTACCTGATTCTCTCCACGtc                                                            <  2:214753/147‑2 (MQ=255)
                                                                                           accATGCCTACATATTTTTCATATTTTACATCCGGCACCCCCCGTTTCCCCCG‑TCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGt                                                        <  1:93422/148‑1 (MQ=255)
                                                                                                gCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCG‑TCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGatat                                                   >  1:92114/1‑150 (MQ=255)
                                                                                                 ccTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCG‑TCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                                                 >  1:645604/1‑151 (MQ=255)
                                                                                                  cTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCG‑TCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                                                 >  1:1083272/1‑150 (MQ=255)
                                                                                                      atatTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCG‑TCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCAc                                            >  1:694151/1‑151 (MQ=255)
                                                                                                       tttttttCATATTTTACACCCGGCACCCCCCGTTTCCCCCG‑CCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCAc                                            <  1:980532/148‑1 (MQ=255)
                                                                                                       tatTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCG‑TCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACt                                           >  1:225500/1‑151 (MQ=255)
                                                                                                          ttttCATATTTTCCACCCGGCACCCCCCTTTTCCCCCG‑TCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTc                                         <  1:383992/150‑1 (MQ=255)
                                                                                                          ttttCATATTTTACATCCGGCAACCACCGTTTACCCCG‑TCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTc                                         >  2:820276/1‑150 (MQ=255)
                                                                                                           tttCTTATTTTACTCCCGCCACCCACCTTTTCCCCCG‑TCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCa                                        <  1:205768/150‑1 (MQ=255)
                                                                                                                tatTTTACATCCGGCAACCACCGTTTACCCCG‑TCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTc                                   >  2:921716/1‑150 (MQ=255)
                                                                                                                 atTTTACATCCGGCAACCACCGTTTACCCCG‑TCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCAt                                 >  2:1063448/1‑151 (MQ=255)
                                                                                                                  ttttACATCCGGCAACCACCGTTTACCCCG‑TCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCAt                                 >  2:491850/1‑150 (MQ=255)
                                                                                                                     tACATCCGGCAACCACCGTTTACCCCG‑TCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATg                             >  2:370801/1‑151 (MQ=255)
                                                                                                                      aCATCCGGCAACCACCGTTTACCCCG‑TCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATgg                            >  1:82343/1‑151 (MQ=255)
                                                                                                                       cacccGGCACCCCCCTTTTCCCCCG‑TCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAAt                              <  1:1066250/145‑1 (MQ=255)
                                                                                                                        aTCCGGCAACCACCGTTTACCCCG‑TCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgc                          >  2:897659/1‑151 (MQ=255)
                                                                                                                        aTCCGGCAACCACCGTTTACCCCG‑TCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgc                          >  1:622696/1‑151 (MQ=255)
                                                                                                                         tCCGGCAACCACCGTTTACCCCG‑TCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcg                         >  2:164269/1‑151 (MQ=255)
                                                                                                                          ccGCCACCCCCCTTTTCCCCGG‑CCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcgc                        <  2:177694/151‑1 (MQ=255)
                                                                                                                          ccGACAACCCCCTTTTCCCCCG‑CCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcgc                        <  2:448556/151‑1 (MQ=255)
                                                                                                                                 cccccGTTTCCCCCG‑CCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaa                  <  1:244004/147‑1 (MQ=255)
                                                                                                                                 ccaccGTTTCCCCCG‑TCCCCCCCGCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaat                 <  1:973633/151‑1 (MQ=255)
                                                                                                                                   accGTTTACCCCG‑TCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaata                >  1:16819/1‑150 (MQ=255)
                                                                                                                                     cGTTTACCCCG‑TCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaataat              >  2:975188/1‑150 (MQ=255)
                                                                                                                                      cttttCCCCGCTCACCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACAGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATAGACGCCACTTCATCGTCATAATGGGCGCGGGc                    <  1:852096/143‑1 (MQ=255)
                                                                                                                                          cccccG‑CCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCaa         <  2:235985/149‑1 (MQ=255)
                                                                                                                                           ccccG‑CCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAata       <  1:763337/151‑1 (MQ=255)
                                                                                                                                           ccccC‑TCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCa          <  1:1046089/148‑1 (MQ=255)
                                                                                                                                            cccG‑CCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAat        <  1:646490/149‑1 (MQ=255)
                                                                                                                                             ccG‑TCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATc     >  2:1144312/1‑151 (MQ=255)
                                                                                                                                              cG‑TCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATcc    >  2:970835/1‑151 (MQ=255)
                                                                                                                                               g‑TCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCg   >  1:781095/1‑151 (MQ=255)
                                                                                                                                                  caccacCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCgg  >  1:503098/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
TTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCG‑TCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGG  >  NC_000913/501426‑501720

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: ‑