Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,155,6390TG42.2% 134.5 / 11.1 71G95G (GGT→GGG) yqhDaldehyde reductase, NADPH‑dependent
Reads supporting (aligned to +/- strand):  ref base T (7/34);  new base G (30/0);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 2.15e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAACTATCCGGAAAATATCGAT‑CCGTGGCACATTCTGCAAACGGGCGGTAAAGAGATTAAAAGCGCCATCCCGATGGGCTGTGTGCTGACGCTGCCAGCAACCGGTTCAGA  >  NC_000913/3155490‑3155788
                                                                                                                                                     |                                                                                                                                                        
gTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGtt                                                                                                                                                         <  2:299901/151‑1 (MQ=255)
     cGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCt                                                                                                                                                       <  1:1037452/148‑1 (MQ=255)
     cGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCt                                                                                                                                                       <  1:887822/148‑1 (MQ=255)
         cAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTgg                                                                                                                                                <  2:972827/151‑1 (MQ=255)
           aGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGGTCTTGACTgg                                                                                                                                                >  2:1125675/1‑149 (MQ=255)
              tCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGA‑Cgg                                                                                                                                           <  2:941603/150‑1 (MQ=255)
              tCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGGTTTGttacgggcgg                                                                                                                                            >  1:980627/1‑141 (MQ=255)
              tCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGGTCTGTAATGGG‑Cgg                                                                                                                                           >  2:1110073/1‑150 (MQ=255)
                   aTGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGGTTTTTACTGGG‑Cgggcac                                                                                                                                       >  1:902354/1‑145 (MQ=255)
                      cccTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGGGGGTTTGTACTTGG‑CGGCAccca                                                                                                                                     >  1:1081592/1‑146 (MQ=255)
                               gCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGGTCTGTACTGGG‑CGGCACCAAAATTATcgcc                                                                                                                           >  2:499521/1‑149 (MQ=255)
                                        tGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGGTCTGGAATGGG‑AGGGACCAAATTTATAGCCGCAGC‑GGCt                                                                                                                 >  1:309884/1‑149 (MQ=255)
                                         gCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTaa                                                                                                               <  2:78213/150‑1 (MQ=255)
                                            ggAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAACt                                                                                                             <  2:436426/149‑1 (MQ=255)
                                             gAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAACTATc                                                                                                          <  2:626785/151‑1 (MQ=255)
                                              aaTTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGGGGGTTTGTAATGGG‑CGGGACCCAAATTTTTGCCGCCGCGGGCTAACAATc                                                                                                          >  1:648842/1‑151 (MQ=255)
                                              aaTTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAACTATcc                                                                                                         <  1:385201/151‑1 (MQ=255)
                                                tttGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGGTTTGTAATGGA‑CGGCACCAAATTTTTCGCCGCAGC‑GGGTAAATATCCgg                                                                                                       >  2:72700/1‑151 (MQ=255)
                                                 ttGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTTGCGGCGGGTCTTTACTGGG‑CGGCACCAAAATTATCGCCGCCGC‑GGCTAACTATTCgg                                                                                                       >  2:963873/1‑150 (MQ=255)
                                                     cGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGG‑CGGCACCCAATTTATCGCCCCAGC‑GGCTAACTATCCCGaaaa                                                                                                   >  2:870555/1‑150 (MQ=255)
                                                       gTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGGTCTGTACTGGG‑CGGCACCCAAATTTTTGCCGCAGC‑GGGTAACTATCCGGAAAata                                                                                                 >  2:414933/1‑150 (MQ=255)
                                                        tATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGGTCTTTCATGGA‑CCGCCCCAAAATTTTCCCCCCAGC‑GGCTAAACATCCCGaaaa                                                                                                   >  1:323304/1‑147 (MQ=255)
                                                            gAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGGTCTGTACTTGA‑CGGCACCAAAATTATTGCCCCAGC‑GGCTAACTATCCGGAAAAAATCGAT‑c                                                                                          >  2:500437/1‑151 (MQ=255)
                                                              gCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAACTATCCGGAAAATATCGAT‑CCg                                                                                        <  2:201394/151‑1 (MQ=255)
                                                                cAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGGTTTGTACTGGG‑CGGGACCAAAATTTTCGCCGCCGC‑GGCTAACTATTCGGAAAAAATCGAttcctg                                                                                       >  2:1074295/1‑146 (MQ=255)
                                                                 aaaCCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGGTCTGTACTGGA‑CGGCACCAAATTTTTTGCCCCAGC‑GGCTAACTATCCCGGAAAAATTGAT‑CCGTgg                                                                                     >  2:822239/1‑151 (MQ=255)
                                                                     ccGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGGTTTGTACTGGG‑CGGCACCAAAATTATTGCCGCAGC‑GGCTAACTATCCCGAAAAAATCGAT‑CCgggg                                                                                     >  2:379630/1‑147 (MQ=255)
                                                                       ggCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAACTATCCGGAAAATATCGAT‑CCGTGGCACAtt                                                                               <  1:480804/151‑1 (MQ=255)
                                                                        gCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAACTATCCGGAAAATATCGAT‑CCGTGGCACATTc                                                                              <  2:720875/151‑1 (MQ=255)
                                                                            aTGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAACTATCCGGAAAATATCGAT‑CCGTGGCACATTCTGCa                                                                          <  1:32522/151‑1 (MQ=255)
                                                                               aaaCGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAACTATCCGGAAAATATCGAT‑CCGTGGCACATTCTGCAAAc                                                                       <  2:206646/151‑1 (MQ=255)
                                                                                aaCGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAACTATCCGGAAAATATCGAT‑CCGTGGCACATTCTGCAAACg                                                                      <  1:900194/151‑1 (MQ=255)
                                                                                   gCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAACTATCCGGAAAATATCGAT‑CCGTGGCACATTCTGCAAACGGGc                                                                   <  2:605055/151‑1 (MQ=255)
                                                                                      gatgaACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGA‑CGGCACCAAATTTATCGCCACAGC‑GGCTAACTATCCGGAAAATATCGAT‑CCGTGGCACATTCTGCAAACGGGCGGt                                                                <  1:767865/151‑1 (MQ=255)
                                                                                         gaACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAACTATCCGGAAAATATCGAT‑CCGTGGCACATTCTGCAAACGGGCGGTaa                                                              <  2:1147984/150‑1 (MQ=255)
                                                                                          aaCGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGGGGCGGGTCTGTACCGGG‑CGGCACCAAATTTTTCGCCGCAGC‑GGCTAAATATCCGGGAAATATCCAT‑CCGTGGCCCACTCCGCAAAAGGGCGGTaaa                                                             >  1:55689/1‑150 (MQ=255)
                                                                                          aaCGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGGTCTGTACTTGA‑CGGCACCCAAATTTTCGCCGCAGC‑GGCTAACTATTCGGAAAATATCGAT‑TCCTGGGACATTCTGCAAAACGGCGGGaa                                                              >  2:379191/1‑149 (MQ=255)
                                                                                          aaCGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGGTCTGTACTGGA‑CGGGCCCCAATTTATTGCCGCAGC‑GGCTAAATATCCCGGAAATATCCAT‑CCGTGGCACATTCTGCAAACGGGCGGGaaa                                                             >  2:571831/1‑150 (MQ=255)
                                                                                           aCGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAACTATCCGGAAAATATCGAT‑CCGTGGCACATTCTGCAAACGGGCGGTAAaga                                                           <  2:878347/151‑1 (MQ=255)
                                                                                            cGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGGTTTGTACTGGG‑CGGCACCAAAATTATCGCCGCCGC‑GGCTAACTATCCGGGAAATATCGAT‑CCGTGGCACATTCTGCAAACGGGCGGGaaa                                                             >  2:845208/1‑148 (MQ=255)
                                                                                             gCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGG‑CGGCCCCAAATTTATTGCCGCAGC‑GGCTAACTATTCGGGAAATATTGAA‑CCCTGGCACATTCTGCAAACGGGCGGGaaa                                                             >  2:113680/1‑147 (MQ=255)
                                                                                             gCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAACTATCCGGAAAATATCGAT‑CCGTGGCACATTCTGCAAACGGGCGGTAAAgaga                                                         <  2:809541/151‑1 (MQ=255)
                                                                                                  gAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGGTCTGTACTGGACCGGCCCCAATTTTTTGCCCCCAGG‑GCGCAACTTCCCGGAAAA‑ATCGATACCGGGGCCCATTCTCCAAACggcggggaagaaaataa                                                      >  1:1102261/1‑135 (MQ=255)
                                                                                                     aCTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAACTATCCGGAAAATATCGAT‑CCGTGGCACATTCTGCAAACGGGCGGTAAAGAGATTAAAAgc                                                 <  2:428049/151‑1 (MQ=255)
                                                                                                     aCTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGGTCTTTACTGGA‑CGGCCCCAAATTTTTCGCCGCAGC‑GGCTAACTATCCGGAAAAAATCGAT‑CCGTGGCACATTTTTCAAACGGGGGGGAAAGAGATTaaaa                                                   >  2:423401/1‑149 (MQ=255)
                                                                                                         gTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGGTTTGTACTGGG‑CGGCACCAAATTTTTCGCCGCCGC‑GGCTAACTATCCGGGAAAAATCGAT‑CCGTGGCACATTCTGCAAAAGGGCGGGAAAAAGATTAAAAACGcc                                              >  2:504943/1‑150 (MQ=255)
                                                                                                          ttCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGGTCTGTACTGGG‑CGGGCCCAAAATTATTGCCGCCGC‑GGGTAACTATCCGGGAAAAATAGAT‑CCGTGGCACATTTTGCAAAAGGGCGGGAAAGAGAATaaaa                                                   >  1:63117/1‑144 (MQ=255)
                                                                                                            cgcgAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAACTATCCGGAAAATATCGAT‑CCGTGGCACATTCTGCAAACGGGCGGTAAAGAGATTAAAAGCGCCATcc                                          <  2:648874/151‑1 (MQ=255)
                                                                                                            cgcgAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAACTATCCGGAAAATATCGAT‑CCGTGGCACATTCTGCAAACGGGCGGTAAAGAGATTAAAAGCGCCATc                                           <  1:1020617/150‑1 (MQ=255)
                                                                                                             gcgAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAACTATCCGGAAAATATCGAT‑CCGTGGCACATTCTGCAAACGGGCGGTAAAGAGATTAAAAGCGCCATccc                                         <  2:770770/151‑1 (MQ=255)
                                                                                                              cgAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAACTATCCGGAAAATATCGAT‑CCGTGGCACATTCTGCAAACGGGCGGTAAAGAGATTAAAAGCGCCATCCCg                                        <  2:102511/151‑1 (MQ=255)
                                                                                                              cgAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAACTATCCGGAAAATATCGAT‑CCGTGGCACATTCTGCAAACGGGCGGTAAAGAGATTAAAAGCGCCATCCCg                                        <  2:120655/151‑1 (MQ=255)
                                                                                                                aaCAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAACTATCCGGAAAATATCGAT‑CCGTGGCACATTCTGCAAACGGGCGGTAAAGAGATTAAAAGCGCCATCCCGAt                                      <  1:116346/151‑1 (MQ=255)
                                                                                                                aaCAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAACTATCCGGAAAATATCGAT‑CCGTGGCACATTCTGCAAACGGGCGGTAAAGAGATTAAAAGCGCCATCCCGAt                                      <  2:1049005/151‑1 (MQ=255)
                                                                                                                 aCAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGGTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAAATATCCGGGAAATATCGAT‑CCGTGGCACATTCTGCAAAAGGGCGGGAAAGAGATTAAAAACGCCATCCCGa                                       >  2:759320/1‑149 (MQ=255)
                                                                                                                  cAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGGTCTGTACTGGA‑CGGGACCCAATTTATCGCCGCAGC‑GGGTAAATATCCGGGAAAAATCGAT‑CCCTGGCACCATCTGCAAACGGGGGGGAAAGAGATTAAAAACGCCATCCCGATgg                                    >  2:630441/1‑151 (MQ=255)
                                                                                                                   aGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAACTATCCGGAAAATATCGAT‑CCGTGGCACATTCTGCAAACGGGCGGTAAAGAGATTAAAAGCGCCATCCCGATggg                                   <  1:1113642/151‑1 (MQ=255)
                                                                                                                   aGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGGTTTGTACTGGA‑CGGGACCAAATTTTTCGCCGCAGC‑GGGCAAATATCCGGGAAATATCGAT‑CCGTGGCACATTCTGCAAACGGGGGGGAAAGAGGTTAAAAACGCCATACCGATggg                                   >  1:299177/1‑151 (MQ=255)
                                                                                                                   aGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGGTCTGTACTGGG‑CGGGACCAAATTTATCGCCGCAGC‑GGGTAAATATCCGGGAAATATCGAT‑CCGTGGCACATTCTGCAAAAGGGCGGGAAAGAGATTAAAAACGCCATCCCGATggg                                   >  1:90612/1‑151 (MQ=255)
                                                                                                                      aaGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAACTATCCGGAAAATATCGAT‑CCGTGGCACATTCTGCAAACGGGCGGTAAAGAGATTAAAAGCGCCATCCCGATGGGCtg                                <  2:1159457/151‑1 (MQ=255)
                                                                                                                           aCTTTCCTGCTGGCGGTTGGCGGCGGGTCTGTACTGGA‑CGGCACCCAATTTATCGCCGCAGC‑GGCTAACTATCCGGAAAATATCGAT‑CCGTGGCCCATTCTGCAAACGGGCGGGAAAGAGATTAAAAACGCCATCCCGATGGGCtgtgtg                            >  2:574182/1‑150 (MQ=255)
                                                                                                                             tttCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGA‑CGGCACCCAATTTTTCGCCGCCGC‑GGCTAACTATCCGGAAAATATCGAT‑CCGTGGCACATTCTGCAAAAGGGCGGGAAAGAGATTAAAAACGCCCTCCCGAtgggctgtgggctg                         >  1:198602/1‑151 (MQ=255)
                                                                                                                               tCCTGCTGGCGGTTGGCGGCGGTTCTTTACTGGA‑CGGGCCCAAAATTATTGCCGGCGG‑GGCTAACTATCCGGAAAAAATCCAT‑CCCTGGCACATTCTGCAAACGGGGGGGAAAAAGATTAAAAGGGCCCTCCCCATGGGCTGTGTGCTGa                        >  2:202099/1‑150 (MQ=255)
                                                                                                                                ccTGCTGGCGGTTGGCGGCGGTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAACTATCCGGAAAATATCGAT‑CCGTGGCACATTCTGCAAACGGGCGGTAAAGAGATTAAAAGCGCCATCCCGATGGGCTGTGTGCTGAc                       <  2:1002531/150‑1 (MQ=255)
                                                                                                                                 ctgctgGCGGTTGGCGGCGGGTCTTTACTGGA‑CGGCCCCAAATTTATCGCCCCCGC‑GGGTAACTATCCGGAAAATATTGAA‑CCGTGGCCCATTTTGCAAACGGGCGGGAAAAAGAATAAAAACGCCCTCCCCATGGGCTGTGTGCTg                         >  2:829181/1‑147 (MQ=255)
                                                                                                                                   gctgGCGGTTGGCGGCGGTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAACTATCCGGAAAATATCGAT‑CCGTGGCACATTCTGCAAACGGGCGGTAAAGAGATTAAAAGCGCCATCCCGATGGGCTGTGTGCTGACGCt                    <  1:499497/150‑1 (MQ=255)
                                                                                                                                      ggCGGTTGGCGGCGGTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAACTATCCGGAAAATATCGAT‑CCGTGGCACATTCTGCAAACGGGCGGTAAAGAGATTAAAAGCGCCATCCCGATGGGCTGTGTGCTGACGCTGcc                 >  2:6771/1‑150 (MQ=255)
                                                                                                                                            tGGCGGCGGTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAACTATCCGGAAAATATCGAT‑CCGTGGCACATTCTGCAAACGGGCGGTAAAGAGATTAAAAGCGCCATCCCGATGGGCTGTGTGCTGACGCTGCCAGCAAcc          >  2:773082/1‑151 (MQ=255)
                                                                                                                                               cggcggTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAACTATCCGGAAAATATCGAT‑CCGTGGCACATTCTGCAAACGGGCGGTAAAGAGATTAAAAGCGCCATCCCGATGGGCTGTGTGCTGACGCTGCCAGCAACCGGt       <  2:609206/151‑1 (MQ=255)
                                                                                                                                                ggcggTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAACTATCCGGAAAATATCGAT‑CCGTGGCACATTCTGCAAACGGGCGGTAAAGAGATTAAAAGCGCCATCCCGATGGGCTGTGTGCTGACGCTGCCAGCAACCGGtt      >  1:1006348/1‑151 (MQ=255)
                                                                                                                                                  cggTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAACTATCCGGAAAATATCGAT‑CCGTGGCACATTCTGCAAACGGGCGGTAAAGAGATTAAAAGCGCCATCCCGATGGGCTGTGTGCTGACGCTGCCAGCAACCGGTTCa    <  2:1191075/151‑1 (MQ=255)
                                                                                                                                                    gTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAACTATCCGGAAAATATCGAT‑CCGTGGCACATTCTGCAAACGGGCGGTAAAGAGATTAAAAGCGCCATCCCGATGGGCTGTGTGCTGACGCTGCCAGCAACCGGTTCAGa  <  1:969161/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                        
GTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGA‑CGGCACCAAATTTATCGCCGCAGC‑GGCTAACTATCCGGAAAATATCGAT‑CCGTGGCACATTCTGCAAACGGGCGGTAAAGAGATTAAAAGCGCCATCCCGATGGGCTGTGTGCTGACGCTGCCAGCAACCGGTTCAGA  >  NC_000913/3155490‑3155788

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