BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,705,978	0	A	C	30.1%	478.6 / 260.7	281	intergenic (‑188/+120)	dppB/dppA	dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand): ref base A (89/83); new base C (74/0); total (198/83)
Fisher's exact test for biased strand distribution p-value = 3.33e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.69e-01

Only 100 of 283 total aligned reads displayed.

GAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCA  >  NC_000913/3705879‑3706076

gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                               <  1:1120067/100‑1 (MQ=255)

 naggaTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTaa                                                                                                              >  2:506530/5‑100 (MQ=255)

     aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                                          <  1:2813743/100‑1 (MQ=255)

     aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                                          <  1:985720/100‑1 (MQ=255)

     aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                                          <  1:1311204/100‑1 (MQ=255)

         ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                      >  2:2130855/1‑100 (MQ=255)

         ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                      >  1:1202138/1‑100 (MQ=255)

          tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                     >  1:651520/1‑100 (MQ=255)

              ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                                 >  2:2765359/1‑100 (MQ=255)

              ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                                 <  2:1228278/100‑1 (MQ=255)

                aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                               <  1:2708959/100‑3 (MQ=37)

                 tgtcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATa                                                                                              <  2:2597359/97‑1 (MQ=25)

                  ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATaa                                                                                             <  1:2187116/100‑2 (MQ=25)

                   gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACa                                                                                            >  1:3175436/1‑100 (MQ=255)

                    cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaa                                                                                           >  1:1146424/1‑100 (MQ=255)

                    cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaa                                                                                           <  1:1932104/100‑1 (MQ=255)

                     gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaa                                                                                          >  1:834200/1‑100 (MQ=255)

                        gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATc                                                                                       >  2:1328782/1‑100 (MQ=255)

                        gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATc                                                                                       >  1:2452839/1‑100 (MQ=255)

                         cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                                      >  1:569584/1‑100 (MQ=255)

                            gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                   <  1:1050419/100‑1 (MQ=255)

                             cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                                  <  1:273720/100‑1 (MQ=255)

                             cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                                  >  2:821864/1‑100 (MQ=255)

                                  gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                             <  1:402780/100‑1 (MQ=255)

                                  gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                             <  1:2840855/100‑1 (MQ=255)

                                   cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                            <  2:1297053/100‑1 (MQ=255)

                                   cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                            >  2:1509856/1‑100 (MQ=255)

                                      aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATa                                                                         <  1:1722279/100‑1 (MQ=255)

                                          ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                                     >  1:2396575/1‑100 (MQ=255)

                                              cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑g                                                             <  2:315224/100‑1 (MQ=255)

                                              cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑g                                                             >  2:2963667/1‑100 (MQ=255)

                                              cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑g                                                             <  1:623907/100‑1 (MQ=255)

                                              cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑g                                                             >  2:738546/1‑100 (MQ=255)

                                                     aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTc                                                 <  2:1349152/100‑1 (MQ=255)

                                                      cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑g                                              >  2:1117587/1‑100 (MQ=255)

                                                       aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑Gt                                             >  2:3192613/1‑100 (MQ=255)

                                                       aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑Gt                                             >  1:879260/1‑100 (MQ=255)

                                                        gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTcn                                                <  2:464951/98‑2 (MQ=255)

                                                        gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTa                                            >  2:30064/1‑100 (MQ=255)

                                                         aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCATGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAg                                           <  1:3094939/100‑1 (MQ=255)

                                                           tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGc                                         <  2:613185/100‑1 (MQ=255)

                                                                 cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATa                                   <  2:1142706/100‑1 (MQ=255)

                                                                 cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATa                                   <  2:3244934/100‑1 (MQ=255)

                                                                 cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATa                                   >  1:2090669/1‑100 (MQ=255)

                                                                  aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATaa                                  <  2:1583110/100‑1 (MQ=255)

                                                                     gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCg                               <  1:1566178/100‑1 (MQ=255)

                                                                      tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCgg                              <  2:2390336/100‑1 (MQ=255)

                                                                       gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGa                             >  1:1339308/1‑100 (MQ=255)

                                                                       gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGa                             >  1:2766322/1‑100 (MQ=255)

                                                                        ngcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGa                             >  2:1417806/2‑99 (MQ=255)

                                                                          cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                          >  1:2258854/1‑100 (MQ=255)

                                                                           tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcgc                         <  1:2920300/100‑1 (MQ=255)

                                                                            ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGCTAAGCGGAGCGCa                        <  2:2019973/100‑1 (MQ=255)

                                                                             tttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCATACAAAAGCGTGCAAATTCAATACATTG‑C‑‑A‑TG‑GG‑C‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCAt                       >  2:1844463/8‑100 (MQ=17)

                                                                             tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTG‑C‑‑A‑CC‑AT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCAt                       >  2:510925/8‑100 (MQ=18)

                                                                              ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTG‑C‑‑A‑GAGAT‑T‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATc                      >  2:2345972/7‑100 (MQ=17)

                                                                               tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                     <  1:725937/100‑1 (MQ=255)

                                                                               tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTC‑TGCAAATTCAACAAATTG‑C‑‑‑‑GA‑ATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                    >  1:3003049/6‑100 (MQ=18)

                                                                               ngtaggCCGGATAAGGCGTTTACGCCGCATCTGGCCTACAAATTCTTCCAAATTCAATATATTG‑C‑‑A‑GA‑AA‑T‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                     >  2:696428/6‑100 (MQ=18)

                                                                               cgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCGTACGAAATCGTGCTAATTCAATATATTG‑T‑‑G‑GG‑CT‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCa                     >  2:2242991/6‑100 (MQ=21)

                                                                                ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                    <  1:1408333/100‑1 (MQ=255)

                                                                                gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTG‑C‑‑A‑GA‑AA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                    >  2:2899330/5‑100 (MQ=18)

                                                                                gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAACCGCGCAAATTCAATATATTG‑CGGA‑GA‑AA‑A‑‑‑T‑‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                   >  1:794244/5‑100 (MQ=18)

                                                                                 taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTG‑C‑‑A‑AT‑TC‑T‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                   >  2:1472395/4‑100 (MQ=18)

                                                                                 taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTG‑C‑‑A‑TG‑AG‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                   >  1:2283585/4‑100 (MQ=17)

                                                                                 taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTG‑C‑‑A‑GG‑AG‑C‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                   >  1:1747510/4‑100 (MQ=18)

                                                                                 taggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                   >  2:2439937/4‑100 (MQ=18)

                                                                                  aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTG‑C‑‑ATGG‑GC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                  >  2:100649/3‑100 (MQ=17)

                                                                                  aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTG‑C‑‑A‑GA‑GT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                  >  2:1550083/3‑100 (MQ=17)

                                                                                  aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTG‑C‑‑A‑GA‑GT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                  >  1:1225830/3‑100 (MQ=17)

                                                                                  aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTG‑C‑‑A‑CC‑AT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                  >  1:1661117/3‑100 (MQ=21)

                                                                                  aggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAGTTCGTGCAAATTCAATAAATTG‑C‑‑A‑‑A‑TA‑T‑‑‑GAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcg                 >  2:2281968/3‑99 (MQ=18)

                                                                                  aggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑AAACCGCGCAAATTCAATATATTG‑CGGA‑GA‑AA‑A‑‑‑T‑‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGct                 >  2:1160865/3‑99 (MQ=18)

                                                                                  aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                  >  2:19963/1‑100 (MQ=255)

                                                                                  aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                  >  1:523481/1‑100 (MQ=255)

                                                                                   ggtcggataaggcgTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTG‑C‑‑‑‑‑A‑AT‑C‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                 >  2:2488120/11‑100 (MQ=18)

                                                                                   ggacggataaggcgTTCACGCCGCATCCGTCCTACAAGTTCATGCAATTTTAATAAGTTGCC‑‑A‑GA‑TA‑T‑‑‑CCT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                  >  1:2423221/11‑100 (MQ=25)

                                                                                   ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTG‑C‑‑ATGG‑GC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                 >  1:3200750/2‑100 (MQ=17)

                                                                                   ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTG‑A‑‑T‑CG‑AT‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                 >  1:35164/2‑100 (MQ=21)

                                                                                   ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTG‑C‑‑A‑CC‑AT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                 >  1:1470640/2‑100 (MQ=21)

                                                                                   ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTG‑C‑‑A‑GG‑AG‑C‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                 >  1:774328/2‑100 (MQ=17)

                                                                                   ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTG‑C‑‑A‑GG‑AG‑C‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                 >  1:2661160/2‑100 (MQ=18)

                                                                                    gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                >  1:2120521/1‑100 (MQ=255)

                                                                                    gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTG‑C‑‑A‑GA‑AA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                >  2:177402/1‑99 (MQ=18)

                                                                                    gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTG‑C‑‑‑‑‑A‑AT‑T‑‑‑CTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGctg                >  1:1821118/1‑98 (MQ=18)

                                                                                    gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTG‑C‑‑‑‑‑A‑AT‑TCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa                 >  1:3137849/1‑100 (MQ=21)

                                                                                     ccggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTG‑C‑‑ATGG‑GC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat               >  1:2537022/9‑98 (MQ=17)

                                                                                      cggataaggcgTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTG‑C‑‑A‑GA‑GA‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAatt              >  2:2807920/8‑100 (MQ=17)

                                                                                       ggataaggcgTTTACGCCGCATCCGGCCTACAAATTC‑TGCAAATTCAACAAATTG‑C‑‑‑‑GA‑ATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc            >  2:2689001/7‑98 (MQ=21)

                                                                                       ggataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTG‑C‑‑A‑GA‑AA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt             >  2:367496/7‑100 (MQ=18)

                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTG‑C‑‑A‑TG‑AG‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt             >  1:2153548/7‑100 (MQ=17)

                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTG‑C‑‑A‑GT‑TT‑A‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct             >  1:375868/7‑96 (MQ=25)

                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAATCATGCAAATCCAAAATATTG‑C‑‑A‑GG‑AG‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat            >  2:3263145/7‑98 (MQ=25)

                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAATCATGCAAATCCAAAATATTG‑C‑‑A‑GG‑AG‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat            >  2:1733069/7‑98 (MQ=21)

                                                                                         aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg           <  2:381869/100‑1 (MQ=255)

                                                                                          taaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTG‑C‑‑AGGG‑AC‑T‑‑‑AT‑‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctttac          >  2:2160413/4‑93 (MQ=21)

                                                                                          taaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTG‑C‑‑AGGA‑AC‑G‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaactacgt          >  1:3191400/4‑93 (MQ=25)

                                                                                            cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct        >  2:1045833/1‑100 (MQ=255)

                                                                                            aggcgTTCACGCCGCATCCGGCCTAC‑AAATCATGCAAATCCAAAATATTG‑C‑‑A‑GG‑AG‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCATTtatgcgtt       >  1:1497246/2‑97 (MQ=25)

                                                                                                gTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc    >  1:2679524/1‑100 (MQ=255)

GAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCA  >  NC_000913/3705879‑3706076

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