BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,002,448	0	T	A	35.7%	108.5 / 229.1	199	noncoding (966/1540 nt)	rrs3	16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base T (68/60); new base A (35/36); total (103/96)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.01e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 206 total aligned reads displayed.

GGTTGCGCTCGTTGCGGGACTTAACCCAACATCTCACGACACGAGCTGACGACAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTCAATTCCTTTGAGTTTCAACCTTGCGGTCGTACTCCCCAGGCGGAGTGCTTAATGCGTTAGCTGCAGCACTAAGGGGCGGAAACCCCCTAACACTTAGC  >  CP000730/2002298‑2002589

ggTTGCGCTCGTTGCGGGACTTAACCCAACATCTCACGACACGAGCTGACGACAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCt                                                                                                                                               <  2:436518/151‑1 (MQ=32)

ggTTGCGCTCGTTGCGGGACTTAACCCAACATCTCACGACACGAGCTGACGACAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCt                                                                                                                                               >  2:167980/1‑151 (MQ=31)

 gTTGCGCTCGTTGCGGGACTTAACCCAACATCTCACGACACGAGCTGACGACAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCtt                                                                                                                                              >  1:748756/1‑151 (MQ=32)

     cgcTCGTTGCGGGACTTAACCCAACATCTCCCGCCAAGAGCTGACGACAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGaa                                                                                                                                          >  2:1061816/1‑151 (MQ=12)

     cgcTCGTTGCGGGACTTAACCCAACATCTCACGACACGAGCTGACGACAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGaa                                                                                                                                          <  1:781156/151‑1 (MQ=32)

      gcTCGTTGCGGGACTTAACCCAACATCTCACGACACGAGCTGACGACAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAAt                                                                                                                                         <  2:236272/151‑1 (MQ=31)

         cGTTGCGGGACTTAACCCAACATCTCACGACACGAGCTGACGACAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTaa                                                                                                                                      >  1:718931/1‑151 (MQ=32)

         cGTTGCGGGACTTAACCCAACATCTCACGACACGAGCTGACGACAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTaa                                                                                                                                      >  1:1056236/1‑151 (MQ=32)

          gTTGCGGGACTTAACCCAACATCTCACGACACGAGCTGACGACAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTaaa                                                                                                                                     <  1:339162/151‑1 (MQ=32)

           ttGCGGGACTTAACCCAACATCTCACGACACGAGCTGACGACAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAAc                                                                                                                                    <  1:389087/151‑1 (MQ=32)

                ggACTTAACCCAACATCTCACGACACGAGCTGACGACAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTACACCACAt                                                                                                                               >  2:444887/1‑151 (MQ=14)

                 gACTTAACCCAACATCTCACGACACGAGCTGACGACAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATg                                                                                                                              >  2:1026370/1‑151 (MQ=32)

                   cTTAACCCAACATCTCACGACACGAGCTGACGACAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCt                                                                                                                            >  1:959866/1‑151 (MQ=32)

                   cTTAACCCAACATCTCACGACACGAGCTGACGACAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCt                                                                                                                            >  2:84886/1‑151 (MQ=32)

                    ttAACCCAACATCTCACGACACGAGCTGACGACAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTc                                                                                                                           <  2:676438/151‑1 (MQ=31)

                       aCCCGACATCTCACGACTCGAGCTGACGACAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGTGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCAc                                                                                                                        <  2:81026/151‑1 (MQ=12)

                         ccAACATCTCACGACACGAGCTGACGACAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCAccc                                                                                                                      >  2:506758/1‑150 (MQ=17)

                         ccAACATCTCACGACACGAGCTGACGACAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCg                                                                                                                      >  1:632651/1‑151 (MQ=32)

                         ccAACATCTCACGACACGAGCTGACGACAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCg                                                                                                                      >  1:706166/1‑151 (MQ=32)

                            aCATCTCACGACACGAGCTGACGACAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCtt                                                                                                                   <  1:1061761/151‑1 (MQ=31)

                             cATCTCACGACACGAGCTGACGACAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTtg                                                                                                                  <  1:427342/151‑1 (MQ=32)

                               tctcACGACACGAGCTGACGACAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTtgtg                                                                                                                >  1:42296/1‑151 (MQ=31)

                               tctcACGACACGAGCTGACGACAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTtgtg                                                                                                                >  1:315309/1‑151 (MQ=31)

                                ctcACGACACGAGCTGACGACAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGc                                                                                                               >  1:162139/1‑151 (MQ=32)

                                ctcACGACACGAGCTGACGACAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGc                                                                                                               >  1:993914/1‑151 (MQ=32)

                                 tcACGATACGAGCTGACTACAACCATGCACCACCTGTCACTTGGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCg                                                                                                              <  2:1015225/151‑1 (MQ=9)

                                 tcACGACACGAGCTGACGACAACCATGCACCACCTGTCACTTTGTCCCTCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCg                                                                                                              <  2:990406/151‑1 (MQ=17)

                                     gacacgaGCTGACGACAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTc                                                                                                          <  1:377852/151‑1 (MQ=32)

                                      acacgaGCTGACGACAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTcc                                                                                                         <  1:461162/151‑1 (MQ=32)

                                         cgaGCTGACGACAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCg                                                                                                      <  2:419704/151‑1 (MQ=32)

                                           aGCTGACGACAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCCCCACCGCTTGTGCGGGTCCCCGTc                                                                                                    >  2:944496/1‑151 (MQ=14)

                                                 cgacAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTCAATTcc                                                                                              >  2:201365/1‑151 (MQ=32)

                                                   acAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTCAATTCCtt                                                                                            >  1:123876/1‑151 (MQ=32)

                                                   acAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGATGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTCAATTCCtt                                                                                            >  2:373688/1‑151 (MQ=17)

                                                    cAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTCAATTCCttt                                                                                           >  1:324202/1‑151 (MQ=32)

                                                        cATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTCAATTCCTTTGAGt                                                                                       >  1:86135/1‑151 (MQ=32)

                                                         aTGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTCAATTCCTTTGAGtt                                                                                      >  2:629645/1‑151 (MQ=32)

                                                          tGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTCAATTcc                                                                                              >  2:197084/1‑142 (MQ=32)

                                                           gCACCACCTGTCACGTTGGCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTCAATTCCTTTGAGTTTc                                                                                    <  2:844066/151‑1 (MQ=14)

                                                               caccTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTCAATTCCTTTGAGTTTCAAcc                                                                                >  1:779529/1‑151 (MQ=32)

                                                                 ccTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTCAATTCCTTTGAGTTTCAACCtt                                                                              <  2:228120/151‑1 (MQ=32)

                                                                   tGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGGTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTCAATTCCTTTGAGTTTCAACCTTGc                                                                            >  1:733470/1‑151 (MQ=14)

                                                                    gTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTCAATTCCTTTGAGTTTCAACCTTGCg                                                                           >  1:562557/1‑151 (MQ=32)

                                                                    gTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTCAATTCCTTTGAGTTTCAACCTTGCg                                                                           >  1:346685/1‑151 (MQ=32)

                                                                    gTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTCAATTCCTTTGAGTTTCAACCTTGCg                                                                           >  1:786671/1‑151 (MQ=32)

                                                                    gTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTCAATTCCTTTGAGTTTCAACCTTGCg                                                                           >  2:296614/1‑151 (MQ=32)

                                                                     tCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTCAATTCCTTTGAGTTTCAACCTTGCgg                                                                          <  1:73199/151‑1 (MQ=32)

                                                                      cACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTCAATTCCTTTGAGTTTCAACCTTGCGGt                                                                         <  2:125098/151‑1 (MQ=31)

                                                                         tttGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTCAATTCCTTTGAGTTTCAACCTTGCGGTCGt                                                                      <  1:27453/151‑1 (MQ=32)

                                                                           tGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTCAATTCCTTTGAGTTTCAACCTTGCGGTCGTAc                                                                    >  1:654864/1‑151 (MQ=31)

                                                                           tGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTCAATTCCTTTGAGTTTCAACCTTGCGGTCGTAc                                                                    >  1:185547/1‑151 (MQ=32)

                                                                            gTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGGCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTCAATTCCTTTGAGTTTCAACCTTGCGGTCGTACt                                                                   >  2:751228/1‑151 (MQ=17)

                                                                              cccccGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTCAATTCCTTTGAGTTTCAACCTTGCGGTCGTACTcc                                                                 <  1:1001052/151‑1 (MQ=32)

                                                                               ccccGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTCAATTCCTTTGAGTTTCAACCTTGCGGTCGTACTccc                                                                >  2:592807/1‑151 (MQ=31)

                                                                                cccGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTCAATTCCTTTGAGTTTCAACCTTGCGGTCGTACTcccc                                                               <  2:517954/151‑1 (MQ=31)

                                                                                cccGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTAAATACCTTTGAGTTTCAACCTTGCGGTCGTACTcccc                                                               >  1:67193/1‑151 (MQ=17)

                                                                                       gggAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTCAATTCCTTTGAGTTTCAACCTTGCGGTCGTACTCCCCAGGCGGa                                                        <  2:901860/151‑1 (MQ=31)

                                                                                       gggAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTCAATTCCTTTGAGTTTCAACCTTGCGGTCGTACTCCCCAGGCGGa                                                        <  2:58758/151‑1 (MQ=32)

                                                                                        ggAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTCAATTCCTTTGAGTTTCAACCTTGCGGTCGTACTCCCCAGGCGGAg                                                       <  1:844827/151‑1 (MQ=32)

                                                                                        ggAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGAGTCCCCGACAATTCCTTTGAGTTACAGCCTTGCGGTCGTACTCCCCAGGCGGAg                                                       <  2:776226/151‑1 (MQ=9)

                                                                                         gAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTCAATTCCTTTGAGTTTCAACCTTGCGGTCGTACTCCCCAGGCGGAGt                                                      >  2:986511/1‑151 (MQ=31)

                                                                                         gAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTAAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTCAATTCCTTTGAGTTTCAACCTTGCGGTCGTACTCCCCAGGCGGAGt                                                      >  1:220146/1‑151 (MQ=14)

                                                                                           aaGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTCAATTCCTTTGAGTTTCAACCTTGCGGTCGTACTCCCCAGGCGGAGTGc                                                    >  1:486633/1‑151 (MQ=31)

                                                                                            aGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTCAATTCCTTTGAGTTTCAACCTTGCGGTCGTACTCCCCAGGCGGAGTGCt                                                   <  1:356299/151‑1 (MQ=11)

                                                                                                      ctgaaaacttccgtggatgtcaagaccaggTAAGGTTCTTCGCGTTGCATCGAATTAAACCACCTGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgtta                                         <  2:868863‑M1/137‑17 (MQ=255)

                                                                                                      ctgaaaacttccgtggatgtcaagaccaggTAAGGTTCTTCGCGTTGCATCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgtta                                         <  2:828033‑M1/137‑17 (MQ=255)

                                                                                                        gaaaacttccgtggatgtcaagaccaGGTAAGGTTCTTCGCGTTGCATCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgttagc                                       >  1:984799‑M1/13‑133 (MQ=255)

                                                                                                           aagttctgtggatgtcaagaCCAGGTAAGGTTCTTCGCGTTGCATCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgttagctcc                                    <  1:411218‑M1/142‑22 (MQ=255)

                                                                                                             cttccgtggatgtcaaGACCAGGTAAGGTTCTTCGCGTTGCATCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgttagctccgg                                  >  1:204954‑M1/8‑128 (MQ=255)

                                                                                                              ttctgtggatgtcaAGACCAGGTAAGGTTCTTCGCGTTGCATCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgtt                                          <  1:527630‑M1/136‑16 (MQ=255)

                                                                                                               tccgtggatgtcAAGACCAGGTAAGGTTCTTCGCGTTGCATCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgttagctccggaa                                >  1:863435‑M1/6‑126 (MQ=255)

                                                                                                                aggtggatgtCAAGACCAGGTAAGGTTCTTCGCGCTGCATCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgttagctccggaag                               >  2:148135‑M1/5‑125 (MQ=255)

                                                                                                                  gtggatGTCAAGACCAGGTAAGGTTCTTCGCGTTGCATCGAATTAAACCACATGCTCCGCCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgttagctccggaagcc                             <  2:459995‑M1/149‑29 (MQ=255)

                                                                                                                  gtggatGTCAAGACCAGGTAAGGTTCTTCGCGTTGCATCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgttagctccggaagcc                             <  2:654624‑M1/149‑29 (MQ=255)

                                                                                                                      atgtcaagaccaggtaaggttcttcgCGTTGCATCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgttagctccggaagccacgc                         <  2:803426‑M1/139‑33 (MQ=255)

                                                                                                                       tgtcaagaccaggtaaggttcttcGCGTTGCATCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgttagctccggaagccacgcc                        >  2:854013‑M1/12‑118 (MQ=255)

                                                                                                                       tgtcaagaccaggtaaggttcttcGCGTTGCATCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgttagctccggaagccacgcc                        >  2:910751‑M1/12‑118 (MQ=255)

                                                                                                                           aagaccaggtaaggttCTTCGCGTTGCATCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgttagctccggaagccacgcctcaa                    >  2:248664‑M1/8‑114 (MQ=255)

                                                                                                                           aagaccaggtaaggttCTTCGCGTTGCATCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgttagcgccggaagccacgcctcan                    <  2:130050‑M1/144‑38 (MQ=255)

                                                                                                                           aagaccaggtaaggttCTTCGCGTTGCATCGAATTAAACCACATGCTCCACCGCTCTCTTAGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgttagctccggaagccacgcctcaa                    <  1:766734‑M1/144‑38 (MQ=255)

                                                                                                                             gaccaggtaaggTTCTTCGCGTTGCATCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgttagctccggaagccacgcctcaagg                  >  1:834836‑M1/6‑112 (MQ=255)

                                                                                                                               ccaggtaaGGTTCTTCGCGTTGCATCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgttagctccggaagccacgcctcaagggc                <  2:927636‑M1/148‑42 (MQ=255)

                                                                                                                               ccaggtaaGGTTCTTCGCGTTGCATCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgttagctccggaagccacgcctcaagggc                <  1:583972‑M1/148‑42 (MQ=255)

                                                                                                                                caggtaAGGTTCTTCGCGTTGCATCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgttagctccggaagccacgcctcaagggca               >  2:566835‑M1/3‑109 (MQ=255)

                                                                                                                                 aggtAAGGTTCTTCGCGTTGCATCGAATTAACCCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgttagctccggaagccacgcctcaagggcac              <  2:352352‑M1/150‑44 (MQ=255)

                                                                                                                                 aggtAAGGTTCTTCGCGTTGCATCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgttagctccggaagccacgcctcaagggcac              <  1:1020166‑M1/150‑44 (MQ=255)

                                                                                                                                  ggTAAGGTTCTTCGCGTTGCATCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgttagctccggaagccacgcctcaagggcaca             <  1:191240‑M1/151‑45 (MQ=255)

                                                                                                                                  ggTAAGGTTCTTCGCGTTGCATCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgttagctccggaagccacgcctcaagggcaca             >  1:254492‑M1/1‑107 (MQ=255)

                                                                                                                                  ggTAAGGTTCTTCGCGTTGCATCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgttagctccggaagccacgcctcaagggcaca             >  1:237169‑M1/1‑107 (MQ=255)

                                                                                                                                  ggTAAGGTTCTTCGCGTTGCATCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgttagctccggaagccacgcctcaagggcaca             >  2:646898‑M1/1‑107 (MQ=255)

                                                                                                                                   gTAAGGTTCTTCGCGGTGCATCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgttagctccggaagccacgcctcaagggcacaa            <  2:532951‑M1/151‑46 (MQ=255)

                                                                                                                                       ggTTCTTCGCGTTGCATCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgttagctccggaagccacgcctcaagggcacaacctc        >  1:531566‑M1/1‑102 (MQ=255)

                                                                                                                                       ggTTCTTCGCGTTGCATCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgttagctccggaagccacgcctcaagggcacaacctc        <  1:660845‑M1/151‑50 (MQ=255)

                                                                                                                                       ggTTCTTCGCGTTGCATCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCAGtcgacttaacgcgttagctccggaagccacgcctcaagggcacaacctc        >  1:907909‑M1/1‑102 (MQ=255)

                                                                                                                                         ttcttcGCGTTGCATCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgttagctccggaagccacgcctcaagggcacaacctcca      >  2:76216‑M1/1‑100 (MQ=255)

                                                                                                                                         ttcttcGCGTTGCATCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGccgacttaacgcgttagctccggaagccacgcctcaagggcacaacctcca      <  2:102080‑M1/151‑52 (MQ=255)

                                                                                                                                           cttcGCGTTTCATCGAATTAAACCACATGCTCCACCGCTTGTTCGGGCCCCCTTCAATTCATTTGAGTTGTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgttagctccggaagccacgcctcaagggcacaacctccaag    <  1:951028‑M1/151‑54 (MQ=255)

                                                                                                                                           cttcGCGTTGCATCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgttagctccggaagccacgcctcaagggcacaacctccaag    <  1:130310‑M1/151‑54 (MQ=255)

                                                                                                                                             tcGCGTTGCATCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgttagctccggaagccacgcctcaagggcacaacctccaagtc  <  1:1043927‑M1/151‑56 (MQ=255)

                                                                                                                                             tcGCGTTGCATCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGtcgacttaacgcgttagctccggaagccacgcctcaagggcacaacctccaagtc  >  1:1033970‑M1/1‑96 (MQ=255)

GGTTGCGCTCGTTGCGGGACTTAACCCAACATCTCACGACACGAGCTGACGACAACCATGCACCACCTGTCACTTTGTCCCCCGAAGGGGAAAGCTCTATCTCTAGAGTTGTCAAAGGATGTCAAGATTTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGTCCCCGTCAATTCCTTTGAGTTTCAACCTTGCGGTCGTACTCCCCAGGCGGAGTGCTTAATGCGTTAGCTGCAGCACTAAGGGGCGGAAACCCCCTAACACTTAGC  >  CP000730/2002298‑2002589

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