BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,814,866	A→C	9.2%	T27P (ACC→CCC)	gmk →	guanylate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,814,866	0	A	C	9.2%	618.9 / 58.0	239	T27P (ACC→CCC)	gmk	guanylate kinase
Reads supporting (aligned to +/- strand): ref base A (115/102); new base C (11/11); total (126/113)
Fisher's exact test for biased strand distribution p-value = 8.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-01

Only 100 of 247 total aligned reads displayed.

CATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAG  >  CP009273/3814719‑3815013

caTGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTaaaaa                                                                                                                                                     >  1:532713/1‑148 (MQ=255)

caTGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTaaaaa                                                                                                                                                     >  1:1213626/1‑148 (MQ=255)

caTGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTaaaaa                                                                                                                                                     >  2:1825065/1‑148 (MQ=255)

   ggCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTAATTAAAAAccc                                                                                                                                                  >  2:1777635/1‑148 (MQ=255)

    gCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCa                                                                                                                                                 <  1:2194339/148‑1 (MQ=255)

     ccGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCaa                                                                                                                                                <  2:2855262/148‑1 (MQ=255)

     ccGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCaa                                                                                                                                                <  1:1634908/148‑1 (MQ=255)

      cgcgCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAAc                                                                                                                                               <  1:109905/148‑1 (MQ=255)

         gcCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGt                                                                                                                                            <  1:1439874/148‑1 (MQ=255)

             tgtATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTa                                                                                                                                        <  2:3202034/148‑1 (MQ=255)

              gtATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTAt                                                                                                                                       >  1:140748/1‑148 (MQ=255)

                ataataAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAACCCCAACCGTTGTATGa                                                                                                                                     >  2:1537560/1‑148 (MQ=255)

                ataataAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAACCCCAACCGTTGTATGa                                                                                                                                     >  2:1537408/1‑148 (MQ=255)

                ataataAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGa                                                                                                                                     >  1:3287764/1‑148 (MQ=255)

                  aataaGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGaca                                                                                                                                   >  2:3444285/1‑148 (MQ=255)

                   ataaGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGacac                                                                                                                                  <  2:1686801/148‑1 (MQ=255)

                    taaGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACAcc                                                                                                                                 >  2:1383283/1‑148 (MQ=255)

                    taaGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACAcc                                                                                                                                 <  1:2749617/148‑1 (MQ=255)

                      aGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAACCCCAACCGTTGTATGACACCCa                                                                                                                               <  2:2822077/148‑1 (MQ=255)

                       gCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAg                                                                                                                              >  1:2697980/1‑148 (MQ=255)

                       gCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAg                                                                                                                              >  2:3367074/1‑148 (MQ=255)

                            tATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTc                                                                                                                         <  1:2162852/148‑1 (MQ=255)

                               gTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGt                                                                                                                      <  1:3474616/148‑1 (MQ=255)

                               gTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGt                                                                                                                      >  1:1932312/1‑148 (MQ=255)

                                   gCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTca                                                                                                                  >  1:1260439/1‑148 (MQ=255)

                                     tttATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTcaca                                                                                                                >  1:3676/1‑148 (MQ=255)

                                      ttATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTcacac                                                                                                               >  1:650048/1‑148 (MQ=255)

                                      ttATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTcacac                                                                                                               >  1:550789/1‑148 (MQ=255)

                                      ttATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTcacac                                                                                                               >  1:1001505/1‑148 (MQ=255)

                                        aTTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCacacac                                                                                                             <  1:1465449/148‑1 (MQ=255)

                                             gCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAACCCCAACCGTTGTATGACACCCAGGTTTCTGTTTCAcacaccacac                                                                                                        >  2:1482388/1‑148 (MQ=255)

                                              cTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACg                                                                                                       >  1:968074/1‑148 (MQ=255)

                                                   cacaTACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAACCCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAAc                                                                                                  <  1:2476118/148‑1 (MQ=255)

                                                   cacaTACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAAc                                                                                                  >  1:190705/1‑148 (MQ=255)

                                                   cacaTACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAAc                                                                                                  >  1:1987265/1‑148 (MQ=255)

                                                   cacaTACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAAc                                                                                                  >  1:2386252/1‑148 (MQ=255)

                                                    acaTACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAACCCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAAcc                                                                                                 <  1:769697/148‑1 (MQ=255)

                                                      aTACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTGGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACcgc                                                                                               >  1:3092355/1‑148 (MQ=255)

                                                       tACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACcgcg                                                                                              >  2:2627741/1‑148 (MQ=255)

                                                         cGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTc                                                                                            <  2:786326/148‑1 (MQ=255)

                                                            aaGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTg                                                                                         >  1:2632092/1‑148 (MQ=255)

                                                            aaGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTg                                                                                         >  1:1573008/1‑148 (MQ=255)

                                                            aaGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTg                                                                                         >  2:31594/1‑148 (MQ=255)

                                                            aaGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTg                                                                                         <  2:1851847/148‑1 (MQ=255)

                                                            aaGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTg                                                                                         >  2:2749700/1‑148 (MQ=255)

                                                                tACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGa                                                                                     >  2:3025816/1‑148 (MQ=255)

                                                                 aCTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAACCCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGaa                                                                                    <  2:1142198/148‑1 (MQ=255)

                                                                 aCTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAACCCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGaa                                                                                    <  2:1142377/148‑1 (MQ=255)

                                                                 aCTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGaa                                                                                    >  1:2534956/1‑148 (MQ=255)

                                                                 aCTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGaa                                                                                    >  2:3220976/1‑148 (MQ=255)

                                                                   tCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGt                                                                                  >  1:3000153/1‑148 (MQ=255)

                                                                    cATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTc                                                                                 >  2:509646/1‑148 (MQ=255)

                                                                    cATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTc                                                                                 >  2:2592592/1‑148 (MQ=255)

                                                                    cATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTc                                                                                 <  1:414760/148‑1 (MQ=255)

                                                                    aaTGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTc                                                                                 >  2:2041076/2‑148 (MQ=255)

                                                                     aTGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTcc                                                                                <  2:968074/148‑1 (MQ=255)

                                                                        gCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAACCCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACg                                                                             >  2:1402756/1‑148 (MQ=255)

                                                                         cTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACgg                                                                            <  2:2002852/148‑1 (MQ=255)

                                                                          tCAAGGCACGCTTTATATTGTTTCTTCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGt                                                                           <  2:3481050/148‑1 (MQ=255)

                                                                           cAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTg                                                                          >  1:1667347/1‑148 (MQ=255)

                                                                            aaGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGa                                                                         <  2:551264/148‑1 (MQ=255)

                                                                              ggCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAACCCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAAc                                                                       <  2:2249129/148‑1 (MQ=255)

                                                                                cACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACat                                                                     <  2:1210286/148‑1 (MQ=255)

                                                                                  cGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACatta                                                                   >  2:2150866/1‑148 (MQ=255)

                                                                                    cTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACattatt                                                                 >  2:2283630/1‑148 (MQ=255)

                                                                                       tataTTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCt                                                              <  1:2028695/148‑1 (MQ=255)

                                                                                           ttGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAACCCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGt                                                          >  2:3262367/1‑148 (MQ=255)

                                                                                           ttGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAACCCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGt                                                          >  1:2308883/1‑148 (MQ=255)

                                                                                             gTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTa                                                        >  1:2858845/1‑148 (MQ=255)

                                                                                             gTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTa                                                        <  2:68074/148‑1 (MQ=255)

                                                                                             gTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTa                                                        <  2:3398906/148‑1 (MQ=255)

                                                                                              tttCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTaa                                                       <  1:3333745/148‑1 (MQ=255)

                                                                                               ttCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAAt                                                      >  2:2026099/1‑148 (MQ=255)

                                                                                                tCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATc                                                     >  1:549459/1‑148 (MQ=255)

                                                                                                    cccccAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCACCCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTCATCatga                                                 >  2:2801745/1‑148 (MQ=255)

                                                                                                      cccAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCatgatg                                               >  1:2870734/1‑148 (MQ=255)

                                                                                                      cccAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCatgatg                                               >  2:2389767/1‑148 (MQ=255)

                                                                                                      cccAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCatgatg                                               >  2:607659/1‑148 (MQ=255)

                                                                                                       ccAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCTCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCAtgatga                                              >  1:1148397/1‑148 (MQ=255)

                                                                                                         aGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAAt                                            <  2:2484158/148‑1 (MQ=255)

                                                                                                               gcgGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAg                                      <  1:1357120/148‑1 (MQ=255)

                                                                                                                cgGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGa                                     <  2:2397981/148‑1 (MQ=255)

                                                                                                                cgGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGa                                     >  2:1233847/1‑148 (MQ=255)

                                                                                                                cgGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTGAAAGa                                     >  1:2252127/1‑148 (MQ=255)

                                                                                                                 gggTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGaa                                    <  2:1312352/148‑1 (MQ=255)

                                                                                                                   gTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAAt                                  >  1:505089/1‑148 (MQ=255)

                                                                                                                        tCCAGCCTGATTCAGGCTTTATTAAAACCCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTa                             <  1:862481/148‑1 (MQ=255)

                                                                                                                        tCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTa                             <  2:3320046/148‑1 (MQ=255)

                                                                                                                        tCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTa                             >  2:3505734/1‑148 (MQ=255)

                                                                                                                          cAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGc                           <  2:95246/148‑1 (MQ=255)

                                                                                                                              cTGATTCAGGATTTATTAAAAACCCAACAGTTGTATGACACCCTGGTTTCTGTTTCACACACAACACGACAACCGCGTCCTGTTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTTAAGAAATGATTAGCatag                       <  2:2930342/148‑1 (MQ=255)

                                                                                                                                gATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGAt                     <  2:2632092/148‑1 (MQ=255)

                                                                                                                                   tCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCg                  <  1:2147550/148‑1 (MQ=255)

                                                                                                                                      ggCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTc               >  1:1776941/1‑148 (MQ=255)

                                                                                                                                      ggCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTc               >  2:776508/1‑148 (MQ=255)

                                                                                                                                        cTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCt             >  1:3184758/1‑148 (MQ=255)

                                                                                                                                          ttattaAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCg           >  2:1741444/1‑148 (MQ=255)

                                                                                                                                             ttaAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAac        <  2:140840/148‑1 (MQ=255)

                                                                                                                                              taAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTCATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAaca       >  2:2739934/1‑148 (MQ=255)

                                                                                                                                                  aaCCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCa   <  1:31415/148‑1 (MQ=255)

CATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAG  >  CP009273/3814719‑3815013

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