BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,814,866	A→C	7.0%	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	7.0%	711.0 / 44.9	256	T27P (ACC→CCC)	gmk	guanylate kinase
Reads supporting (aligned to +/- strand): ref base A (130/108); new base C (9/9); total (139/117)
Fisher's exact test for biased strand distribution p-value = 8.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

Only 100 of 271 total aligned reads displayed.

CATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAG  >  CP009273/3814719‑3815013

caTGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTaaaaa                                                                                                                                                     <  2:1550281/148‑1 (MQ=255)

   ggCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAAccc                                                                                                                                                  <  2:2268612/148‑1 (MQ=255)

      cgcgCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAAc                                                                                                                                               >  1:1677172/1‑148 (MQ=255)

       gcgcCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAACCCCAAcc                                                                                                                                              >  1:206359/1‑148 (MQ=255)

        cgcCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCg                                                                                                                                             >  2:2162156/1‑148 (MQ=255)

           cGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTg                                                                                                                                          >  1:2753791/1‑148 (MQ=255)

             tgtATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTa                                                                                                                                        >  2:2976094/1‑148 (MQ=255)

               tataATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATg                                                                                                                                      <  1:367349/148‑1 (MQ=255)

               tataATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATg                                                                                                                                      >  1:2272051/1‑148 (MQ=255)

               tataATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATg                                                                                                                                      >  1:3395044/1‑148 (MQ=255)

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

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

                         tCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGt                                                                                                                            >  1:398534/1‑148 (MQ=255)

                         tCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGt                                                                                                                            <  2:476910/148‑1 (MQ=255)

                         tCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGt                                                                                                                            >  1:1314556/1‑148 (MQ=255)

                          cGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGtt                                                                                                                           <  2:2227405/148‑1 (MQ=255)

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

                             aTGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCt                                                                                                                        <  1:838222/148‑1 (MQ=255)

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

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

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

                                 aGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGttt                                                                                                                    <  1:908651/148‑1 (MQ=255)

                                 aGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGttt                                                                                                                    <  2:2773950/148‑1 (MQ=255)

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

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

                                         tttCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCAcacacc                                                                                                            <  1:1867525/148‑1 (MQ=255)

                                         tttCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCAcacacc                                                                                                            >  1:3355452/1‑148 (MQ=255)

                                              cTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTGTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCAGAGTCAGGCTTTAGTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACg                                                                                                       >  2:1950373/1‑148 (MQ=255)

                                               tAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGc                                                                                                      >  1:2665645/1‑148 (MQ=255)

                                                 aTCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCa                                                                                                    <  2:1858556/148‑1 (MQ=255)

                                                 aTCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCa                                                                                                    >  1:1910571/1‑148 (MQ=255)

                                                   cacaTACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAAc                                                                                                  <  1:1650329/148‑1 (MQ=255)

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

                                                          gAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTcc                                                                                           <  2:3240553/148‑1 (MQ=255)

                                                           aaaGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCt                                                                                          >  1:826918/1‑148 (MQ=255)

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

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

                                                             aGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTgg                                                                                        >  1:1058847/1‑148 (MQ=255)

                                                              gATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGt                                                                                       >  2:1126360/1‑148 (MQ=255)

                                                              gATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGt                                                                                       >  1:1339142/1‑148 (MQ=255)

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

                                                                   tCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGt                                                                                  <  1:731415/148‑1 (MQ=255)

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

                                                                   tCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGt                                                                                  <  1:3089061/148‑1 (MQ=255)

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

                                                                        gCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACg                                                                             >  2:2760124/1‑148 (MQ=255)

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

                                                                           cAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTg                                                                          <  1:2274529/148‑1 (MQ=255)

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

                                                                              ggCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAAc                                                                       >  2:2837401/1‑148 (MQ=255)

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

                                                                                cACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACat                                                                     >  1:702663/1‑148 (MQ=255)

                                                                                cACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACat                                                                     >  1:2620880/1‑148 (MQ=255)

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

                                                                                    cTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACattatt                                                                 <  2:3104202/148‑1 (MQ=255)

                                                                                      ttATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTc                                                               <  1:1336003/148‑1 (MQ=255)

                                                                                      ttATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTc                                                               <  1:2308777/148‑1 (MQ=255)

                                                                                       tataTTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCt                                                              >  1:1209195/1‑148 (MQ=255)

                                                                                        atatTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCtt                                                             <  2:2424406/148‑1 (MQ=255)

                                                                                        atatTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCtt                                                             >  2:2119791/1‑148 (MQ=255)

                                                                                          atTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTg                                                           <  1:1752843/148‑1 (MQ=255)

                                                                                           ttGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGt                                                          >  1:761292/1‑148 (MQ=255)

                                                                                            tGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGtt                                                         >  2:1294631/1‑148 (MQ=255)

                                                                                             gTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAACCCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTa                                                        <  2:366428/148‑1 (MQ=255)

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

                                                                                              tttCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAACCCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTaa                                                       <  1:328014/148‑1 (MQ=255)

                                                                                              tttCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTaa                                                       >  1:1172859/1‑148 (MQ=255)

                                                                                                 cTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCa                                                    >  2:790438/1‑148 (MQ=255)

                                                                                                   gCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCatg                                                  <  2:1314556/148‑1 (MQ=255)

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

                                                                                                      cccAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCatgatg                                               <  2:3253087/148‑1 (MQ=255)

                                                                                                        cAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGaa                                             >  1:1412220/1‑148 (MQ=255)

                                                                                                         aGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAAt                                            >  1:750349/1‑148 (MQ=255)

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

                                                                                                           tGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGCTGAAttt                                          >  2:1228154/1‑148 (MQ=255)

                                                                                                           tGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAAttt                                          >  1:2362307/1‑148 (MQ=255)

                                                                                                            ggCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTa                                         <  2:711048/148‑1 (MQ=255)

                                                                                                            ggCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTa                                         <  2:3372205/148‑1 (MQ=255)

                                                                                                                  ggTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGaaa                                   >  1:3011802/1‑148 (MQ=255)

                                                                                                                    tAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATg                                 >  1:2576798/1‑148 (MQ=255)

                                                                                                                    tAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATg                                 <  2:1737884/148‑1 (MQ=255)

                                                                                                                       aTCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGAtt                              >  1:661278/1‑148 (MQ=255)

                                                                                                                       aTCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGAtt                              >  1:661472/1‑148 (MQ=255)

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

                                                                                                                           aGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCa                          >  1:1512400/1‑148 (MQ=255)

                                                                                                                           aGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCa                          >  1:3124035/1‑148 (MQ=255)

                                                                                                                               tGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAgaga                      <  1:462757/148‑1 (MQ=255)

                                                                                                                                    cAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGt                 <  1:3223605/148‑1 (MQ=255)

                                                                                                                                    cAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTCAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGt                 <  1:2162156/148‑1 (MQ=255)

                                                                                                                                     aGGCTTTATTAAAACCCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGtt                >  2:1643793/1‑148 (MQ=255)

                                                                                                                                     aGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGtt                >  2:578564/1‑148 (MQ=255)

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

                                                                                                                                        cTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCt             <  2:3091825/148‑1 (MQ=255)

                                                                                                                                           tattaAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGa          <  2:398534/148‑1 (MQ=255)

                                                                                                                                            attaAAACCCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGaa         <  1:939188/148‑1 (MQ=255)

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

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

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

                                                                                                                                                 aaaCCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGc    >  2:118437/1‑148 (MQ=255)

                                                                                                                                                   aCCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAg  >  2:1055257/1‑148 (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: `‑`