BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,382,399	A→C	8.0%	intergenic (+35/‑1)	glyX → / → glyY	tRNA‑Gly/tRNA‑Gly

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,382,399	0	A	C	8.0%	436.2 / 7.0	174	intergenic (+35/‑1)	glyX/glyY	tRNA‑Gly/tRNA‑Gly
Reads supporting (aligned to +/- strand): ref base A (101/54); new base C (3/14); total (104/70)
Fisher's exact test for biased strand distribution p-value = 3.09e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 177 total aligned reads displayed.

AAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACT  >  CP009273/4382252‑4382547

aaaaTTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAATTGATTTAATTCACACACCACCCAAga                                                                                                                                                    >  2:186558/1‑149 (MQ=17)

  aaTTTGAAAAGTGCTGCAAATCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGAACGAACTTTCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATATG‑AAAGTGATATAAGGCACACACCACCCAAGCgg                                                                                                                                                  >  2:743749/1‑150 (MQ=18)

   aTTTGAAAAGTGCTCCAAAGCACAGACCACCCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCggg                                                                                                                                                 <  1:3401262/150‑1 (MQ=255)

    tttGAAAAGTGCTTCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGa                                                                                                                                                <  2:1592993/150‑1 (MQ=255)

    tttGAAAAGTGCTGCAAAGCACTGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGa                                                                                                                                                <  2:1023603/150‑1 (MQ=255)

     ttGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTCGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTGG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGaa                                                                                                                                               >  1:2023690/1‑150 (MQ=255)

      tGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCACGCGGGAAt                                                                                                                                              >  1:857872/1‑150 (MQ=255)

      tGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAAt                                                                                                                                              <  2:1973268/150‑1 (MQ=255)

      tGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAAt                                                                                                                                              >  1:3083553/1‑150 (MQ=255)

        aaaaGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAg                                                                                                                                            >  2:2893645/1‑150 (MQ=255)

           aGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATATG‑AAAGTGCTGTAAGGCACACACCACCCAAGCGCGAATAGcac                                                                                                                                         >  1:667638/1‑148 (MQ=255)

            gTGCTGCACATCACTGACCACCCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCa                                                                                                                                        <  1:127898/150‑1 (MQ=255)

             tgctgcAAAGCACAGACCACCCAAGCGGGATTAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAg                                                                                                                                       <  2:6455/150‑1 (MQ=255)

               ctgcAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCACGCGGGAATAGCTCAggt                                                                                                                                     >  1:2557223/1‑148 (MQ=255)

                  cAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGGTCGAGTCTCGTTTCCCGCTCCAAAATATT‑AAAGAGATATAATAAACACACCGACCAAACGGGAATAGCACAGTTgtt                                                                                                                                  >  1:1996185/1‑148 (MQ=14)

                  aacaccacaGAACTACCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGt                                                                                                                                  <  2:3921658/145‑1 (MQ=11)

                   aaaGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTa                                                                                                                                 >  1:3300953/1‑150 (MQ=255)

                       cacaGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGATGTAAGGCGCAGACCACCCAAGCGGGAATAGCTGAGTTGGTAGAGc                                                                                                                             >  2:1961501/1‑150 (MQ=255)

                       cacaGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATGTT‑AAAGTACTGTAATGCACAGACCGCCCCAGAGGGAATAGCGCAGTTGTTAGAGc                                                                                                                             >  2:553203/1‑150 (MQ=14)

                         caGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAATGCACAGACCACCCCATCGGGAATAGCTCAGTTGGTAGAGAAc                                                                                                                           >  2:2807663/1‑150 (MQ=255)

                              ctcccAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCt                                                                                                                      <  2:63587/148‑1 (MQ=255)

                              caccCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGTGCACGAccg                                                                                                                      >  1:13980/1‑149 (MQ=255)

                                cccAAGCGGGACTAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGAGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTGGAGCACGACCTTg                                                                                                                    >  2:422073/1‑150 (MQ=255)

                                  acaGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGAAAGGCACAGACCACCCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGcc                                                                                                                  <  2:2622054/148‑1 (MQ=255)

                                   acgCGGGAATAGCTCAGTTGGTTGTGCTGGTCAGTGCCAAGGTCGGGGTCGCGAGTTCGCGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCa                                                                                                                 <  2:303667/148‑1 (MQ=11)

                                   aaGTGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCGCGGTTGGTAGAGCCCGACCTTGCCa                                                                                                                 >  2:2981561/1‑150 (MQ=255)

                                   aaGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGTATAGCTCAGTTGTTAGAGCACGACCTAGCCa                                                                                                                 >  2:2451385/1‑150 (MQ=255)

                                     gCGGGAATAGCTCAGTTGGTAGAGCACGCCCTTGCCTAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAg                                                                                                               <  1:1752108/150‑1 (MQ=255)

                                      cGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAgg                                                                                                              <  1:3683682/150‑1 (MQ=38)

                                      cGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAgg                                                                                                              >  2:1142159/1‑150 (MQ=38)

                                       gggAATAGCTCAGTTGGTAGAGCTCGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGAACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCACGGt                                                                                                             >  2:1080129/1‑150 (MQ=255)

                                       gggAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCCGTTGGTAGAGCACGACCGTGCCAAGGt                                                                                                             >  1:1547356/1‑150 (MQ=255)

                                         gAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCg                                                                                                           >  1:3642626/1‑150 (MQ=38)

                                          aaTAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgg                                                                                                          >  1:3351677/1‑150 (MQ=38)

                                           aTAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCggg                                                                                                         <  1:4084811/150‑1 (MQ=38)

                                             aGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAAAAGCTCAGTTGGTAGAGCACGACCATGCCAAGGTCGGGGt                                                                                                       >  1:45477/1‑150 (MQ=255)

                                              gCTCAGTGTGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGCGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTc                                                                                                      <  2:543623/150‑1 (MQ=255)

                                               cTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAATGAACAGACCACCCAAGCGGGAATATGTCAGTTGGTAGAGCACGACCGTGACAAGTTCggcggcg                                                                                                     >  2:1977628/1‑145 (MQ=11)

                                                 cAGTTCGTAGAGGCCGACCTGCCCAAGGTTGGGTTCGCGAGTTCGAGTCTCGTTTCGCCCTAGAAATTCTC‑AAAGTTCTGTAAGACAAAGACGACCCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTcgcg                                                                                                   <  2:3154030/150‑1 (MQ=11)

                                                  aGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGa                                                                                                  >  2:176576/1‑150 (MQ=38)

                                                        tAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAATCGGGAATAGCTCAGTTGGTAGAGCACGACCGTGCCAAGGTCGGGGTCGCGAGTTCGa                                                                                            >  1:1012573/1‑150 (MQ=255)

                                                        tAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGa                                                                                            >  1:1012322/1‑150 (MQ=38)

                                                        tAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGa                                                                                            >  1:1891141/1‑150 (MQ=38)

                                                          gtgCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGCTTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGt                                                                                          <  1:1743330/148‑1 (MQ=255)

                                                             cACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtctc                                                                                       <  1:2717029/150‑1 (MQ=38)

                                                                 aCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGttt                                                                                   >  2:1782318/1‑150 (MQ=38)

                                                                   cTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTcc                                                                                 >  1:677358/1‑150 (MQ=38)

                                                                     tGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTACGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCg                                                                               <  1:3948257/150‑1 (MQ=38)

                                                                      gCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGc                                                                              >  1:1508397/1‑150 (MQ=38)

                                                                        cAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTc                                                                            >  2:571550/1‑150 (MQ=38)

                                                                         aaGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTcc                                                                           >  2:1482223/1‑150 (MQ=38)

                                                                              cGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAtt                                                                      <  1:3237534/150‑1 (MQ=40)

                                                                              cGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAtt                                                                      <  1:387158/150‑1 (MQ=40)

                                                                               ggggTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAttc                                                                     >  2:3758034/1‑150 (MQ=40)

                                                                                 ggTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAttctt                                                                   >  1:3886983/1‑150 (MQ=255)

                                                                                        aGTTCGAGTATCGTTTCCCGCTCCAAAATTTC‑AAAGTGCTGTACGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCaa                                                            <  2:237123/150‑1 (MQ=255)

                                                                                               gTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAata                                                     >  2:2670300/1‑150 (MQ=255)

                                                                                               gTCTCGTTTCCCGCTCCAAAATTGG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAata                                                     <  2:3340372/150‑1 (MQ=255)

                                                                                                   cGTTTCCCGCTCCAAAATTTT‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCca                                                 <  2:492610/150‑1 (MQ=255)

                                                                                                   cGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCca                                                 >  2:375993/1‑150 (MQ=255)

                                                                                                   cGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCca                                                 >  2:2213395/1‑150 (MQ=255)

                                                                                                      ttCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAg                                              >  2:2179065/1‑150 (MQ=255)

                                                                                                      ttCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAg                                              >  1:469957/1‑150 (MQ=255)

                                                                                                      ttCCCGCTACAAATTTTT‑AAAGTGCTGTAAGTCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAg                                              <  1:2201434/150‑1 (MQ=255)

                                                                                                         ccGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGa                                           >  2:3352127/1‑150 (MQ=255)

                                                                                                         ccGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGa                                           >  2:3352098/1‑150 (MQ=255)

                                                                                                          cggtagaaaaTTTT‑AAAGTGCTGTAAGTCACTGACCTCCTAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGAc                                          <  2:2532286/144‑1 (MQ=18)

                                                                                                            cTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGAcgc                                        >  2:88977/1‑150 (MQ=255)

                                                                                                            cTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGAcgc                                        >  1:3971979/1‑150 (MQ=255)

                                                                                                             tCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGAcgcg                                       >  1:507174/1‑150 (MQ=255)

                                                                                                              ccAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGa                                      >  1:3845120/1‑150 (MQ=255)

                                                                                                               cAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGAt                                     >  2:578921/1‑150 (MQ=255)

                                                                                                               acaaaTTTGTAAATTTCTTTAACTCACAGACCACCAACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGa                                      <  2:1152224/148‑1 (MQ=17)

                                                                                                                  aaTTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCg                                  >  1:1236214/1‑150 (MQ=255)

                                                                                                                    tttG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGtt                                >  2:3247202/1‑150 (MQ=255)

                                                                                                                    tgtG‑ACAGTGCTGTAAGGCACAGACCACCCCAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGtt                                <  1:2664689/148‑1 (MQ=255)

                                                                                                                      tG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTAt                              >  2:1671050/1‑150 (MQ=255)

                                                                                                                      tG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTAt                              >  2:209121/1‑150 (MQ=255)

                                                                                                                         aaatttcTGTACGGCAAAGACGACGCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTg                            <  2:3173434/144‑1 (MQ=21)

                                                                                                                         aaaGTGCTGTAAGGCACAGACCACCCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTg                            <  1:2622869/150‑1 (MQ=255)

                                                                                                                           aGTGCTGCAACGCACAGACCCCCCCAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCt                          <  1:3662591/150‑1 (MQ=255)

                                                                                                                            gTGCTGTAAGGCACAGACCACCCACGCGGGAAGAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTg                         <  1:1893655/150‑1 (MQ=255)

                                                                                                                              gCTGTAACGCACAGACGACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCGTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGt                       <  2:2626658/150‑1 (MQ=255)

                                                                                                                                   aaGTCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTttg                  <  2:3238462/150‑1 (MQ=255)

                                                                                                                                   aaGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTttg                  >  1:2445510/1‑150 (MQ=255)

                                                                                                                                       cacaGACCACCTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTtgttgt              <  2:175654/150‑1 (MQ=255)

                                                                                                                                       cacaGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTtgttgt              >  2:3066533/1‑150 (MQ=255)

                                                                                                                                       cacaGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAATTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTtgttgt              >  1:707844/1‑150 (MQ=255)

                                                                                                                                        acaGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGtc             >  1:223980/1‑150 (MQ=255)

                                                                                                                                         caGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGtct            >  2:2839550/1‑150 (MQ=255)

                                                                                                                                          aGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGtctc           >  1:2248338/1‑150 (MQ=255)

                                                                                                                                          aGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGtctc           >  2:6753/1‑150 (MQ=255)

                                                                                                                                           gACCACCCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTctct          <  2:1440271/150‑1 (MQ=255)

                                                                                                                                           gACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTctct          <  2:3449628/150‑1 (MQ=255)

                                                                                                                                           gACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTctct          <  1:949684/150‑1 (MQ=255)

                                                                                                                                             ccaccCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGa        <  1:3675131/150‑1 (MQ=255)

                                                                                                                                               accCAGGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAATTATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACa      <  1:1141941/150‑1 (MQ=255)

                                                                                                                                               accCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACa      <  1:1142204/150‑1 (MQ=255)

                                                                                                                                                   aaGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAAct  >  1:3772805/1‑150 (MQ=255)

                                                                                                                                                   aaGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAAct  >  1:2855359/1‑150 (MQ=255)

AAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACT  >  CP009273/4382252‑4382547

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: `‑`