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	78.9%	151.6 / 52.4	126	noncoding (50/116 nt)	RPPX_08540	5S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base A (11/15); new base G (37/60); total (51/75)
Fisher's exact test for biased strand distribution p-value = 8.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 135 total aligned reads displayed.

GTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTGTT  >  CP009974/1915094‑1915304

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

   aGCTATGAGTTCTAAGCGTTATTGAAGGCTGTGAAAAACAGAGCCCCGGTCGGAATCGGAGACGGGGGCTTTG‑GGTTTGAATGTAGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGGTTCActgctg                                                          <  2:4197779/151‑1 (MQ=9)

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

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

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

    gCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGa                                                         <  1:1454328/151‑1 (MQ=21)

     cTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAg                                                        >  1:491433/1‑151 (MQ=21)

     cTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCAGCGCTGCAGCGTTTCACTGCTGat                                                        >  1:980741/1‑150 (MQ=21)

     cTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAg                                                        >  1:4883493/1‑151 (MQ=21)

     cTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAg                                                        <  2:1108213/151‑1 (MQ=21)

      tATCAGTTCTAGCCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTAAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGt                                                       <  2:2713782/151‑1 (MQ=18)

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

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

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

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

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

        tCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTAAGTTc                                                     >  1:6991114/1‑151 (MQ=21)

        tCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTc                                                     >  1:6953564/1‑151 (MQ=21)

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

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

          aGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGGTTCACTGCTGAGTTCgg                                                   >  2:6254836/1‑151 (MQ=21)

            ttCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCggga                                                 <  2:3975516/151‑1 (MQ=18)

            ttCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCggga                                                 >  2:7121044/1‑151 (MQ=18)

             tcaaaCCGATGTTGAAACCCGGGAACGCCAAAGCCCCTGTCTGCATCGCAGACAAGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCgggat                                                <  2:6256954/148‑1 (MQ=17)

             tCTAACCTATATTACCAACCGTGAAACACAACGCCCCTCTCTGCACCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCgggat                                                <  2:1581458/151‑1 (MQ=14)

             tCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTCCTGAGTTCgggat                                                >  2:1747699/1‑151 (MQ=18)

              cTAACCGATATTGAAACCCGTGAAAGACAAGGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCgggatg                                               <  2:5994912/151‑1 (MQ=18)

              cTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCgggatg                                               <  2:7140519/151‑1 (MQ=18)

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

                aaacGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGAATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCgggatggg                                             <  2:5497732/148‑1 (MQ=17)

                aaCCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCgggatggg                                             >  1:2188406/1‑151 (MQ=18)

                aaCCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCgggatggg                                             >  1:2188676/1‑151 (MQ=18)

                aaCCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCgggatggg                                             >  2:4499325/1‑151 (MQ=18)

                aaCCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCgggatggg                                             >  2:2872830/1‑151 (MQ=18)

                aaCCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTCGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCgggatggg                                             <  2:6930614/151‑1 (MQ=18)

                  ccGATATTGCAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCgggatgggat                                           <  2:2998471/151‑1 (MQ=17)

                  ccGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCgggatgggat                                           >  1:1192644/1‑151 (MQ=18)

                  ccGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCgggatgggat                                           <  2:2984508/151‑1 (MQ=18)

                  ccGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCgggatgggat                                           <  1:3430910/151‑1 (MQ=18)

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

                    gATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCa                                         <  1:4499325/151‑1 (MQ=17)

                    gATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCa                                         <  2:3720159/151‑1 (MQ=18)

                     atatTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAg                                        <  2:1658320/151‑1 (MQ=17)

                       atTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAggt                                      <  1:3404850/151‑1 (MQ=17)

                       atTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAggt                                      <  1:6378460/151‑1 (MQ=17)

                       atTGAAACCCGAAAAAGAAAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GATATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAggt                                      <  2:5172841/151‑1 (MQ=17)

                         tgactcccGTGAGAGACAACGCCCCTGTCTGCATCGAGGTCAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAggtgg                                    <  2:142392/146‑1 (MQ=14)

                         tGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAggtgg                                    <  2:2869926/151‑1 (MQ=17)

                         tGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAggtgg                                    <  1:787506/151‑1 (MQ=17)

                         tGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAggtgg                                    >  2:2789582/1‑151 (MQ=17)

                          gAAGCCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAggtggt                                   <  2:3265799/151‑1 (MQ=17)

                           actaccGTGAAAGACAACCCCCCTCTCTGCATCCCAGACAGGGTCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGtt                                  <  2:2188676/147‑1 (MQ=14)

                           aaaCCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGtt                                  <  2:2231844/151‑1 (MQ=17)

                           aaaCCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGtt                                  <  2:3083219/151‑1 (MQ=17)

                           aaaCCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCGggtgggc                                  >  1:1916499/1‑149 (MQ=17)

                            aaCCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTc                                 <  1:539679/151‑1 (MQ=17)

                            aaCCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGGTTGAATCTTGACGATGTCCTACTCTAACATGGGGAAACCCCACACTACCATCGGCGATGTATCGATTCACTTCTGAGTTCGGGATGGGATTtggtgggta                                 >  2:5130355/1‑148 (MQ=17)

                              cccGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGATTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCa                               <  2:556779/151‑1 (MQ=14)

                               ccGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCaa                              <  2:5849564/151‑1 (MQ=14)

                                cGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAAt                             <  2:2587539/151‑1 (MQ=14)

                                cGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAAt                             <  2:6231927/151‑1 (MQ=14)

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

                                 gTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATg                            <  2:2307394/151‑1 (MQ=17)

                                    aaaGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATGTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGctc                         <  2:6953564/151‑1 (MQ=14)

                                     aaGACAACGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGctct                        <  2:5022016/151‑1 (MQ=14)

                                     aaGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGctct                        <  2:1974178/151‑1 (MQ=14)

                                     aaGACAAAGCCCCTGTCTGCATCGCAGACAGGGGATTTG‑GGTTTGAATCTTGACGATGACATACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGctct                        <  2:5023501/151‑1 (MQ=14)

                                      aGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTGTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTa                       <  2:4346491/151‑1 (MQ=14)

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

                                        aCAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATg                     <  2:580122/151‑1 (MQ=14)

                                        aCAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGACGGGATCAGGTGGTTCCAATGctctctt                     >  1:422386/1‑148 (MQ=14)

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

                                          aaaGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGgt                   >  2:2523596/1‑151 (MQ=14)

                                           aaGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGgtc                  <  2:5173765/151‑1 (MQ=14)

                                           aaGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGGTCAGGTGGTTCCAATGCTCTATGgtc                  >  1:2469546/1‑151 (MQ=14)

                                           aaGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGgtc                  >  1:5670969/1‑151 (MQ=14)

                                           aaGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGgtc                  >  2:595048/1‑151 (MQ=14)

                                           aaGCCCCTGTCGGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGgtc                  <  2:4298150/151‑1 (MQ=14)

                                             gCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGgtcgt                >  1:4444770/1‑151 (MQ=14)

                                              gcccTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGgtcgtc               <  1:6828478/150‑1 (MQ=14)

                                              ccccTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGggcgtc               >  1:5949857/1‑151 (MQ=14)

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

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

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

                                               cccTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCAGCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCa              <  2:6079239/151‑1 (MQ=14)

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

                                                ccTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCaa             >  1:4755894/1‑151 (MQ=11)

                                                 cTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAg            <  2:5361010/151‑1 (MQ=14)

                                                 cTGTCTGCATCGCAGAAAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAg            <  1:731878/151‑1 (MQ=11)

                                                   gTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGaa          <  2:2836090/151‑1 (MQ=11)

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

                                                    tCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGaaa         <  1:672825/151‑1 (MQ=11)

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

                                                      tGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCgggatgg                                              <  1:308300/112‑1 (MQ=11)

                                                      acaagagCAGACACGGGTTTTGCGATATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTGCAATGCTCTATGGTCGTCAAGAAAt        >  2:2760952/7‑151 (MQ=255)

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

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

                                                         aTCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTg    <  2:4889081/151‑1 (MQ=11)

                                                          tCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTGt   <  2:4464425/151‑1 (MQ=11)

                                                            gCAGACAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTGtt  >  2:6865096/1‑151 (MQ=31)

GTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTGTT  >  CP009974/1915094‑1915304

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