BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,706,017	0	C	T	40.8%	318.4 / inf	320	intergenic (‑227/+81)	dppB/dppA	dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand): ref base C (73/101); new base T (117/3); total (216/104)
Fisher's exact test for biased strand distribution p-value = 1.59e-26
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.11e-01

Only 100 of 323 total aligned reads displayed.

GCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTT  >  NC_000913/3705918‑3706116

gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATAc                                                                                                               <  2:2978446/100‑1 (MQ=255)

   ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                                                                            <  1:3172799/100‑1 (MQ=255)

   ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                                                                            <  1:1326881/100‑1 (MQ=255)

    cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                                                                           >  2:3603274/1‑100 (MQ=255)

    cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                                                                           <  2:3235191/100‑1 (MQ=255)

     agcGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGc                                                                                                          >  1:2630652/3‑100 (MQ=255)

         ggTATACAGATTGTTGACAAAGGGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑Ga‑a                                                                                                  >  2:2254631/1‑100 (MQ=255)

             taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑At‑t‑‑‑tt                                                                                          <  2:922424/100‑1 (MQ=255)

               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑g                                                                                      <  2:1918488/100‑1 (MQ=255)

               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑g                                                                                      <  1:1017674/100‑1 (MQ=255)

                            aaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAg                                                                         <  2:1430235/100‑1 (MQ=255)

                               tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCg                                                                       <  2:663491/99‑1 (MQ=255)

                                gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCAGAtcggaaga                                                                     >  2:1885358/1‑93 (MQ=18)

                                 cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                                                                    >  1:3508771/1‑100 (MQ=255)

                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                                                                  <  1:2945038/100‑1 (MQ=255)

                                    tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcgc                                                                 >  1:2189719/1‑100 (MQ=255)

                                    tccgatctggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑A‑GG‑AC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAgcgc                                                                 <  2:3005078/91‑1 (MQ=18)

                                    cttttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGC‑‑A‑GA‑GT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAgcgc                                                                 >  1:3261556/10‑100 (MQ=14)

                                     ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCa                                                                <  1:843644/100‑1 (MQ=255)

                                     gtttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAACCGCGCAAATTCAATATATTGCGGA‑GA‑AA‑A‑‑‑T‑‑‑‑GTAGGCCTGATAAGCGTAGCGCAt                                                               >  2:1128730/9‑100 (MQ=18)

                                      cttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGC‑‑A‑TG‑AG‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCAt                                                               >  2:59665/8‑100 (MQ=17)

                                       ttgtaggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                                                              >  2:3033894/7‑100 (MQ=18)

                                        tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑AAGTTCGTGCAAATTCAATAAATTGC‑‑A‑‑A‑TA‑T‑‑‑GAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                            >  2:355080/6‑100 (MQ=18)

                                        cgtaggCCGGATAAGGCGGTCACGCCGCATCCGGCCTAC‑GAAACCGCAAGAATTCAATATATTGC‑‑A‑GG‑AG‑C‑‑‑GGT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                                                             >  1:471476/6‑100 (MQ=18)

                                         ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                                                            <  2:1211273/100‑1 (MQ=255)

                                         gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑‑A‑GA‑AA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                            >  1:584944/5‑100 (MQ=18)

                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGC‑‑A‑CC‑AT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                           >  1:1533441/4‑100 (MQ=21)

                                          tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                                                           <  2:2441015/100‑1 (MQ=255)

                                          tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                                                           <  2:2699527/100‑1 (MQ=255)

                                           aggCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCCAAATTCAATATATTGC‑‑A‑GA‑AA‑T‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                          >  1:957809/3‑100 (MQ=18)

                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTCTTCCAAATTCAATATATTGC‑‑A‑GA‑AA‑T‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                          >  2:1626610/3‑100 (MQ=17)

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

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

                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAACAGTGCAAATTCAATATATTGC‑‑A‑AT‑CT‑T‑‑‑CTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGc                                                          >  1:2736501/3‑100 (MQ=18)

                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAACAGTGCAAATTCAATATATTGC‑‑A‑AT‑CT‑T‑‑‑CTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGc                                                          >  2:3266462/3‑100 (MQ=18)

                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAAAGCATGCAAATTCAATATATTGC‑‑A‑GA‑GT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                          >  1:3042469/3‑100 (MQ=17)

                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGC‑‑A‑TG‑AG‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCAtaa                                                             >  1:1441892/3‑95 (MQ=17)

                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGC‑‑A‑TG‑AG‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                          >  2:1962226/3‑100 (MQ=17)

                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGA‑‑T‑CG‑AT‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                          >  1:1304816/3‑100 (MQ=21)

                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGC‑‑A‑GG‑AA‑A‑‑‑CAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGc                                                          >  1:685846/3‑100 (MQ=18)

                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTAC‑‑AGGG‑AC‑T‑‑‑AT‑‑‑GTAGGCCTGATAAGCATAGCGCATCAGGc                                                          >  2:2729459/3‑100 (MQ=18)

                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGC‑‑A‑GG‑AG‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  1:3414957/3‑100 (MQ=21)

                                           aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                                                          >  1:945137/1‑100 (MQ=255)

                                           aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                                                          >  1:1572263/1‑100 (MQ=255)

                                            ggtcggataaggcgTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTGC‑‑A‑GG‑AA‑C‑‑‑GAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  2:3463217/11‑100 (MQ=21)

                                            ggcaggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAATTCAATATATTGC‑‑A‑GA‑GA‑T‑‑‑GGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  2:1670160/11‑100 (MQ=1)

                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑‑GA‑ATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt                                                        >  1:234876/2‑100 (MQ=18)

                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGC‑‑A‑GA‑GT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  1:2227626/2‑100 (MQ=17)

                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑A‑GG‑AC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  1:3440470/2‑100 (MQ=18)

                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGC‑‑A‑AG‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  1:2388329/2‑100 (MQ=17)

                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGC‑‑A‑AG‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  1:1902340/2‑100 (MQ=17)

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

                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑AAGTTCGTGCAAATTCAATAAATTGC‑‑A‑‑A‑TA‑T‑‑‑GAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcga                                                        >  2:860427/2‑98 (MQ=21)

                                             ggcggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑‑‑A‑AT‑C‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                        >  2:499341/10‑99 (MQ=18)

                                             gCCGGATAAGNCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑‑A‑GA‑AA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                        >  1:1383562/1‑99 (MQ=17)

                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGC‑‑A‑GA‑GT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                        >  1:2074816/1‑99 (MQ=17)

                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGC‑‑A‑GA‑GT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                        >  1:3563343/1‑99 (MQ=17)

                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑‑‑A‑AT‑TCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa                                                         >  1:1161746/1‑100 (MQ=21)

                                              ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                       >  1:2819595/1‑100 (MQ=255)

                                               cggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGC‑‑A‑GAGAT‑T‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                                                      >  1:3005037/8‑97 (MQ=17)

                                                ggataagtcgTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGC‑‑A‑GG‑AG‑C‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  1:2041409/9‑100 (MQ=17)

                                                ggataaggcgTTTACGCCTTATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑‑GA‑ATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                                                    >  2:1512796/7‑98 (MQ=21)

                                                ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑‑A‑GA‑AA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  1:3640349/7‑100 (MQ=18)

                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑A‑GG‑AC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  1:2827041/7‑100 (MQ=18)

                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑A‑GG‑AC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  1:1993758/7‑100 (MQ=18)

                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGC‑‑A‑CC‑AT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaattc                                                     >  2:2392177/7‑96 (MQ=21)

                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGC‑‑A‑AG‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacgct                                                     >  1:2494059/7‑96 (MQ=17)

                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGC‑‑A‑GG‑AG‑C‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  1:2305955/7‑100 (MQ=17)

                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGC‑‑A‑GG‑AA‑A‑‑‑CAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCAgttt                                                     >  1:3657859/7‑100 (MQ=18)

                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGC‑‑A‑GG‑GA‑C‑‑‑TAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctt                                                     >  1:613220/7‑96 (MQ=21)

                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGC‑‑A‑GG‑AG‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat                                                    >  1:486544/7‑98 (MQ=25)

                                                ggataaggcgTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTGC‑‑AGGA‑AC‑G‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaacta                                                     >  2:1022566/7‑96 (MQ=25)

                                                ggataaggcgTTCACGCCGCATCCAGCCTAC‑AAAATTGTGCAAAATCAATGGATTGC‑‑A‑CA‑GC‑T‑‑‑AGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagatt                                                     >  2:333001/7‑96 (MQ=37)

                                                 gataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGA‑‑T‑CG‑AT‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                                                    >  1:428618/6‑100 (MQ=21)

                                                  aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATTGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGTATTtttg                                                   <  2:2044357/100‑1 (MQ=255)

                                                  aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                   <  1:1240887/100‑1 (MQ=255)

                                                           tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                          <  2:3790875/100‑1 (MQ=255)

                                                            aaCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt                                         >  2:760196/1‑100 (MQ=255)

                                                              ngCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa                                       >  2:1392847/2‑100 (MQ=255)

                                                               gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg                                      <  2:572845/100‑1 (MQ=255)

                                                                ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑TTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg                                      <  1:166732/99‑1 (MQ=255)

                                                                  ttATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt                                   <  2:253234/100‑1 (MQ=255)

                                                                  ttATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt                                   <  2:2614530/100‑1 (MQ=255)

                                                                    aTCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                 <  1:3532228/100‑1 (MQ=255)

                                                                    aTCCGCCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AA‑A‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                 >  2:2204057/1‑100 (MQ=255)

                                                                       cGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                              >  1:3365034/1‑100 (MQ=255)

                                                                       cGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                              <  2:2458682/100‑1 (MQ=255)

                                                                       cGGCCTAC‑AAAATCGTGCAAATCCAGTACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                              <  2:207506/100‑1 (MQ=255)

                                                                          ccTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                           <  2:431452/100‑1 (MQ=255)

                                                                                    tCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAatta                  <  2:3676762/100‑1 (MQ=255)

                                                                                     cGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattat                 >  1:3812896/1‑100 (MQ=255)

                                                                                       tGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc               <  2:2775197/100‑1 (MQ=255)

                                                                                         cAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCTTTTGTTATCAGTCTTGTATGGCTTTTAATTATTCGa             >  2:2844337/1‑100 (MQ=255)

                                                                                         cAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGa             >  2:3729719/1‑100 (MQ=255)

                                                                                         cAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGa             <  1:2436213/100‑1 (MQ=255)

                                                                                            aTTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATag          <  1:760198/100‑1 (MQ=255)

                                                                                            aTTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATag          <  2:1990638/100‑1 (MQ=255)

                                                                                            aTTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATag          <  2:1055197/100‑1 (MQ=255)

                                                                                              tCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATagag        >  1:1822352/1‑100 (MQ=255)

                                                                                               cAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGCTAAGTGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga       <  2:2328558/100‑1 (MQ=255)

GCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTT  >  NC_000913/3705918‑3706116

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