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	40.2%	320.9 / 312.5	263	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 (72/77); new base T (98/2); total (184/79)
Fisher's exact test for biased strand distribution p-value = 1.18e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

Only 100 of 272 total aligned reads displayed.

TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACC  >  NC_000913/3705949‑3706146

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

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

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

      ttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCa                                                                                              >  2:2120934/9‑100 (MQ=18)

      ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑TAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCa                                                                                              <  1:2432192/100‑1 (MQ=255)

       tttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATC‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATc                                                                                            >  1:1584114/8‑100 (MQ=17)

        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                                                                                            <  1:1944220/100‑1 (MQ=255)

         tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCa                                                                                           >  2:2587878/6‑100 (MQ=21)

         tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGCA‑GG‑AT‑‑A‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                                                                                           >  1:863251/6‑100 (MQ=18)

          ttATGCCGGATGCGGCGTTAACGCCTTATCCTGCCTAC‑AAAATCGTGCAAATTCAATACATGGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGAGAAGCGGAGCGCATCAg                                                                                          >  1:2027745/1‑100 (MQ=255)

          gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑A‑AT‑‑TCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCa                                                                                           >  2:2191371/5‑100 (MQ=21)

          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCGTAC‑GAAATCGTGCTAATTCAATATATTGTG‑GG‑CT‑‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                                          >  1:2269360/5‑100 (MQ=21)

          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                                          >  2:2786412/5‑100 (MQ=21)

           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑A‑AT‑‑C‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                         >  1:2481115/4‑100 (MQ=18)

            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGTG‑GG‑CT‑‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                        >  2:3200515/3‑100 (MQ=18)

            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AG‑‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                       >  2:1129671/3‑100 (MQ=21)

            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AG‑‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                       >  2:1230055/3‑100 (MQ=21)

            aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                                                                                        <  2:265218/100‑1 (MQ=255)

            aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                                                                                        >  2:409009/1‑100 (MQ=255)

            aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑AA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                                                                                        >  2:2437221/1‑100 (MQ=255)

             ggtcggataaggcgTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTGCAGGA‑AC‑‑G‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                       >  2:2732408/11‑100 (MQ=21)

             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGCA‑TG‑GG‑‑C‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                       >  2:1145500/2‑100 (MQ=17)

             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AG‑‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt                                                                                      >  2:77638/2‑100 (MQ=21)

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

              gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                                      >  2:3015468/1‑99 (MQ=17)

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

              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCA‑TG‑AG‑‑C‑‑‑CAC‑‑GTAGGCCTGATAATCGTAGCGAATAAgaaaa                                                                                      >  2:1368949/1‑96 (MQ=17)

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

               ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCA‑TG‑AG‑‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                                                     >  2:2333946/9‑98 (MQ=17)

               ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TT‑‑A‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaat                                                                                     >  2:628674/9‑98 (MQ=21)

               ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGAT‑CG‑AT‑‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                                                     >  1:1359262/9‑98 (MQ=18)

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

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

               ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGCA‑AG‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGTAGCGCAGCAgg                                                                                         >  1:3045634/9‑96 (MQ=14)

               ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCA‑GG‑AG‑‑C‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                                                     >  2:2280278/9‑98 (MQ=18)

                ctgataaggcgTTTACGCCTTATCAGGCCTAC‑AAGATCGTGCAAATTCAACATATTGC‑‑CA‑CT‑‑C‑‑‑ACCCAGTAGGCCTGATAAGCGCAGCGCATCAGGCaat                                                                                     >  2:2460863/8‑98 (MQ=25)

                 ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGCA‑‑A‑TTCAA‑‑‑CTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatt                                                                                    >  2:821523/7‑97 (MQ=21)

                 ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGCA‑‑A‑TTCAA‑‑‑CTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatt                                                                                    >  2:1937670/7‑97 (MQ=21)

                 ggataaggcgTTCACGCCGCATCCGGCGTAC‑GAAATCGTGCTAATTCAATATATTGTG‑GG‑CT‑‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                   >  2:923982/7‑100 (MQ=18)

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

                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCA‑GA‑GA‑‑T‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                                                                                   >  1:324604/7‑100 (MQ=17)

                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCAGGG‑AC‑‑T‑‑‑AT‑‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctt                                                                                   >  1:2680959/7‑96 (MQ=18)

                 ggATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                                                                                   >  2:2951293/1‑100 (MQ=255)

                  gataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCA‑GAGAT‑‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattta                                                                                  >  1:2575390/6‑99 (MQ=17)

                  gataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCAGGG‑AC‑‑T‑‑‑AT‑‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagcttt                                                                                  >  2:2693744/6‑95 (MQ=21)

                   ataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAGTtttg                                                                                 >  1:642990/5‑100 (MQ=17)

                     gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                                                                               >  2:3460464/1‑100 (MQ=255)

                     gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                                                                               <  2:2561252/100‑1 (MQ=255)

                      cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCTGGCATTtttgcgt                                                                              <  2:2359807/100‑3 (MQ=255)

                       ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGAGAATCAGGCATTtttgc                                                                                >  1:3231777/1‑97 (MQ=255)

                         cgTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCA‑GA‑GA‑‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAATtttgcgtttg                                                                           >  1:131119/1‑100 (MQ=17)

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

                          gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc                                                                          >  2:950791/1‑100 (MQ=255)

                          gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc                                                                          <  1:2225914/100‑1 (MQ=255)

                            tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                                                        >  2:1653748/1‑100 (MQ=255)

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

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

                              aCGCCTTATTCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                                      <  1:96064/100‑1 (MQ=255)

                              aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                                      <  2:3412789/100‑1 (MQ=255)

                                     aTCCGGCCTAC‑AAAATCGTGTAAATTCAATACATTGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                                               <  2:1961075/100‑1 (MQ=255)

                                     aTCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                                               <  1:1813262/100‑1 (MQ=255)

                                      tCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGt                                                              <  1:1411125/100‑1 (MQ=255)

                                       ccGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTTTTGTa                                                             >  1:1289863/1‑100 (MQ=255)

                                        cGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                                                            >  1:2201309/1‑100 (MQ=255)

                                         ggCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                                                           >  2:693601/1‑100 (MQ=255)

                                          gCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATgg                                                          >  1:2942960/1‑100 (MQ=255)

                                               c‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttt                                                     <  1:561201/100‑1 (MQ=255)

                                               c‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttt                                                     <  2:3076756/100‑1 (MQ=255)

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

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

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

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

                                                             aTTCAATACATTGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATag                                        <  1:3336470/100‑1 (MQ=255)

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

                                                              ttCAATACATTGCA‑GA‑AT‑‑T‑‑‑TTC‑‑CTAGGCCTGATAAGCGGAGCGCATAATGCATATTTGCCTTTGCCATCAGGCTTGTATGGCTTTTAATTATTCGATa                                         >  1:706362/1‑98 (MQ=37)

                                                                cAATACATTGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga                                     <  1:1059153/100‑1 (MQ=255)

                                                                 aaTACATTGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGAc                                    >  2:449050/1‑100 (MQ=255)

                                                                  aTACATTGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACg                                   <  1:495953/100‑1 (MQ=255)

                                                                    aCATTGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGtt                                 >  1:1541030/1‑100 (MQ=255)

                                                                    aCATTGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGtt                                 <  1:1683071/100‑1 (MQ=255)

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

                                                                        tGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCg                             <  1:2303627/100‑1 (MQ=255)

                                                                        tGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCg                             <  2:467690/100‑1 (MQ=255)

                                                                        tGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCg                             <  1:3423959/100‑1 (MQ=255)

                                                                         gCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                            >  1:389314/1‑100 (MQ=255)

                                                                           a‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg                          <  1:608858/100‑1 (MQ=255)

                                                                              a‑aT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCGTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTg                        <  1:2904003/100‑1 (MQ=255)

                                                                                 t‑‑t‑‑‑ttC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGAt                      >  1:218669/1‑100 (MQ=255)

                                                                                        ttC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGt                    <  2:2922792/100‑1 (MQ=255)

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

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

                                                                                               aGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGAAGTTTGcc               <  2:1826174/100‑1 (MQ=255)

                                                                                                     gATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATgg         <  1:2951505/100‑1 (MQ=255)

                                                                                                      aTAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGa        >  2:2386859/1‑100 (MQ=255)

                                                                                                       tAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTCTTAGCTAGAGACGTTGTCGAAGTGATGTTTGCCTAATgggt       >  1:2399184/1‑98 (MQ=38)

                                                                                                          gCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCaa    <  2:3169765/100‑1 (MQ=255)

                                                                                                           cGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAc   <  2:1578174/100‑1 (MQ=255)

                                                                                                            ggAGCGCATCAGGCATTTTTTCTTTTTTTTTTAGTCTTGTATGGCTTTTAATTATTCGATAGTGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAcc  >  1:1518780/1‑100 (MQ=255)

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

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

TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACC  >  NC_000913/3705949‑3706146

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