BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP009974	1,915,245	A→G	noncoding (50/116 nt)	RPPX_08540 ←	5S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009974	1,915,245	0	A	G	89.7%	144.3 / 11.0	87	noncoding (50/116 nt)	RPPX_08540	5S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base A (6/3); new base G (39/39); total (45/42)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTGT  >  CP009974/1915094‑1915303

gTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCacta                                                            >  1:2593894/1‑151 (MQ=21)

 tCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCactac                                                           <  2:1218770/151‑1 (MQ=21)

  cAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCactgct                                                          <  2:1220578/151‑1 (MQ=21)

  cAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCactact                                                          >  2:5025421/1‑151 (MQ=21)

   aGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCActgctg                                                         <  1:1496993/151‑1 (MQ=21)

   aGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCActgctg                                                         <  2:3047724/151‑1 (MQ=21)

   aGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCActgctg                                                         >  2:899087/1‑151 (MQ=21)

    gCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGa                                                        >  2:4719881/1‑151 (MQ=21)

      tATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGt                                                      <  2:2565155/151‑1 (MQ=21)

       aTCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGtt                                                     <  2:906068/151‑1 (MQ=21)

       aTCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGtt                                                     >  1:4697500/1‑151 (MQ=21)

        tCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTc                                                    <  2:5107954/151‑1 (MQ=21)

         cAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCg                                                   <  2:1871302/151‑1 (MQ=21)

         cAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCg                                                   <  2:1871116/151‑1 (MQ=21)

         cAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCg                                                   >  1:3772321/1‑151 (MQ=21)

          aGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCgg                                                  >  2:2561087/1‑151 (MQ=21)

          aGTTCGAACCGCTATTGAAACCCGTGAAATACAAAGGCGCTGTGTGCATCGCAGACAGGGGCTGTG‑GGTATGAATCTTGAAGATGACCTACTCACACAGGGGGAAACCCCACACTACCATCGGCGATACATCGTTTCACTATTGCGCTCgg                                                  >  2:4121604/1‑151 (MQ=17)

           gTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCggg                                                 >  1:3746018/1‑151 (MQ=18)

           gTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGACTTCggg                                                 >  1:4957152/1‑151 (MQ=18)

            ttCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCggga                                                >  1:4980868/1‑151 (MQ=18)

               tAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCgggatgg                                             >  1:3209776/1‑151 (MQ=18)

                aaCCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCCGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCgggatggg                                            <  2:4569838/151‑1 (MQ=18)

                aaCCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGACGATGCATCGTTTCACTGCTGAATTCgggatggg                                            >  1:4191243/1‑151 (MQ=18)

                   ctatatTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGACATCGGGGTGGTATc                                         >  1:1346947/3‑151 (MQ=18)

                   ctatatTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTCAGTTCTGGATGGTATc                                         >  1:1347251/3‑151 (MQ=18)

                   cGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATc                                         <  2:912075/151‑1 (MQ=18)

                   cGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATc                                         <  2:4682660/151‑1 (MQ=18)

                    tgtatTTAATCCGGTGAGAGACAAAGCCGCTGTCTACATGGCGGACGGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCa                                        <  2:3772321/149‑1 (MQ=14)

                    gATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCa                                        <  2:938766/151‑1 (MQ=17)

                    gATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCa                                        >  2:265438/1‑151 (MQ=18)

                    gATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGAGCa                                        >  1:2127244/1‑151 (MQ=18)

                     atatTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAg                                       <  1:2761973/151‑1 (MQ=18)

                      tatTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAgg                                      >  1:4581941/1‑151 (MQ=18)

                       cctgcaaCCCGTGAAAGACAAAGCCCCTCTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACGCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAggt                                     <  1:3794996/146‑1 (MQ=14)

                       atTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAggt                                     >  2:417320/1‑151 (MQ=17)

                        ttGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAggtg                                    >  1:2537106/1‑151 (MQ=17)

                            aaCCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTc                                <  1:436911/151‑1 (MQ=17)

                            aaCCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGCTCAGGTGGTTc                                >  2:1519210/1‑151 (MQ=17)

                            aaCCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCAGGATGGGATCAGGTGGTTc                                >  1:3177461/1‑151 (MQ=17)

                             aCCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGTGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTcc                               >  1:1625215/1‑151 (MQ=17)

                             aCCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTcc                               <  2:4146972/151‑1 (MQ=17)

                              cccGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCa                              >  1:1860505/1‑151 (MQ=17)

                                cGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAAt                            >  1:1451994/1‑151 (MQ=17)

                                cGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAAt                            <  1:4732896/151‑1 (MQ=17)

                                 gTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACTATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATg                           >  1:2419672/1‑151 (MQ=14)

                                 gTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATg                           >  1:944160/1‑151 (MQ=17)

                                  tGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGc                          >  2:1179018/1‑151 (MQ=14)

                                   gAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGct                         >  1:4539503/1‑151 (MQ=14)

                                    aaaGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGctc                        >  1:3367321/1‑151 (MQ=14)

                                     aaGCCATAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCGCATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGctct                       <  2:1228599/151‑1 (MQ=14)

                                     aaGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGctct                       >  2:3768484/1‑151 (MQ=14)

                                       gACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTAt                     <  1:4673621/151‑1 (MQ=14)

                                       gACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTAt                     <  1:1640094/151‑1 (MQ=14)

                                        aCAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATg                    <  1:1442463/151‑1 (MQ=14)

                                        aCAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATg                    <  2:922014/151‑1 (MQ=14)

                                        aCAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATg                    >  1:5014245/1‑151 (MQ=14)

                                        aCAAAGCCCCTGTCTGCATCGCAGACAGAGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATg                    <  1:302368/151‑1 (MQ=14)

                                         cAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATgg                   <  2:611540/151‑1 (MQ=14)

                                         cAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATgg                   <  2:611481/151‑1 (MQ=14)

                                         cAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATgg                   >  1:544301/1‑151 (MQ=14)

                                         cAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTACgg                   >  1:2128609/1‑151 (MQ=14)

                                           aaGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGgtc                 <  1:230514/151‑1 (MQ=14)

                                           aaGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATAGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGgtc                 >  2:3280837/1‑151 (MQ=14)

                                            aGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGgtcg                <  2:3612983/151‑1 (MQ=14)

                                               cccTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCa             <  1:2752445/151‑1 (MQ=14)

                                               cccTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCa             <  2:779987/151‑1 (MQ=14)

                                               cccTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCa             >  1:4265153/1‑151 (MQ=14)

                                               cccTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCa             <  2:779145/151‑1 (MQ=14)

                                               cccTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCa             >  2:2114828/1‑151 (MQ=14)

                                               cacTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCa             <  2:5112398/149‑1 (MQ=14)

                                                ccTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCaa            <  1:2930093/151‑1 (MQ=14)

                                                ccTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCaa            <  2:2595089/151‑1 (MQ=11)

                                                ccTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCaa            <  2:3414013/151‑1 (MQ=14)

                                                 cTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAg           >  1:4587235/1‑151 (MQ=11)

                                                  tGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGa          >  1:1035301/1‑151 (MQ=11)

                                                  tGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGgtc                 <  2:183351/144‑1 (MQ=11)

                                                  tGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGgtc                 >  1:183351/1‑144 (MQ=11)

                                                   gTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCa             <  1:4251452/147‑1 (MQ=11)

                                                   gTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCa             >  2:4251452/1‑147 (MQ=11)

                                                   gTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGaa         >  2:3053801/1‑151 (MQ=11)

                                                    tCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGaaa        <  1:1289559/151‑1 (MQ=11)

                                                     cTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAAt       >  1:3255914/1‑151 (MQ=11)

                                                      tGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAAtt      <  2:1560741/151‑1 (MQ=11)

                                                      tGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAAtt      <  1:2561087/151‑1 (MQ=11)

                                                      tGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAAtt      >  2:2286624/1‑151 (MQ=11)

                                                        cATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCt    <  2:750341/151‑1 (MQ=11)

                                                         aTCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTg   <  2:2014424/151‑1 (MQ=11)

                                                         aTCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGCCCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTg   <  2:3581622/151‑1 (MQ=11)

                                                         aTCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTg   >  2:4766495/1‑151 (MQ=11)

                                                         aTCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTg   >  1:3711624/1‑151 (MQ=11)

                                                          tCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTGt  <  1:3535669/151‑1 (MQ=11)

                                                          tCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTGt  <  1:2883318/151‑1 (MQ=11)

                                                          tCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTGt  >  1:2012034/1‑151 (MQ=11)

                                                          tCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTGt  >  1:1793148/1‑151 (MQ=11)

GTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTGT  >  CP009974/1915094‑1915303

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