Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,423,2970GC41.8% 149.3 / 51.1 98intergenic (+103/+126)yhdZ/rrfFputative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base G (19/38);  new base C (41/0);  total (60/38)
Fisher's exact test for biased strand distribution p-value = 4.83e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 102 total aligned reads displayed.

CCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGC  >  NC_000913/3423147‑3423447
                                                                                                                                                       |                                                                                                                                                           
ccTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCg                                                                                                                                                             <  1:532195/151‑1 (MQ=255)
   aaaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCggg                                                                                                                                                           <  1:602278/150‑1 (MQ=255)
   aaaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCgg                                                                                                                                                            <  4:50268/149‑1 (MQ=255)
        agagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTg                                                                                                                                                     <  4:436329/151‑1 (MQ=255)
          agCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCCTCCTCCCCCCCC‑CCCTCCCGGc                                                                                                                                                          >  3:86785/1‑144 (MQ=255)
           gCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑t                                                                                                                                                 <  1:769387/150‑1 (MQ=255)
                         ttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCC‑CCCCCCCGGGCGGTG‑C‑TTCGCAA‑AGTT‑Caaa                                                                                                                                 >  1:434428/1‑150 (MQ=25)
                           ttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGAGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATcc                                                                                                                              >  4:637380/1‑151 (MQ=255)
                             aTCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGc                                                                                                                            <  4:694735/151‑1 (MQ=255)
                             aTCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCCCCGGGCCGTTG‑C‑TTCGCAA‑AGTT‑CAAATCCGc                                                                                                                            >  1:191091/1‑151 (MQ=255)
                             aTCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCC‑CCCCCCCGGCCGGTG‑C‑TTCGCAA‑CGAT‑CAAATcccc                                                                                                                            >  1:1147633/1‑149 (MQ=255)
                              tCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCC‑CCCCCCGGGCCGTTG‑T‑TTTGCAA‑CATT‑CAAAAcccc                                                                                                                            >  1:280796/1‑148 (MQ=21)
                                gCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCTCCCCCCC‑CCCCCCCGGCCGCTG‑C‑TTCGCAA‑CGTT‑CAAATccccccc                                                                                                                         >  1:785432/1‑146 (MQ=21)
                                 cAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCTGCCCCCCC‑CCCCCCCGGCCCGTG‑C‑TTCGCAA‑CCTT‑TAAATCCGCTccc                                                                                                                        >  4:7278/1‑151 (MQ=255)
                                 cAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTccc                                                                                                                        <  1:1350166/151‑1 (MQ=255)
                                 cAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCC‑CCCCCCCGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTccc                                                                                                                        >  1:1284053/1‑151 (MQ=255)
                                  aGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCcg                                                                                                                       <  2:65904/151‑1 (MQ=255)
                                    gTAATCCATTAATTGAATGTTATTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCATTTCGCTTGATGTGATGTCGGGCGGTTTTTGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCcggc                                                                                                                     <  4:1283924/151‑1 (MQ=21)
                                    gTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCC‑CCCCCCCGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATcccccccc                                                                                                                        >  1:68290/1‑148 (MQ=21)
                                    gTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCCCCCCCCC‑CCCCCCCGGCCGGTG‑C‑TTCGCCA‑CCTT‑CAAAACCGCTccccg                                                                                                                      >  4:76960/1‑148 (MQ=21)
                                       aTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCC‑CCCCCCCGGCCGTTG‑C‑TTCTCAA‑ACTT‑TAAATCCCCCCCccgcgg                                                                                                                   >  1:699548/1‑150 (MQ=21)
                                         ccATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCC‑CCCCCCCGGCCGGTG‑G‑TTTGCAA‑AATT‑CAAAACCGCCCCCCGCGGGtt                                                                                                                >  4:1403866/1‑151 (MQ=18)
                                          cATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCC‑CCCCCCCGGCCGTTG‑C‑TTCCCCA‑CGTT‑TAAAACCGCcccccccgggtt                                                                                                                >  4:1157167/1‑142 (MQ=255)
                                           aTTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCA‑CCCTCCCGGCCGTTG‑C‑TTCCCAA‑CGTT‑CAAATCCGCCCCCCGCGGGttt                                                                                                               >  1:1187447/1‑150 (MQ=18)
                                                 tGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCC‑CCCCCCCGGCCGGTT‑C‑TTCGCAA‑CATT‑TAAAACCCCCCCCCGCGGGTTTGttc                                                                                                           >  1:1141929/1‑146 (MQ=18)
                                                 tGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCGCCC‑CCCCCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCCGCGGGTTTGTCCTAc                                                                                                        >  4:505071/1‑151 (MQ=255)
                                                 tGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCC‑CCCCCCCGGGCGGTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCCCCCCGCGGGTTTGttc                                                                                                           >  1:177367/1‑146 (MQ=18)
                                                 tGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCTCCCGCCC‑CCCCCCCGGCCGGTG‑C‑TTTGCCA‑CCTT‑TCAATCCGCCCCCCGCGGGGTTTTCCt                                                                                                          >  4:316133/1‑149 (MQ=255)
                                                   aaTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTc                                                                                                      <  2:188401/151‑1 (MQ=255)
                                                   aaTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTc                                                                                                      <  1:365263/151‑1 (MQ=255)
                                                    aTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCg                                                                                                     <  4:1012772/151‑1 (MQ=255)
                                                        tAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCC‑CCCCCCCGGCCGGTT‑C‑TTCTCAA‑AGTT‑TAAATCCCCCCCCCGCGGGTTTGttcctctccgg                                                                                                   >  1:662599/1‑139 (MQ=17)
                                                         aGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGgaga                                                                                                <  1:1018497/151‑1 (MQ=255)
                                                         aGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGgag                                                                                                 <  1:1426769/150‑1 (MQ=255)
                                                         aGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCC‑CCCCCCCCGCCGGTG‑C‑TTCGCAA‑CGTTCCAAATCCCCCCCCGCGGGGTTTTTcccctccgggg                                                                                                  >  4:361096/1‑141 (MQ=255)
                                                           ttCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCC‑CCCCCCCGGGCGTTG‑C‑TTTGCAA‑CGCT‑CTAATACGCCCCCCGCGGGTTTTTCCTTCTCGGGggagg                                                                                              >  4:701619/1‑150 (MQ=17)
                                                            tCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCC‑CCCCCCCGGCCGGTG‑C‑TTCCCAA‑AGTT‑CAAAACCCCCCCCGGCGGGTTTTTTCTTCTCCGGagaggg                                                                                             >  1:1076298/1‑149 (MQ=17)
                                                            tCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTAAGGCGGGCGTCCTGCCCGCCG‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAgt                                                                                              <  4:672468/150‑1 (MQ=255)
                                                              gAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCCTCCCCCCC‑CCCCCCCGGCCGTTG‑C‑TTTTCCA‑CCCT‑TTAATCCGCCCCCCGCGGGTTTTTTCTtctcccgggaggg                                                                                             >  4:1032012/1‑140 (MQ=255)
                                                              gAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCC‑CCCCCCCGGCCGGTG‑C‑TTTGCAA‑CGTT‑CCAATCCGCTCCCCGCGGGTTTGTCCTACTCCGGggaggggt                                                                                           >  4:874658/1‑147 (MQ=17)
                                                              gAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCGCCA‑CCCCCCCGGCCGTTG‑C‑TTTGCAA‑CGTT‑TAAATCCGCCCCCCGCGGGGTTTTTCTACTCCGGAGAgtg                                                                                             >  1:644808/1‑149 (MQ=255)
                                                               aaaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCC‑CCCCCCCGGCCGTTG‑C‑TTCCCAA‑CGTT‑CAAAACCCCCCCCGGCGGGGTTTTCCTACTCCGGAGAgggt                                                                                            >  1:714536/1‑149 (MQ=17)
                                                               aaaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCTCCCGCCC‑CCCCCCCGGCCCGTG‑C‑TTCGCAA‑AGTT‑TAAATCCCCCCCCCGCGGGTTTTTCCTTCTCCGGAGAGTGtt                                                                                           >  1:429743/1‑150 (MQ=255)
                                                                aaaGCAAAAAGGCCATCTTTTCGGATGCCCTTTCGTTTGATTTAATGTCTGGCAGTTTATGGCGGGCGTCTTCCCCGCCG‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGtt                                                                                           <  4:461391/149‑1 (MQ=17)
                                                                 aaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCAc                                                                                        <  1:218710/151‑1 (MQ=255)
                                                                  aGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCAcc                                                                                       <  2:27386/151‑1 (MQ=255)
                                                                  aGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCCCCGGGGCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCCGCGGGTTTTTTCTACTCCGGGGAGGGGTCAcc                                                                                       >  1:336241/1‑151 (MQ=14)
                                                                    cAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCTGCCCCCCC‑CCCCCCCGGCCCTTG‑C‑TTCGCAA‑AGTT‑CAAATCCGCTCCCCGCGGGTTTTTCCTACTCCGGggaggggttacc                                                                                       >  2:7406/1‑141 (MQ=255)
                                                                     aaaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGac                                                                                    <  1:1070815/151‑1 (MQ=255)
                                                                       aaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTTCCCGCCC‑CCCCCCCGGCCGGTG‑C‑TTCCCAA‑CCTT‑CAAAAACCCTCCCCGCCGGTTTGTTCTCCTCGGGGGGGGGTTTAccca                                                                                     >  4:811842/1‑146 (MQ=12)
                                                                        aaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCTCCCCCCC‑CCCCCCCCGCCGCCT‑C‑TTCTCAA‑CGTT‑CCAATTCCCCCCCCGCGGGGTTTTCCTTCCCGGGGGAGGGTTCACCCacca                                                                                  >  4:1216328/1‑148 (MQ=12)
                                                                        aaGGCCATCCTGTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaa                                                                                 <  4:598817/151‑1 (MQ=255)
                                                                          ggCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTTCCCCCCC‑CCCCCCCGGCCCTTG‑C‑TTTGCAA‑AGTT‑CAAAACCCCCCCCCGCGGGTTTTTCCTCCTCGGGGGGGGGTTCACCCacaaa                                                                                 >  1:518546/1‑149 (MQ=17)
                                                                          ggCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCC‑CCCCCCCCGCCGGTG‑C‑TTCTCAACCGTC‑CAAAT‑CGCTCCCCGCGGCTTTGTCCTCCTCGGGGGAGGGTTTACCCaccaa                                                                                 >  4:151361/1‑149 (MQ=12)
                                                                              aTCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAg                                                                           <  1:801980/151‑1 (MQ=255)
                                                                              aTCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAg                                                                           <  4:647723/151‑1 (MQ=255)
                                                                              aTCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCC‑CCCCCCCGGCCGGTG‑C‑TTTGCAA‑CGTT‑CAAATCCCCCCCCCGCGGGTTTGTTCTACTCCGGGGAGGGGTCACCCacaaaaaa                                                                              >  1:714801/1‑148 (MQ=14)
                                                                               tCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCC‑CCCCCCCGGCCGGTG‑C‑TTCGCAA‑CGTT‑CAAAATCCCCCCCCGCGGGTTTTTCCTACTCCGGGGAGGGGTCACCCccaaacaa                                                                              >  1:521421/1‑147 (MQ=255)
                                                                                   ttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAG‑ATaa                                                                      <  2:33099/150‑1 (MQ=255)
                                                                                   ttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAG‑ATaa                                                                      <  2:182137/150‑1 (MQ=255)
                                                                                   ttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCGCCACCCCCCCGGCCCGTTG‑‑‑TTCGCAA‑CGTT‑TAAATCCCCTCCCCGCCGGTTTGTCCTACTCCGGGGAGTGTTTACCCGCAAACAAAAA‑ATaa                                                                      >  1:519526/1‑150 (MQ=255)
                                                                                    tCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAG‑ATaaa                                                                     <  4:177358/150‑1 (MQ=255)
                                                                                     cGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCT‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAAc                                                                   <  1:1454716/151‑1 (MQ=255)
                                                                                      ggATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGAGGGCTGCCCGCCT‑CCCTCCGGGCCGCTG‑G‑TTTGCCA‑CGGT‑TAAATCTGCCCCCGGCGGGGTTGTCCCTCTCGGGGGGGGGGTTACCGCCAAACAAAAC‑ATaaaa                                                                    >  4:174219/1‑149 (MQ=255)
                                                                                          ggCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACaaaa                                                               <  4:679132/150‑1 (MQ=255)
                                                                                             cTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCC‑CCCTCCCGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCCCCCCGCGGATTTTTCCTACTCCGGAGAGGGTTCACCGACAAACAACAG‑ATAAAACaaaa                                                               >  1:1417213/1‑147 (MQ=14)
                                                                                              tttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTCCCCCCCA‑CCCTCCCGGCCGTTG‑C‑CTTGCAA‑CGTT‑TCAAACCGCTCCCCGCGGCTTTTTCCTTCTCGGGGGAGTGTTTCCCCCCAAACAACAG‑ATAAAACAAAAAGccc                                                          >  1:1357974/1‑151 (MQ=14)
                                                                                               ttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCC‑CCCCCCCGGCCGTTG‑C‑TTTGCAA‑AGTT‑TCAAACCGCCCCCGGCGGGTTTTTTCTACTCCGGAGAGGGGTCACCCACCAACAACAa‑aaaaaaaaaaaaaccc                                                          >  1:404214/1‑145 (MQ=2)
                                                                                               ttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCA‑CCCCCCGGGCCGGTT‑C‑TTCTCAA‑AGTT‑TCAAACCCCCCCCCGCGGGTTTTTTCTACTCCGGAGAGGGGTTTCCCACAAAAAAAAG‑AAAAAAAAAAAAGccc                                                          >  2:7882/1‑150 (MQ=9)
                                                                                                tCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCA‑CCCTCCCGGCCGCTG‑C‑TTCGCCA‑CGTT‑CAAATCCCCCCCCCGCGGGTTTGTCCTACTCCGGAGAGGGGTCACCGACAAACAACAA‑ATAAAACAAAAAGccc                                                          >  1:454998/1‑149 (MQ=11)
                                                                                                 cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAg                                                        <  1:582636/150‑1 (MQ=255)
                                                                                                 cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAg                                                        <  1:1348672/150‑1 (MQ=255)
                                                                                                 cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAGt                                                       <  4:131823/151‑1 (MQ=40)
                                                                                                   cTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCTTCCCCCA‑CCCTCCGGGCCGGTGTC‑TCCGC‑A‑GGTT‑AAAAAACGCCCCCCGCGGCGTTTTGCTGACCGGGCGAGGGTTTTCCGCGCAACAAAAG‑ATAAAAAAAAAGGacc                                                          >  2:11290/1‑146 (MQ=2)
                                                                                                     tgatttgatGTCTGGCAGTTTATGGCGGGCGGCCTTCCCCCCA‑CCCCCCCGGCCGGTG‑CTTTCCCAA‑CGTC‑CAAA‑CCGCCCCCCGCGGGGTTTTTCTCCTCCGGGGGGGGGTTACCCACaaaaaaaaa‑aaaaaaaaaaaagaccccgcct                                                    >  4:803015/1‑139 (MQ=255)
                                                                                                       atttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCCGGCCGGTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCCGCGGATTTTTCCTACTCGGGAGAGTGTTCACCCACAAACAACAA‑ATAAAAAAAAAAGCCCAGTCTTcc                                                  >  1:1114742/1‑150 (MQ=14)
                                                                                                           gatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAGTCTTCCGACt                                              <  1:476457/150‑1 (MQ=38)
                                                                                                            atGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAGTCTTCCGACTGa                                            <  2:185603/151‑1 (MQ=38)
                                                                                                              gccTGGCAGTTTATGGCGGGCGTCCTGCCCGCCC‑CCCCCCCGGCCCGTT‑C‑TTCGCAA‑CCTT‑CAAAACCCCCCCCCGCGGGTTTTTTCCACTCACGAGAGCGGTCACCCACAAAAAAAAGAAAAAAAAAAAAGGcccgacttt                                                    >  1:689225/3‑136 (MQ=0)
                                                                                                              gTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCA‑CCCCCCCGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAAACCCCCCCCCGCGGGTTTGTTCTACTCCGGAGAGGGTTCACCGACAAAAAACAG‑ATAAAAAAAAAAGGCCAGTCTTCCGACTGAGc                                          >  1:811286/1‑151 (MQ=11)
                                                                                                                  ggCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtt                                       <  4:108828/150‑1 (MQ=35)
                                                                                                                      gTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtt                                  <  1:1031654/151‑1 (MQ=38)
                                                                                                                        ttATGGCGGGCGTCCTGCCCCCCA‑CCCCCCCGGCCGCTG‑C‑TTCGCAA‑AGTT‑CAAATACGCCCCCCGCGGGTTTTTTCCACTCCGGGGAGGGGTTACCGACAAAAAACAG‑AAAAAAAAAAAAGCCCAGGCTTTCGACCGAGGCTtttttttt                                >  4:415903/1‑151 (MQ=2)
                                                                                                                            ggcgggggTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTCCTCGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAttt                            <  4:1291843/151‑1 (MQ=255)
                                                                                                                             gcgggcgTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTg                           >  4:640812/1‑151 (MQ=38)
                                                                                                                             gcgggcgTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTg                           >  4:1336748/1‑151 (MQ=38)
                                                                                                                               gggcgTCCTGCCCGTCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGa                          <  4:508058/150‑1 (MQ=38)
                                                                                                                               gggcgTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGAt                         >  4:199909/1‑151 (MQ=38)
                                                                                                                                   gTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCt                     >  4:219216/1‑151 (MQ=40)
                                                                                                                                   gTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCt                     <  1:1148090/151‑1 (MQ=40)
                                                                                                                                     ccTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTgg                   <  4:975351/151‑1 (MQ=40)
                                                                                                                                       tGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc                  >  2:87163/1‑150 (MQ=40)
                                                                                                                                        gCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc                  >  4:1018139/1‑149 (MQ=40)
                                                                                                                                        gCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCa                 <  1:1182003/150‑1 (MQ=40)
                                                                                                                                         cccGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGt               >  4:975474/1‑151 (MQ=40)
                                                                                                                                           cGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAGAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTc             >  4:1245419/1‑151 (MQ=255)
                                                                                                                                           cGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTc             >  4:1043051/1‑151 (MQ=40)
                                                                                                                                              cA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCt          >  1:520311/1‑151 (MQ=40)
                                                                                                                                              cA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCt          <  4:1360689/151‑1 (MQ=40)
                                                                                                                                                       gggCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGc  >  4:1180141/1‑151 (MQ=40)
                                                                                                                                                       |                                                                                                                                                           
CCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCA‑CCCTCCGGGCCGTTG‑C‑TTCGCAA‑CGTT‑CAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGC  >  NC_000913/3423147‑3423447

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