BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,019,203	G→T	7.5%	E27* (GAG→TAG)	guaD →	guanine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,019,203	0	G	T	7.5%	691.7 / 21.5	187	E27* (GAG→TAG)	guaD	guanine deaminase
Reads supporting (aligned to +/- strand): ref base G (79/94); new base T (9/5); total (88/99)
Fisher's exact test for biased strand distribution p-value = 2.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 193 total aligned reads displayed.

CTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTAC  >  CP009273/3019056‑3019350

cTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAGCCCGGAag                                                                                                                                                     <  2:2030885/148‑1 (MQ=255)

cTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAag                                                                                                                                                     <  2:1088752/148‑1 (MQ=255)

  aCCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAagag                                                                                                                                                   >  1:2416219/1‑148 (MQ=255)

  aCCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAagag                                                                                                                                                   >  1:2226833/1‑148 (MQ=255)

    cGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGAt                                                                                                                                                 >  2:662711/1‑148 (MQ=255)

      aCAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTg                                                                                                                                               <  2:1888933/148‑1 (MQ=255)

      aCAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTg                                                                                                                                               <  2:2856208/148‑1 (MQ=255)

       cAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGc                                                                                                                                              <  1:1500649/148‑1 (MQ=255)

       cAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGc                                                                                                                                              >  2:1640348/1‑148 (MQ=255)

           aaaaCGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCctct                                                                                                                                          >  1:2925271/1‑148 (MQ=255)

           aaaaCGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCctct                                                                                                                                          >  2:775907/1‑148 (MQ=255)

               cGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTctgcgc                                                                                                                                      <  2:2245853/148‑1 (MQ=255)

                 ttCTGCGTGGTTTTTCCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCTATGATGTCATGAGAACACACGTTAAAAGCGGTACGAGGCTGTTTTTTTGATGTCACCCGTACGATCGATTACCCGGAAGAGATTGCCTctgcgctg                                                                                                                                    <  2:3093282/148‑1 (MQ=255)

                   cTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTctgcgctgcg                                                                                                                                  <  2:1827955/148‑1 (MQ=255)

                    tGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCgg                                                                                                                                 <  1:2684792/148‑1 (MQ=255)

                        tGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTa                                                                                                                             <  1:3115650/148‑1 (MQ=255)

                            aTTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGa                                                                                                                         >  1:2038124/1‑148 (MQ=255)

                             ttACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAg                                                                                                                        <  1:1389112/148‑1 (MQ=255)

                             ttACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAg                                                                                                                        >  1:997340/1‑148 (MQ=255)

                              tACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAgg                                                                                                                       <  1:2006018/148‑1 (MQ=255)

                              tACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAgg                                                                                                                       <  1:500533/148‑1 (MQ=255)

                                  tCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATgg                                                                                                                   >  1:324346/1‑148 (MQ=255)

                                     gCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGttt                                                                                                                <  2:459500/148‑1 (MQ=255)

                                      cGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTtta                                                                                                               <  1:1595483/148‑1 (MQ=255)

                                         aaGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTttatta                                                                                                            <  1:3195974/148‑1 (MQ=255)

                                         aaGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTttatta                                                                                                            <  1:1084876/148‑1 (MQ=255)

                                         aaGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTttatta                                                                                                            >  2:1324983/1‑148 (MQ=255)

                                         aaGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTttatta                                                                                                            <  1:1033121/148‑1 (MQ=255)

                                               tGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCAtt                                                                                                      >  1:594257/1‑148 (MQ=255)

                                                 gATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTaa                                                                                                    >  2:1657883/1‑148 (MQ=255)

                                                  aTTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTaaa                                                                                                   >  2:2465097/1‑148 (MQ=255)

                                                   ttAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAAc                                                                                                  <  2:1072131/148‑1 (MQ=255)

                                                    tAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACa                                                                                                 >  1:1426913/1‑148 (MQ=255)

                                                     aaCAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAg                                                                                                <  2:3220226/148‑1 (MQ=255)

                                                      aCAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAgg                                                                                               <  1:1052192/148‑1 (MQ=255)

                                                        aGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGa                                                                                             <  1:580234/148‑1 (MQ=255)

                                                        aGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGa                                                                                             >  1:889036/1‑148 (MQ=255)

                                                         gTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGaa                                                                                            >  1:2706215/1‑148 (MQ=255)

                                                                 gTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGaa                                                                                    <  1:1924499/148‑1 (MQ=255)

                                                                  tCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGaag                                                                                   >  2:1968002/1‑147 (MQ=255)

                                                                  tCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAAt                                                                                   >  2:299841/1‑148 (MQ=255)

                                                                  tCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAAt                                                                                   <  2:1234356/148‑1 (MQ=255)

                                                                    aaTGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATgg                                                                                 <  2:1313033/148‑1 (MQ=255)

                                                                        atgTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTg                                                                             >  2:1387077/1‑148 (MQ=255)

                                                                        atgTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTg                                                                             >  1:1665095/1‑148 (MQ=255)

                                                                        atgTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTg                                                                             >  1:414394/1‑148 (MQ=255)

                                                                         tttCAGGAGAACACACGTTAAAAGGTTTACGAGGCATTTTTATTGATGTCACCCTTACGATCTATAACCCTGAAGAGATTGACTCTGCGCTGCGGTTTATTGAGGTTGGTTTATTACTCATTAATGAGGGAAAAGTGTAATGGATTgg                                                                            <  2:3077297/146‑1 (MQ=255)

                                                                              ggAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAAt                                                                       <  1:2561227/148‑1 (MQ=255)

                                                                               gagaACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGACTg                                                                      >  2:578865/1‑148 (MQ=255)

                                                                               gagaACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATg                                                                      >  2:176649/1‑148 (MQ=255)

                                                                                 gaACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATggg                                                                    >  1:99217/1‑148 (MQ=255)

                                                                                  aaCACACGTTATAAGCGGTACGAGGAAGTTTTATTGATGTCTCCCGTACGATCGATAACCCTGAAGAGATGGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGa                                                                   <  1:176188/148‑1 (MQ=255)

                                                                                    cacacGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGaaa                                                                 >  2:1433152/1‑148 (MQ=255)

                                                                                          ttAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGaa                                                           >  2:403521/1‑148 (MQ=255)

                                                                                          ttAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGaa                                                           >  2:2915661/1‑148 (MQ=255)

                                                                                           tAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGaaa                                                          >  2:372279/1‑148 (MQ=255)

                                                                                            aaaaGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAg                                                         <  2:2846918/148‑1 (MQ=255)

                                                                                              aaGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCa                                                       <  2:3230936/148‑1 (MQ=255)

                                                                                              aaGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCa                                                       <  2:82444/148‑1 (MQ=255)

                                                                                                gCGGTACGAGGCAGTTTGATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATc                                                     <  2:2214853/148‑1 (MQ=255)

                                                                                                  ggTACGAGTCAGTTTTATTGATGTCACCCGTACGATCTATAACCCGTAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGTAAGCATCaa                                                   <  2:2344247/148‑1 (MQ=255)

                                                                                                     aCGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAAtt                                                <  1:805598/148‑1 (MQ=255)

                                                                                                      cccggCAGGTGTACTGAGGACACCCGTACGATCGAGCACCCGCAAGAGAGTGCCTCTGCGCCGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAACAGTGCAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTc                                               <  2:1772196/145‑1 (MQ=255)

                                                                                                       gAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTcc                                              >  1:2086374/1‑148 (MQ=255)

                                                                                                        aGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCt                                             >  1:23029/1‑148 (MQ=255)

                                                                                                         ggCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGTTTTATTGAGGATGGTTTATGACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTg                                            <  2:2229356/148‑1 (MQ=255)

                                                                                                         ggCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTg                                            >  2:1837467/1‑148 (MQ=255)

                                                                                                         ggCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTg                                            >  1:2478589/1‑148 (MQ=255)

                                                                                                         ggCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTg                                            <  2:1768337/148‑1 (MQ=255)

                                                                                                         ggCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTg                                            >  2:2834575/1‑148 (MQ=255)

                                                                                                          gCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGa                                           >  2:507734/1‑148 (MQ=255)

                                                                                                          gCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGa                                           >  2:2740169/1‑148 (MQ=255)

                                                                                                            aGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGaca                                         <  1:2828328/148‑1 (MQ=255)

                                                                                                              ttttATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACAcc                                       >  2:2245535/1‑148 (MQ=255)

                                                                                                              ttttATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACAcc                                       >  2:540950/1‑148 (MQ=255)

                                                                                                               tttATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCa                                      <  1:1588593/148‑1 (MQ=255)

                                                                                                                    tGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTcgc                                 >  1:773857/1‑148 (MQ=255)

                                                                                                                          cACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTgcgc                           >  1:1699474/1‑148 (MQ=255)

                                                                                                                          cACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTgcgc                           >  1:2446140/1‑148 (MQ=255)

                                                                                                                            cccGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGa                         >  2:263263/1‑148 (MQ=255)

                                                                                                                               gTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTTAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTa                      <  2:1912215/148‑1 (MQ=255)

                                                                                                                               gTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTa                      <  2:1603584/148‑1 (MQ=255)

                                                                                                                               gTACGATAGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTTCTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGctacta                      <  2:3157753/148‑1 (MQ=255)

                                                                                                                                tACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGTTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTAt                     <  2:317173/148‑1 (MQ=255)

                                                                                                                                tACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTAt                     <  2:1450961/148‑1 (MQ=255)

                                                                                                                                tACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTAt                     <  1:299841/148‑1 (MQ=255)

                                                                                                                                 aCGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATc                    >  1:2888957/1‑148 (MQ=255)

                                                                                                                                  cgatcgatAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATcg                   <  2:2431896/148‑1 (MQ=255)

                                                                                                                                   gatcgatAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATcgc                  >  1:2152322/1‑148 (MQ=255)

                                                                                                                                    atcgatAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGACAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATcgcg                 >  2:356810/1‑148 (MQ=255)

                                                                                                                                     tcgatAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGACAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCgg                >  2:2104033/1‑148 (MQ=255)

                                                                                                                                     tcgatAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCgg                >  1:2858609/1‑148 (MQ=255)

                                                                                                                                        atAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCaa             <  1:2915661/148‑1 (MQ=255)

                                                                                                                                        atAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCaa             >  1:705737/1‑148 (MQ=255)

                                                                                                                                         tAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCaaa            >  1:2942540/1‑148 (MQ=255)

                                                                                                                                          aaCCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAAc           <  2:2741933/148‑1 (MQ=255)

                                                                                                                                          aaCCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAAc           >  1:2825420/1‑148 (MQ=255)

                                                                                                                                               ggAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATa      <  2:43096/148‑1 (MQ=255)

                                                                                                                                                 aaGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGt    <  2:2852859/148‑1 (MQ=255)

                                                                                                                                                   gagaTTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTAc  <  2:267670/148‑1 (MQ=255)

CTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTAC  >  CP009273/3019056‑3019350

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