BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,706,048	0	G	T	31.3%	356.8 / 169.0	189	intergenic (‑258/+50)	dppB/dppA	dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand): ref base G (61/62); new base T (55/1); total (126/63)
Fisher's exact test for biased strand distribution p-value = 3.83e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01

Only 100 of 194 total aligned reads displayed.

TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACCA  >  NC_000913/3705949‑3706147

tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCgg                                                                                                     <  1:1138425/100‑1 (MQ=255)

tGCGCTTTGTTTATGCCGGATGCAGCATTAACGCCTTATCCGGCATAC‑AAAATCGTGCAAATTCGAGAAATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCgg                                                                                                     <  1:1783048/100‑1 (MQ=17)

 gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGa                                                                                                    >  2:1422895/1‑100 (MQ=255)

  cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCGATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                                                                                                   <  1:659977/100‑1 (MQ=255)

   gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                                                                                                  <  1:2398827/100‑1 (MQ=255)

   gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                                                                                                  <  1:2239602/100‑1 (MQ=255)

   gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                                                                                                  <  1:252980/100‑1 (MQ=255)

    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                                                                                                 <  1:1821906/100‑1 (MQ=255)

    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                                                                                                 <  1:1422889/100‑1 (MQ=255)

    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                                                                                                 >  2:763133/1‑100 (MQ=255)

    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgag                                                                                                 >  2:972100/1‑98 (MQ=255)

        ttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑A‑AT‑TCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCAt                                                                                              >  1:1993749/7‑100 (MQ=21)

        ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TT‑A‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATc                                                                                             >  2:2350799/7‑100 (MQ=21)

          ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                                                                                           <  2:2941901/100‑1 (MQ=255)

          ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                                                                                           <  1:2478080/100‑1 (MQ=255)

          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGTTCGTGCAAATTCAATAAATTGCA‑‑A‑TA‑T‑‑‑GAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                          >  1:1945164/5‑100 (MQ=18)

          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCA‑GG‑AG‑C‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                                           >  1:1629410/5‑100 (MQ=18)

          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AG‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                          >  2:253099/5‑100 (MQ=21)

           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGAT‑CG‑AT‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                          >  2:373817/4‑100 (MQ=21)

           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCAGGG‑AC‑T‑‑‑AT‑‑‑GTAGGCCTGATAAGCATAGCGCATCAgg                                                                                          >  1:2565217/4‑100 (MQ=18)

            aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  2:2653790/3‑100 (MQ=18)

            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCA‑GAGAT‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                         >  2:2804013/3‑100 (MQ=17)

            aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                                                                                         >  2:2463611/3‑100 (MQ=18)

             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCA‑GA‑GA‑T‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  1:2921783/2‑100 (MQ=17)

             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAACAGTGCAAATTCAATATATTGCA‑AT‑CT‑T‑‑‑CTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGCa                                                                                        >  1:472855/2‑100 (MQ=18)

             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGTG‑GG‑CT‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCaagatcgg                                                                                        >  1:781994/2‑93 (MQ=21)

             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑A‑AT‑C‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  2:576834/2‑100 (MQ=18)

             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCA‑GG‑AA‑A‑‑‑CAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCa                                                                                        >  2:545733/2‑100 (MQ=18)

              gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGGC‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATCAGCAGAGCGCATCAAGcct                                                                                       >  1:283062/1‑98 (MQ=17)

              gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                                                       >  2:285078/1‑100 (MQ=255)

              gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                                       <  1:695210/100‑2 (MQ=17)

              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                                       >  2:944004/1‑99 (MQ=18)

              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAACTTCAATGTATTGCA‑GT‑TT‑A‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa                                                                                       >  1:1369201/1‑99 (MQ=18)

              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AG‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                                                                                      >  1:990152/1‑100 (MQ=21)

               ccggataaggcgTTCACGCCGCATCCGGCTTAC‑AAAACTTTGCAAATTCAATATATTGCA‑AT‑CT‑C‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                                                      >  2:2391848/9‑98 (MQ=18)

               ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                                                      >  2:1167708/9‑98 (MQ=21)

               ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGCA‑AG‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacg                                                                                      >  2:2012178/9‑98 (MQ=17)

               ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCA‑GAGAT‑T‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                                                      >  2:719438/9‑98 (MQ=17)

                cggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGCATGG‑GC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                                                                                     >  2:2476859/8‑97 (MQ=18)

                cggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAACAGTGCAAATTCAATATATTGCA‑AT‑CT‑T‑‑‑CTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGCaagc                                                                                     >  2:2575016/8‑97 (MQ=21)

                cggataaggcgTTCACGCCTTATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTGCA‑GG‑AA‑C‑‑‑GAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaact                                                                                     >  1:1586176/8‑97 (MQ=21)

                cggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCA‑GG‑AG‑C‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagtt                                                                                     >  1:2138807/8‑97 (MQ=17)

                cggataaggcgTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTGCAGGA‑AC‑G‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaact                                                                                     >  2:1874620/8‑97 (MQ=25)

                 ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                    >  2:1085290/7‑100 (MQ=18)

                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCACCAGGCAgttt                                                                                    >  2:1600501/7‑100 (MQ=18)

                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGTG‑GG‑CT‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                    >  1:1695713/7‑100 (MQ=21)

                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGCA‑AG‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacgct                                                                                    >  1:1181880/7‑96 (MQ=17)

                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCA‑GG‑AG‑C‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                    >  2:498307/7‑100 (MQ=18)

                 ggataaggcgTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTGCAGGA‑AC‑G‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaacta                                                                                    >  2:1798008/7‑96 (MQ=25)

                 ggataaggcgTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                                                                                    >  2:2082386/7‑100 (MQ=18)

                  gataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                                                                                   >  2:1668316/6‑100 (MQ=18)

                    tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                                                                                 <  2:1510431/100‑1 (MQ=255)

                    tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                                                                                 >  2:608428/1‑100 (MQ=255)

                      cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                                               >  2:2869997/1‑100 (MQ=255)

                      cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                                               >  2:1828657/1‑100 (MQ=255)

                       ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt                                                                              >  2:2099499/1‑100 (MQ=255)

                         cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                                            <  1:1543176/100‑1 (MQ=255)

                          gTTAACTCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                                            <  1:2065753/99‑1 (MQ=255)

                             aaCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt                                                                        >  2:1278732/1‑100 (MQ=255)

                                 ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt                                                                    >  2:2636226/1‑100 (MQ=255)

                                 ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt                                                                    >  1:358117/1‑100 (MQ=255)

                                  cTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTc                                                                   >  2:1687948/1‑100 (MQ=255)

                                    tATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt                                                                 >  2:2970679/1‑100 (MQ=255)

                                     aTCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCTTTTGCTATCAGt                                                                    >  1:583376/1‑96 (MQ=255)

                                       ccGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                                              >  2:14808/1‑100 (MQ=255)

                                           ccTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                                                          <  2:2187829/100‑1 (MQ=255)

                                             tAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                                                        >  2:1927199/1‑100 (MQ=255)

                                             tAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                                                        <  1:2750765/100‑1 (MQ=255)

                                             tAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                                                        >  2:830033/1‑100 (MQ=255)

                                                  aaaTCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTaa                                                    >  2:432423/1‑100 (MQ=255)

                                                  aaaTCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTaa                                                    <  1:862190/100‑1 (MQ=255)

                                                   aaTCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                                                   <  2:13296/100‑1 (MQ=255)

                                                       gTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt                                               >  1:371993/1‑100 (MQ=255)

                                                        tGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                              <  1:2869980/100‑1 (MQ=255)

                                                        tGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                              <  2:896466/100‑1 (MQ=255)

                                                          cAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGa                                            <  2:1620272/100‑1 (MQ=255)

                                                           aaaTTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt                                           <  2:2545178/100‑1 (MQ=255)

                                                             aTTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATag                                         >  2:2429042/1‑100 (MQ=255)

                                                              ttCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTTCTTTTTCCATCAGTCTTGTATGGCTTTTAATTATTCGATaga                                        >  1:2235914/1‑100 (MQ=255)

                                                                     cATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGttt                                 <  1:2688906/100‑1 (MQ=255)

                                                                      aTTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGtttt                                <  2:781397/100‑1 (MQ=255)

                                                                      aTTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGAGGAGAGCATAAGGCATCTATGCCTTTGCCATTCGTCTTGTATGGCTTTTAATTATTCGATAGAGACGtttt                                >  1:124849/1‑100 (MQ=255)

                                                                       ttGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTc                               >  1:31853/1‑100 (MQ=255)

                                                                                aT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGa                        >  1:268522/1‑100 (MQ=255)

                                                                                         c‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGttt                   >  1:1086301/1‑100 (MQ=255)

                                                                                            gTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTg                  >  2:2270342/1‑100 (MQ=255)

                                                                                              aGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGcc                <  1:945892/100‑1 (MQ=255)

                                                                                              aGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGcc                <  1:2775229/100‑1 (MQ=255)

                                                                                              aGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGcc                <  1:2463600/100‑1 (MQ=255)

                                                                                               ggCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCt               <  1:1614796/100‑1 (MQ=255)

                                                                                               ggCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCt               >  2:676321/1‑100 (MQ=255)

                                                                                               ggCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCt               >  2:2348151/1‑100 (MQ=255)

                                                                                                gCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTa              >  1:1870067/1‑100 (MQ=255)

                                                                                                  cTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAAt            <  1:2282417/100‑1 (MQ=255)

                                                                                                  cTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAAt            >  1:1110311/1‑100 (MQ=255)

                                                                                                    gATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATgg          >  2:2677589/1‑100 (MQ=255)

                                                                                                      tAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGAt        >  2:1389996/1‑100 (MQ=255)

                                                                                                        aGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCa      <  1:1131156/100‑1 (MQ=255)

                                                                                                           ggAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAcc   <  1:2590638/100‑1 (MQ=255)

                                                                                                            gAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACca  <  2:358118/100‑1 (MQ=255)

TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACCA  >  NC_000913/3705949‑3706147

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 22 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`