BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,264,654	C→A	30.7%	A26A (GCG→GCT)	ldrA ←	toxic polypeptide, small

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,264,654	0	C	A	30.7%	224.5 / 94.7	140	A26A (GCG→GCT)	ldrA	toxic polypeptide, small
Reads supporting (aligned to +/- strand): ref base C (51/45); new base A (18/26); total (69/71)
Fisher's exact test for biased strand distribution p-value = 2.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 143 total aligned reads displayed.

GTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCAAGATTAGCCCCCG  >  CP009273/1264506‑1264802

gTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATc                                                                                                                                                      >  1:1151952/1‑149 (MQ=255)

 ttCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCg                                                                                                                                                     >  1:3098888/1‑149 (MQ=255)

   cTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCt                                                                                                                                                   <  2:1514475/149‑1 (MQ=255)

   cTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCt                                                                                                                                                   <  2:2276319/149‑1 (MQ=255)

     gCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGc                                                                                                                                                 <  1:4425266/149‑1 (MQ=255)

     gCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATAGCTGc                                                                                                                                                 <  1:3185935/149‑1 (MQ=255)

       gTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCgg                                                                                                                                               >  2:3007995/1‑149 (MQ=255)

       gTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATAGCTGCgg                                                                                                                                               <  1:1538473/149‑1 (MQ=255)

        ttAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTTCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCggg                                                                                                                                              >  1:508558/1‑148 (MQ=255)

        ttAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGt                                                                                                                                              >  2:578188/1‑149 (MQ=255)

        ttAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGt                                                                                                                                              >  2:78181/1‑149 (MQ=255)

        ttAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTAGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTACGGTTACACCACCAGCTGACAATAGCAGCGGt                                                                                                                                              >  1:2484517/1‑149 (MQ=17)

         tAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGta                                                                                                                                             <  2:3090129/149‑1 (MQ=255)

          aaGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGtaa                                                                                                                                            <  1:4120676/149‑1 (MQ=255)

          aaGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGtaa                                                                                                                                            >  1:2306724/1‑149 (MQ=255)

           aGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCCTATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGtaat                                                                                                                                           >  2:4443255/1‑149 (MQ=255)

            gACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGtaata                                                                                                                                          >  2:362563/1‑149 (MQ=255)

             aCTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATAGCTGCGGtaataa                                                                                                                                         <  1:2303330/149‑1 (MQ=255)

             aCTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATAGCTGCGGtaataa                                                                                                                                         <  1:4500510/149‑1 (MQ=255)

               ttttttCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAAtt                                                                                                                                       >  2:1786280/1‑149 (MQ=255)

                tttttCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATAGCTGCGGTAATAATTc                                                                                                                                      <  1:2125270/149‑1 (MQ=255)

                tttttCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATAGCTGCGGTAATAATTc                                                                                                                                      >  1:2771020/1‑149 (MQ=255)

                 ttttCATGGGTGCCGGATACAAATAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCTGCCTGACGCACGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATAGCTGCGGTAATAATTcc                                                                                                                                     <  2:2336348/149‑1 (MQ=255)

                      aTGGGTGCCGGATACAAAATAGGCCGCAGGCTGTTACCCCTGCGGCCTGTTTCGGGCGCATATTGCCATCACGGCATCCTGACGCCCGTTTTCACCTTACTTCCGGTTTCGACACCAGCTGACAATAGCTGCGGTAATAATTCCCGCCa                                                                                                                                <  1:4688733/149‑1 (MQ=14)

                      aTGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATAGCTGCGGTAATAATTCCCGCCa                                                                                                                                <  2:1314857/149‑1 (MQ=255)

                       tGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAg                                                                                                                               <  2:3293343/149‑1 (MQ=255)

                       tGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAg                                                                                                                               <  1:3845629/149‑1 (MQ=255)

                       tGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATAGCTGCGGTAATAATTCCCGCCAg                                                                                                                               <  1:1437723/149‑1 (MQ=255)

                       tGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATAGCTGCGGTAATAATTCCCGCCAg                                                                                                                               <  2:3007788/149‑1 (MQ=255)

                       tGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCAGGTTACGCCACCAGCTGACAATCGCTGCGGTACTAATTCCCGCCAg                                                                                                                               >  2:3972490/1‑149 (MQ=255)

                        gggTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAgg                                                                                                                              >  1:4264371/1‑149 (MQ=40)

                        gggTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAgg                                                                                                                              >  1:1543443/1‑149 (MQ=40)

                        gggTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCAAAATGCCATCACGGCAGCCAGACGCCCGTTTTCACCGTACTTCCGGTTACGCCCCCAGCTGACAATAGCTGCGGTAAAAAATACCGCCAgg                                                                                                                              >  1:1747420/1‑149 (MQ=14)

                        gggTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGGTTCGGGCGCATATTGCCATCACGGCAGCCAGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACACTAGCTGCGGTAATAATTCCCGCCAgg                                                                                                                              >  1:4136675/1‑149 (MQ=255)

                         ggTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGa                                                                                                                             <  1:215063/149‑1 (MQ=40)

                         ggTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGa                                                                                                                             <  1:215060/149‑1 (MQ=40)

                          gTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGAt                                                                                                                            <  1:2718784/149‑1 (MQ=39)

                          gTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGAt                                                                                                                            <  1:2718796/149‑1 (MQ=39)

                           tGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATc                                                                                                                           >  2:2958718/1‑149 (MQ=38)

                           tGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATc                                                                                                                           >  2:3500165/1‑149 (MQ=38)

                              cGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGt                                                                                                                        >  1:4237402/1‑149 (MQ=35)

                              cGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATAGCTGCGGTAATAATTCCCGCCAGGATCGGt                                                                                                                        >  1:4376811/1‑149 (MQ=25)

                              cGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCAACTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGAACGGt                                                                                                                        >  1:4237392/1‑149 (MQ=25)

                               ggATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGtg                                                                                                                       >  2:1455202/1‑149 (MQ=35)

                               ggATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGtg                                                                                                                       >  2:1766619/1‑149 (MQ=35)

                                  tACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATAACGGCAGCCTGACGCCCGTTTTCACCTTAATTACGGTTAAGCCACCAGCTGACAAAAGCTGCGGTCATAATTACCGCCAGGATCGGtgctg                                                                                                                    >  1:2005235/1‑149 (MQ=11)

                                   aCAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGtgctgc                                                                                                                   <  1:4439385/149‑1 (MQ=32)

                                    cAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGcc                                                                                                                  <  2:3698447/149‑1 (MQ=32)

                                     aaaaaaGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTACCGCCAGGATCGGTGCTGCCa                                                                                                                 >  1:172897/1‑149 (MQ=17)

                                     aaaaaaGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTAAGCCACCACCTGACAAAAGCTGCGGTAATAAATCCCGCCCGGATCGGTGCTGCCa                                                                                                                 >  1:3730476/1‑149 (MQ=11)

                                     aaaaaaGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTAGGGCGCATATTGCCATAACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATAGCTGCGGTAATAATTACCGCCCGGATAGGTGCTGCCa                                                                                                                 >  1:3730421/1‑149 (MQ=11)

                                      aaaaaGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAg                                                                                                                <  2:2592779/149‑1 (MQ=32)

                                      aaaaaGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAg                                                                                                                >  2:2963930/1‑149 (MQ=32)

                                      aaaaaGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAg                                                                                                                <  2:1120734/149‑1 (MQ=32)

                                      aaaaaGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAg                                                                                                                >  1:3557119/1‑149 (MQ=32)

                                      aaaaaGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGGTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAg                                                                                                                <  2:3508985/149‑1 (MQ=17)

                                       aaaaGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAgg                                                                                                               <  1:3160651/149‑1 (MQ=32)

                                         aaGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCCGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCCGGTc                                                                                                             >  1:4245232/1‑149 (MQ=14)

                                         aaGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTc                                                                                                             <  2:967848/149‑1 (MQ=31)

                                                                                                                           tCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAt                           <  2:2425833/149‑1 (MQ=11)

                                                                                                                           tCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAt                           <  2:32928/149‑1 (MQ=11)

                                                                                                                           tCCGGTTACGCCACCAGCTGACAATAGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAt                           <  2:811376/149‑1 (MQ=11)

                                                                                                                            ccGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAtt                          >  2:2212147/1‑149 (MQ=31)

                                                                                                                             cGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTg                         >  1:2041563/1‑149 (MQ=32)

                                                                                                                              ggTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGt                        <  2:1970786/149‑1 (MQ=32)

                                                                                                                              ggTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATACAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGt                        <  2:1860704/149‑1 (MQ=17)

                                                                                                                              ggTTACGCCACCAGCTGACAATAGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGt                        <  2:810298/149‑1 (MQ=17)

                                                                                                                               gTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGtt                       <  1:1140064/149‑1 (MQ=32)

                                                                                                                               gTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGtt                       <  1:1481695/149‑1 (MQ=32)

                                                                                                                               gTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGtt                       >  1:3921809/1‑149 (MQ=32)

                                                                                                                               gTTACGCCACCAGCTGACAATAGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGGTGTAATGACAACATTTTGCGGCTATTCTTGAATTGtt                       <  2:4136675/149‑1 (MQ=17)

                                                                                                                                ttACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTc                      <  2:508558/149‑1 (MQ=33)

                                                                                                                                gtACGCCACCAGCTGACAAGCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTNc                      <  2:102007/148‑1 (MQ=17)

                                                                                                                                 tACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCt                     >  1:1525966/1‑149 (MQ=34)

                                                                                                                                  aCGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTg                    >  2:1692079/1‑149 (MQ=35)

                                                                                                                                  aCGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTg                    >  2:1692022/1‑149 (MQ=35)

                                                                                                                                  aCGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTg                    <  1:635737/149‑1 (MQ=35)

                                                                                                                                   cGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTgg                   <  2:1849542/149‑1 (MQ=35)

                                                                                                                                   cGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTgg                   <  2:1849648/149‑1 (MQ=35)

                                                                                                                                    gCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGt                  >  1:3515697/1‑149 (MQ=37)

                                                                                                                                      caccaGCTGACAATAGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTc                >  2:1925409/1‑149 (MQ=38)

                                                                                                                                       accaGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCa               <  2:682102/149‑1 (MQ=38)

                                                                                                                                        ccaGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCaa              >  1:4620936/1‑149 (MQ=39)

                                                                                                                                        ccaGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCaa              >  1:3125847/1‑149 (MQ=39)

                                                                                                                                        ccaGCTGACAATAGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCaa              <  1:1811933/149‑1 (MQ=255)

                                                                                                                                        ccaGCTGACAATAGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCaa              <  1:1811876/149‑1 (MQ=255)

                                                                                                                                         caGCTGACAATAGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCAAg             <  2:128863/149‑1 (MQ=255)

                                                                                                                                         caGCTGACAATAGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTACTGGTTCAAg             >  2:3586005/1‑149 (MQ=255)

                                                                                                                                          aGCTGACAATCGCTGCGGTAATAAGTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCAAGa            <  1:4007382/149‑1 (MQ=255)

                                                                                                                                           gCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTCTTCTTGAATTGTTCTGGTTCAAGAt           >  1:3343423/1‑149 (MQ=255)

                                                                                                                                           gCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCAAGAt           >  1:954583/1‑149 (MQ=40)

                                                                                                                                           gCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCAAGAt           >  1:605520/1‑149 (MQ=40)

                                                                                                                                           gCTGACAATAGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCAAGAt           >  2:928950/1‑149 (MQ=255)

                                                                                                                                            cTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCAAGAtt          <  1:251911/149‑1 (MQ=40)

                                                                                                                                             tGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCAAGATTa         >  2:76567/1‑149 (MQ=40)

                                                                                                                                               aCAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCAAGATTAGc       >  2:2832381/1‑149 (MQ=40)

                                                                                                                                               aCAATAGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCAAGATTAGc       >  2:611214/1‑149 (MQ=255)

                                                                                                                                                cAATAGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCAAGATTAGcc      >  2:2529937/1‑149 (MQ=255)

                                                                                                                                                 aaTCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCAAGATTAGccc     >  1:4666624/1‑149 (MQ=40)

                                                                                                                                                    cGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCAAGATTAGCCCCCg  <  1:1848053/149‑1 (MQ=40)

GTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCAAGATTAGCCCCCG  >  CP009273/1264506‑1264802

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