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	42.2%	252.1 / 272.0	214	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 (51/64); new base T (81/3); total (146/68)
Fisher's exact test for biased strand distribution p-value = 1.63e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-01

Only 100 of 224 total aligned reads displayed.

TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACCA  >  NC_000913/3705949‑3706147

tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCgg                                                                                                     <  1:1726720/100‑1 (MQ=255)

tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCgg                                                                                                     >  2:1836634/1‑100 (MQ=255)

tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCgg                                                                                                     <  1:820293/100‑1 (MQ=255)

 gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGa                                                                                                    <  1:2240555/100‑1 (MQ=255)

  cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAg                                                                                                   <  1:1779697/100‑1 (MQ=255)

   gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                                                                                                  <  2:2044160/100‑1 (MQ=255)

   gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                                                                                                  <  1:216799/100‑1 (MQ=255)

    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                                                                                                 >  1:433220/1‑100 (MQ=255)

    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                                                                                                 >  1:1876506/1‑100 (MQ=255)

     tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcgc                                                                                                <  2:73948/100‑1 (MQ=255)

      ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                                                                                               <  2:2571594/100‑1 (MQ=255)

       tttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGC‑‑A‑GA‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCAt                                                                                              >  1:1261428/8‑100 (MQ=14)

       tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑‑‑A‑AT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCAt                                                                                              >  1:1181570/8‑100 (MQ=18)

       tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGC‑‑AGGG‑AC‑TAT‑‑‑GTAGGCCTGATAAGCATAGCGCAt                                                                                              >  1:2202814/8‑100 (MQ=18)

       tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCACTACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                                                                                              >  2:2594305/1‑100 (MQ=255)

        ttgtaggCCGGATAAGGCGTTCGCGCCGCATCCGGCATAC‑AAAAGTTTGCAAATTCAATAAATTGC‑‑A‑GA‑AT‑TCAT‑‑GTAGGCCTGATAAGCGAAGCGCAt                                                                                              >  2:2453010/7‑99 (MQ=17)

         tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                                                                                            >  1:2072212/1‑100 (MQ=255)

         tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                                                                                            >  1:2695557/1‑100 (MQ=255)

          ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                                                                                           <  2:2337828/100‑1 (MQ=255)

           taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAACAGTGCAAATTCAATATATTGC‑‑A‑AT‑CT‑TCTG‑‑GTAGGCCTGATAAGCATCGCGCATCAgg                                                                                          >  1:1280926/4‑100 (MQ=18)

           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGC‑‑A‑GG‑AT‑ACGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                          >  1:2459294/4‑100 (MQ=18)

           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGC‑‑A‑GG‑AT‑ACGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                          >  2:2018793/4‑100 (MQ=17)

            aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGC‑‑A‑GA‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                         >  1:680191/3‑100 (MQ=17)

            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGC‑‑A‑GT‑TT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGc                                                                                         >  1:2344368/3‑100 (MQ=25)

            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGC‑‑A‑CC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                         >  1:2434943/3‑100 (MQ=21)

            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑‑‑A‑AT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                         >  1:2839048/3‑100 (MQ=18)

            aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AC‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                                                                                         >  2:1243322/3‑100 (MQ=17)

             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑A‑GG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  1:2578198/2‑100 (MQ=18)

             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑A‑GG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  1:2647303/2‑100 (MQ=18)

             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGC‑‑A‑TG‑AG‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  1:1055559/2‑100 (MQ=17)

             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGC‑‑A‑GT‑TT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa                                                                                        >  2:1327292/2‑100 (MQ=21)

             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGA‑‑T‑CG‑AT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  2:1028915/2‑100 (MQ=18)

             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGC‑‑A‑GAGAT‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  1:1217599/2‑100 (MQ=17)

             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGC‑‑A‑GA‑GA‑TCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  2:2614083/2‑100 (MQ=17)

             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑A‑AT‑CT‑CCGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  1:1586249/2‑100 (MQ=18)

             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGC‑‑A‑GG‑AT‑ACGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  2:600718/2‑100 (MQ=18)

              gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                                                       >  1:310663/1‑100 (MQ=255)

              gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGC‑‑A‑GA‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                                       >  2:1200681/1‑99 (MQ=17)

              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGT‑‑G‑GG‑CT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                                       >  2:638621/1‑99 (MQ=21)

              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGC‑‑A‑GG‑AA‑ACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCag                                                                                       >  1:879250/1‑99 (MQ=18)

              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGC‑‑A‑GG‑GA‑CTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCag                                                                                       >  2:1508364/1‑99 (MQ=18)

              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAACCGCGCAAATTCAATATATTGCGGA‑GA‑AA‑AT‑‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGctgt                                                                                      >  2:2659404/1‑97 (MQ=21)

               ccggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttt                                                                                     >  2:2089084/9‑100 (MQ=18)

               ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGC‑‑A‑GT‑TT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaat                                                                                      >  1:345152/9‑98 (MQ=25)

                cgatctaggcgTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTGC‑‑AGGA‑AC‑GAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaact                                                                                     <  1:2597716/93‑4 (MQ=25)

                 ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑‑A‑GA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                    >  1:2001385/7‑100 (MQ=18)

                 ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGC‑‑A‑GA‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                    >  1:2747323/7‑100 (MQ=17)

                 ggataaggcgTTTACGCCGCATCCGGCATAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                                                                                   >  1:2136497/7‑98 (MQ=18)

                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGC‑‑A‑GA‑GA‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                                                                                    >  2:2358888/7‑100 (MQ=17)

                  gataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGC‑‑A‑CC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaattct                                                                                   >  1:1640765/6‑95 (MQ=25)

                  gataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGC‑‑A‑GG‑AG‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttg                                                                                   >  1:1861659/6‑99 (MQ=18)

                   aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                  <  2:2841885/100‑1 (MQ=255)

                    taaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGC‑‑A‑CC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattctgt                                                                                 >  2:1276733/4‑96 (MQ=21)

                     gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                                                                                <  2:2307376/100‑1 (MQ=255)

                      cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                                               >  1:1135449/1‑100 (MQ=255)

                       ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt                                                                              >  1:828002/1‑100 (MQ=255)

                        gcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgccttt                                                                             <  2:221906/100‑1 (MQ=255)

                          gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc                                                                           <  2:1221602/100‑1 (MQ=255)

                          gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc                                                                           <  2:1876493/100‑1 (MQ=255)

                           ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                                                          <  2:1596941/100‑1 (MQ=255)

                           ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                                                          <  2:1460154/100‑1 (MQ=255)

                             aaCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                                                         <  1:1407571/99‑1 (MQ=255)

                             aaCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt                                                                        >  1:1967762/1‑100 (MQ=255)

                                 ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                                       >  2:1759215/1‑97 (MQ=255)

                                 ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt                                                                    >  1:1331398/1‑100 (MQ=255)

                                  cTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTc                                                                   <  2:1127792/100‑1 (MQ=255)

                                  cTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTc                                                                   <  2:769853/100‑1 (MQ=255)

                                      tCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGt                                                               >  2:1636189/1‑100 (MQ=255)

                                         ggCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                                                            >  2:278703/1‑100 (MQ=255)

                                          gCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATgg                                                           <  2:1708506/100‑1 (MQ=255)

                                              aC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                                                       <  1:605392/100‑1 (MQ=255)

                                               c‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttt                                                      >  2:333693/1‑100 (MQ=255)

                                                   aaTCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                                                   <  2:377670/100‑1 (MQ=255)

                                                      cGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattat                                                >  2:1085040/1‑100 (MQ=255)

                                                        tGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                              <  2:1808652/100‑1 (MQ=255)

                                                            aaTTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATa                                          <  1:97915/100‑1 (MQ=255)

                                                               tCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATagag                                       <  2:271996/100‑1 (MQ=255)

                                                               tCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATagag                                       <  2:825323/100‑1 (MQ=255)

                                                                cAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga                                      <  1:1010727/100‑1 (MQ=255)

                                                                cAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga                                      <  2:2669874/100‑1 (MQ=255)

                                                                     cATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGGCTTGTATGGCTTTTAATTATTCGATAGAGACGttt                                 >  1:1566551/1‑100 (MQ=255)

                                                                       ttGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTc                               >  1:2357347/1‑100 (MQ=255)

                                                                         gC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                             >  2:140945/1‑100 (MQ=255)

                                                                             a‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg                           >  1:1935516/1‑100 (MQ=255)

                                                                               gA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt                          <  1:1322064/100‑1 (MQ=255)

                                                                                 naT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTGGCCTTTGCCATCAGTCTTGAATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTg                         >  2:1580879/2‑100 (MQ=255)

                                                                                           gTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGTCGTTTTCGAAGTGATGTTTg                  <  2:1771758/100‑1 (MQ=255)

                                                                                           gTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTg                  <  2:1916141/100‑1 (MQ=255)

                                                                                               gCCTGATAAGCGGAGCGCATCAGGCATTTTTTTTTTTTTTTTTTTTCTTGTATGGCTTTTAATTATTTGATAGAGACGTTTTCGAAGTGATGTTTGCCt               >  2:583898/1‑99 (MQ=25)

                                                                                               gCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTa              <  2:310662/100‑1 (MQ=255)

                                                                                               gCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTa              <  1:764676/100‑1 (MQ=255)

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

                                                                                                  tGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATg           >  2:2251988/1‑100 (MQ=255)

                                                                                                  tGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATg           >  1:476510/1‑100 (MQ=255)

                                                                                                  tGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATg           <  1:2189135/100‑1 (MQ=255)

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

                                                                                                      aaGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATc       <  1:278704/100‑1 (MQ=255)

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

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

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

TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACCA  >  NC_000913/3705949‑3706147

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