BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP009974	1,915,172	3 bp→GTT	intergenic (+61/+5)	RPPX_08535 → / ← RPPX_08540	ferrioxamine B receptor/5S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009974	1,915,172	0	A	G	83.5%	503.5 / 143.4	213	intergenic (+61/+7)	RPPX_08535/RPPX_08540	ferrioxamine B receptor/5S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (73/104); minor base . (17/18); total (91/122)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01
*	CP009974	1,915,173	0	A	T	100.0%	629.9 / NA	213	intergenic (+62/+6)	RPPX_08535/RPPX_08540	ferrioxamine B receptor/5S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base A (0/0); new base T (91/122); total (91/122)
*	CP009974	1,915,174	0	A	T	99.5%	626.7 / ‑3.1	213	intergenic (+63/+5)	RPPX_08535/RPPX_08540	ferrioxamine B receptor/5S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base A (0/0); major base T (90/122); minor base . (1/0); total (91/122)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 213 total aligned reads displayed.

CCTCATGTGCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTGTT  >  CP009974/1915035‑1915304

|||

ccTCGTGTGCCGGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGGTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGGCAGGGGCTTTGGG‑TTTGAATCTTGACg                                                                                                                        <  2:5146054/151‑1 (MQ=255)

  tCATGTGCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGAt                                                                                                                      >  1:764299/1‑151 (MQ=255)

   cATGTGCCAGCCTGGACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATg                                                                                                                     <  2:6338603/151‑1 (MQ=255)

     tgtgCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGAc                                                                                                                   >  1:3728235/1‑151 (MQ=255)

     tgtgCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGAc                                                                                                                   >  2:429056/1‑151 (MQ=255)

     tgtgCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGAc                                                                                                                   <  2:1855874/151‑1 (MQ=255)

     tgtgCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGAc                                                                                                                   <  1:6461780/151‑1 (MQ=255)

      gtgCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGAcc                                                                                                                  >  2:6383605/1‑151 (MQ=255)

      gtgCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGAcc                                                                                                                  >  1:1951348/1‑151 (MQ=255)

      gtgCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGAcc                                                                                                                  >  1:2322556/1‑151 (MQ=255)

        gCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTa                                                                                                                >  1:3413225/1‑151 (MQ=255)

         ccAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTAc                                                                                                               <  1:1055153/151‑1 (MQ=255)

         ccAGCCTGAAATACTGCTAGATGGGTGAAGAGCGCAAAGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAGCCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTAc                                                                                                               <  2:613192/151‑1 (MQ=255)

           aGCCTGAACTACTGGTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTActc                                                                                                             <  2:4742791/151‑1 (MQ=255)

              cTGCACTACTGCTACATGGGGGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTca                                                                                                          <  2:1912399/151‑1 (MQ=255)

                gAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTcaca                                                                                                        <  2:3485903/151‑1 (MQ=255)

                  actactGCTAGATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATg                                                                                                      <  2:2863024/151‑1 (MQ=255)

                  actactGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATg                                                                                                      >  1:2305638/1‑151 (MQ=39)

                  actactGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATg                                                                                                      <  2:2310286/151‑1 (MQ=255)

                  actactGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCGGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATg                                                                                                      <  1:5822861/151‑1 (MQ=255)

                   ctactGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATgg                                                                                                     >  2:1784643/1‑151 (MQ=255)

                     actGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTTTAACCGTTATTGAAACCCGTGAAAGACAAAGACCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATgggg                                                                                                   >  1:2565102/1‑151 (MQ=38)

                     actGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATgggg                                                                                                   <  1:5883831/151‑1 (MQ=255)

                         cTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAAc                                                                                               >  2:3818044/1‑151 (MQ=38)

                         cTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAAc                                                                                               >  2:5492926/1‑151 (MQ=38)

                          tACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAAcc                                                                                              >  1:4354434/1‑151 (MQ=38)

                           aCATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAAccc                                                                                             <  1:3491168/151‑1 (MQ=38)

                            cATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAAcccc                                                                                            >  1:5971423/1‑151 (MQ=38)

                            cATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAAcccc                                                                                            <  1:2190191/151‑1 (MQ=38)

                             aTGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCca                                                                                           >  2:4012177/1‑151 (MQ=38)

                              tGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCcac                                                                                          <  2:790948/151‑1 (MQ=38)

                              tGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCcac                                                                                          >  1:1613302/1‑151 (MQ=38)

                               gggTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCcaca                                                                                         >  1:2230223/1‑151 (MQ=38)

                                    aaGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTAcc                                                                                    >  2:6070627/1‑151 (MQ=37)

                                      gagCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCAt                                                                                  <  1:803489/151‑1 (MQ=37)

                                        gcgcAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCg                                                                                <  1:2417711/151‑1 (MQ=37)

                                          gcAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGc                                                                              >  1:3692711/1‑151 (MQ=37)

                                          gcAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGc                                                                              >  2:4244009/1‑151 (MQ=37)

                                           cAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCg                                                                             <  2:4967376/151‑1 (MQ=37)

                                             aTGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGAt                                                                           <  1:4583526/151‑1 (MQ=37)

                                              tGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATg                                                                          <  2:5279690/151‑1 (MQ=37)

                                                tAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCa                                                                        <  2:5773183/151‑1 (MQ=25)

                                                   gcgcTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCg                                                                     <  2:2444398/151‑1 (MQ=25)

                                                    cgcTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGt                                                                    >  2:3033395/1‑151 (MQ=25)

                                                      cTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGttt                                                                  <  1:2704437/151‑1 (MQ=25)

                                                      cTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGttt                                                                  >  1:392979/1‑151 (MQ=25)

                                                        aCGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGc                                                                         >  2:5317814/1‑142 (MQ=25)

                                                         cGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCac                                                               <  1:891486/151‑1 (MQ=25)

                                                          ggTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCact                                                              <  1:5177099/151‑1 (MQ=21)

                                                          ggTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCact                                                              >  2:148785/1‑151 (MQ=25)

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

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

                                                                 tATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGt                                                       >  2:3449864/1‑151 (MQ=21)

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

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

                                                                   tCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTc                                                     >  1:4213362/1‑151 (MQ=21)

                                                                    cAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCg                                                    >  1:2968509/1‑151 (MQ=18)

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

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

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

                                                                        tCTAACCGATATTGAAACCCGTGAAAGACAACGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCgggat                                                <  2:4636445/151‑1 (MQ=18)

                                                                        tCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCgggat                                                <  2:2322556/151‑1 (MQ=18)

                                                                         cTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCgggatg                                               <  1:1178540/151‑1 (MQ=18)

                                                                         cTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCgggatg                                               <  1:850386/151‑1 (MQ=18)

                                                                          tAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCgggatgg                                              <  2:6111055/151‑1 (MQ=18)

                                                                           aaCCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCgggatggg                                             <  2:2029274/151‑1 (MQ=18)

                                                                           aaCCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCgggatggg                                             <  2:5321699/151‑1 (MQ=18)

                                                                           aaCCGATATAGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCgggatggg                                             <  2:2204767/151‑1 (MQ=17)

                                                                            aCCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCgggatggga                                            <  2:4213362/151‑1 (MQ=18)

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

                                                                                  atTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGCGGAAACCCCACACTACCATCGGCGATGAATCATTTCACTACTGAGTTCaggatgggacccagt                                      >  2:6441393/1‑145 (MQ=17)

                                                                                   ttGAATCCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAggtg                                     <  2:781982/151‑1 (MQ=17)

                                                                                    tGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAggtgg                                    <  2:5175107/151‑1 (MQ=17)

                                                                                     gAAACCCGGGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAggtggt                                   <  2:5070526/151‑1 (MQ=17)

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

                                                                                       aaCCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTc                                 >  2:2094405/1‑151 (MQ=17)

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

                                                                                        aCCCGTTAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTcc                                <  1:2436106/151‑1 (MQ=14)

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

                                                                                          ccGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCaa                              <  1:147716/151‑1 (MQ=17)

                                                                                           cGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAAt                             <  2:876776/151‑1 (MQ=17)

                                                                                             tGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGc                           <  2:195283/151‑1 (MQ=14)

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

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

                                                                                                 aGACAAAGCCCCTGTCAGCATCGCAGACAGGGGCTTTGGG‑TTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTa                       <  2:1831899/151‑1 (MQ=14)

                                                                                                     aaaGCCCCTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGgt                   >  2:3310138/1‑151 (MQ=14)

                                                                                                      aggCCCGTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGgtc                  <  2:3758430/149‑1 (MQ=14)

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

                                                                                                      aaGCCCCTGTCTGCATCGCAGACAGGGG‑CTTGGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGgtc                  <  1:6435073/151‑1 (MQ=14)

                                                                                                          cccTTTCTGCATCGCAGACAGGGGCTTTG‑GGTTTTAATCTTGACGATGACCTACTCTCACATGGGTAAACCCCACACTACCATCGGCGATGCATCGGTTCACTGCTGAGTTCGGGATGTGATCAGGTGGTTCCAATGCTCTATGATCGTCa              >  2:512003/1‑151 (MQ=11)

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

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

                                                                                                            cTGTCTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAg            >  2:3667526/1‑151 (MQ=14)

                                                                                                                cTGCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAAt        <  2:72938/151‑1 (MQ=11)

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

                                                                                                                  gCATCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTc      <  1:1853456/151‑1 (MQ=11)

                                                                                                                     tCGCAGACAGGGGCTTTG‑GGTTTGAATCTTGACGATGCCCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTGt   <  2:438245/151‑1 (MQ=11)

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

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

                                                                                                                      cGCAGACAGGGGCTTTGGGTT‑TGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTGTg  >  1:1756024‑M2/1‑150 (MQ=255)

|||

CCTCATGTGCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTGTT  >  CP009974/1915035‑1915304

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