BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,970,026	C→A	8.7%	A130D (GCC→GAC)	wzxE →	O‑antigen translocase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,970,026	0	C	A	8.7%	904.0 / 38.2	252	A130D (GCC→GAC)	wzxE	O‑antigen translocase
Reads supporting (aligned to +/- strand): ref base C (122/108); new base A (11/11); total (133/119)
Fisher's exact test for biased strand distribution p-value = 8.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 261 total aligned reads displayed.

GCACTTCATCAGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGTACTACGTCAGTTACCGTTTGGGCGGTTATGAAGGGGCGTTGCTGGGTCT  >  CP009273/3969879‑3970173

gCACTTCATCAGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGc                                                                                                                                                     >  1:2546951/1‑148 (MQ=255)

 cACTTCATCAGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGcc                                                                                                                                                    >  2:58117/1‑148 (MQ=255)

   cTTCATCAGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTg                                                                                                                                                  <  2:888763/148‑1 (MQ=255)

     tcatcaGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTggg                                                                                                                                                >  2:826221/1‑148 (MQ=255)

     tcatcaGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTggg                                                                                                                                                >  2:665126/1‑148 (MQ=255)

      catcaGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTgggg                                                                                                                                               >  1:1073529/1‑148 (MQ=255)

         caGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCaa                                                                                                                                            >  1:1933944/1‑148 (MQ=255)

         caGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCaa                                                                                                                                            >  1:1933969/1‑148 (MQ=255)

          aGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAAc                                                                                                                                           <  2:2437280/148‑1 (MQ=255)

             gATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTg                                                                                                                                        <  1:2353349/148‑1 (MQ=255)

                ggTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTa                                                                                                                                     <  1:2889531/148‑1 (MQ=255)

                ggTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTa                                                                                                                                     <  1:1697267/148‑1 (MQ=255)

                 gTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGACTGGGGCAACCTGTTAc                                                                                                                                    >  1:3182076/1‑148 (MQ=255)

                  tACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACt                                                                                                                                   >  1:1557540/1‑148 (MQ=255)

                   aCTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTg                                                                                                                                  >  2:374125/1‑148 (MQ=255)

                      tGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGgcg                                                                                                                               <  2:1376890/148‑1 (MQ=255)

                      tGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGgcg                                                                                                                               <  2:1230800/148‑1 (MQ=255)

                            ctctACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCtgatg                                                                                                                         <  1:1055941/148‑1 (MQ=255)

                            ctctACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCtgatg                                                                                                                         <  2:3203071/148‑1 (MQ=255)

                              ctACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGaa                                                                                                                       >  2:2615572/1‑148 (MQ=255)

                               tACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGaaa                                                                                                                      <  1:57247/148‑1 (MQ=255)

                               tACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGaaa                                                                                                                      <  1:398060/148‑1 (MQ=255)

                               tACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGaaa                                                                                                                      <  1:2485207/148‑1 (MQ=255)

                                    tGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGACTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCtt                                                                                                                 <  1:2012893/148‑1 (MQ=255)

                                     gATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCttt                                                                                                                <  2:970639/148‑1 (MQ=255)

                                       tGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTcg                                                                                                              >  2:2327314/1‑148 (MQ=255)

                                           gcTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGAt                                                                                                          >  1:607335/1‑148 (MQ=255)

                                             tGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGc                                                                                                        >  1:3224218/1‑148 (MQ=255)

                                               gTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCg                                                                                                      >  2:2951428/1‑148 (MQ=255)

                                                tttttGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGTATCGACTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGc                                                                                                     >  1:121145/1‑148 (MQ=255)

                                                 ttttGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCa                                                                                                    >  2:10548/1‑148 (MQ=255)

                                                     gtgCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTa                                                                                                <  1:908816/148‑1 (MQ=255)

                                                      tgCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTaa                                                                                               <  2:1506675/148‑1 (MQ=255)

                                                         tGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGc                                                                                            >  1:2846431/1‑148 (MQ=255)

                                                         tGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGc                                                                                            <  2:1596892/148‑1 (MQ=255)

                                                            cAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGGGTTTAGTGGCGCTGGTGCAAATGGTGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGtt                                                                                         <  1:2004387/148‑1 (MQ=255)

                                                            cAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGtt                                                                                         <  2:2011908/148‑1 (MQ=255)

                                                              gCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTAt                                                                                       >  2:1090840/1‑148 (MQ=255)

                                                               cTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATc                                                                                      >  1:2977844/1‑148 (MQ=255)

                                                                tGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCg                                                                                     <  2:1313209/148‑1 (MQ=255)

                                                                tGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCg                                                                                     <  2:3082222/148‑1 (MQ=255)

                                                                 gcgcCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGc                                                                                    <  2:1795278/148‑1 (MQ=255)

                                                                    ccAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGa                                                                                 >  1:434277/1‑148 (MQ=255)

                                                                     cAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGAt                                                                                >  1:3286139/1‑148 (MQ=255)

                                                                     cAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGAt                                                                                <  1:2966064/148‑1 (MQ=255)

                                                                         cAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTc                                                                            >  1:2680610/1‑148 (MQ=255)

                                                                            ccAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGc                                                                         <  1:910160/148‑1 (MQ=255)

                                                                            ccAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGc                                                                         <  1:2736691/148‑1 (MQ=255)

                                                                              aGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAg                                                                       >  1:792809/1‑148 (MQ=255)

                                                                              aGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAg                                                                       >  1:1488004/1‑148 (MQ=255)

                                                                               gggATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAAGGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGc                                                                      <  2:1173451/148‑1 (MQ=255)

                                                                                ggATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCt                                                                     <  2:2071603/148‑1 (MQ=255)

                                                                                  aTTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTg                                                                   >  2:439502/1‑148 (MQ=255)

                                                                                   ttGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGa                                                                  >  2:823746/1‑148 (MQ=255)

                                                                                    tGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGAt                                                                 >  1:3287259/1‑148 (MQ=255)

                                                                                     gTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGAtt                                                                <  2:2661167/148‑1 (MQ=255)

                                                                                      tttGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTg                                                               >  1:1419196/1‑148 (MQ=255)

                                                                                       ttGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGACTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTgg                                                              >  2:1104335/1‑148 (MQ=255)

                                                                                           tAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGtt                                                          >  1:949407/1‑148 (MQ=255)

                                                                                           tAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGtt                                                          >  2:2070782/1‑148 (MQ=255)

                                                                                            aaTACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTtc                                                         <  2:1920267/148‑1 (MQ=255)

                                                                                            aaTACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTtc                                                         <  2:3239583/148‑1 (MQ=255)

                                                                                              tACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTtctc                                                       >  2:489696/1‑148 (MQ=255)

                                                                                               aCCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTcg                                                      >  2:3272562/1‑148 (MQ=255)

                                                                                                 cGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTcgcg                                                    >  1:2062517/1‑148 (MQ=255)

                                                                                                   aCTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGta                                                  >  1:2235857/1‑148 (MQ=255)

                                                                                                     tATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGtact                                                >  2:2480031/1‑148 (MQ=255)

                                                                                                     tATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGtact                                                >  2:2481083/1‑148 (MQ=255)

                                                                                                        cAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGTACTACg                                             >  2:275667/1‑148 (MQ=255)

                                                                                                        cAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGACTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGTACTACg                                             >  2:189447/1‑148 (MQ=255)

                                                                                                         aGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGTACTACGt                                            >  1:36404/1‑148 (MQ=255)

                                                                                                         aGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGTACTACGt                                            >  1:3044484/1‑148 (MQ=255)

                                                                                                           gggCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGTACTACGTCa                                          <  2:2796691/148‑1 (MQ=255)

                                                                                                              cTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGTACTACGTCAGtt                                       <  2:3127270/148‑1 (MQ=255)

                                                                                                                ggTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGTACTACGTCAGTTAc                                     >  1:2336286/1‑148 (MQ=255)

                                                                                                                   gCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGTACTACGTCAGTTACCGt                                  >  1:3209732/1‑148 (MQ=255)

                                                                                                                     gTTTAGTGGCGCTGGTGCAAATGGGTATCGACTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGTACTACGTCAGTTACCGttt                                <  2:121145/148‑1 (MQ=255)

                                                                                                                      tttAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGTACTACGTCAGTTACCGTTTg                               >  1:1593037/1‑148 (MQ=255)

                                                                                                                       ttAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGTACTACGTCAGTTACCGTTTgg                              >  1:1255547/1‑148 (MQ=255)

                                                                                                                       ttAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGTACTACGTCAGTTACCGTTTgg                              >  2:2851563/1‑148 (MQ=255)

                                                                                                                       ttAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGTACTACGTCAGTTACCGTTTgg                              >  1:973984/1‑148 (MQ=255)

                                                                                                                           tGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGTACTACGTCAGTTACCGTTTGGGCgg                          <  1:628122/148‑1 (MQ=255)

                                                                                                                            ggCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGTACTACGTCAGTTACCGTTTGGGCGGt                         <  2:1441417/148‑1 (MQ=255)

                                                                                                                              cgcTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGTACTACGTCAGTTACCGTTTGGGCGGTTa                       >  1:679581/1‑148 (MQ=255)

                                                                                                                              cgcTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGTACTACGTCAGTTACCGTTTGGGCGGTTa                       <  1:3023970/148‑1 (MQ=255)

                                                                                                                                 tGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGTACTACGTCAGTTACCGTTTGGGCGGTTATGa                    >  1:1294061/1‑148 (MQ=255)

                                                                                                                                     gCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGTACTACGTCAGTTACCGTTTGGGCGGTTATGAAggg                >  2:169138/1‑148 (MQ=255)

                                                                                                                                     gCAAATGGGGATCGCCTGGGGCAACCTGATACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGTACTACGTCAGTTACCGTTTGGGCGGTTATGAAggg                >  2:1839668/1‑148 (MQ=255)

                                                                                                                                       aaaTGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGTACTACGTCAGTTACCGTTTGGGCGGTTATGAAGGGGc              <  1:495897/148‑1 (MQ=255)

                                                                                                                                        aaTGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGTACTACGTCAGTTACCGTTTGGGCGGTTATGAAGGGGCg             <  1:1328150/148‑1 (MQ=255)

                                                                                                                                        aaTGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGTACTACGTCAGTTACCGTTTGGGCGGTTATGAAGGGGCg             <  1:2157172/148‑1 (MQ=255)

                                                                                                                                         aTGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGTACTACGTCAGTTACCGTTTGGGCGGTTATGAAGGGGCGt            <  2:3091620/148‑1 (MQ=255)

                                                                                                                                         aTGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTGTCGCGATGCCGCAGGTAATGGGTTATCGGTGATTGTCGGCAGGGTGATTGGCGTTCTCGCGTACTAGGTCAGTGACGGGTTGGGGGGTTTTGTAggggggg            >  2:2169645/1‑145 (MQ=255)

                                                                                                                                          tGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGTACTACGTCAGTTACCGTTTGGGCGGTTATGAAGGGGCGtt           >  1:2819874/1‑148 (MQ=255)

                                                                                                                                             ggATCGACTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGTACTACGTCAGTTACCGTTTGGGCGGTTATGAAGGGGCGTTGCt        <  1:1104335/148‑1 (MQ=255)

                                                                                                                                                tCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGTACTACGTCAGTTACCGTTTGGGCGGTTATGAAGGGGCGTTGCTggg     <  1:881216/148‑1 (MQ=255)

                                                                                                                                                tCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGTACTACGTCAGTTACCGTTTGGGCGGTTATGAAGGGGCGTTGCTggg     >  2:917407/1‑148 (MQ=255)

                                                                                                                                                tCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGTACTACGTCAGTTACCGTTTGGGCGGTTATGAAGGGGCGTTGCTggg     >  2:917448/1‑148 (MQ=255)

                                                                                                                                                  gCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGTACTACGTCAGTTACCGTTTGGGCGGTTATGAAGGGGCGTTGCTGGGTc   <  2:729472/148‑1 (MQ=255)

                                                                                                                                                   ccTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGTACTACGTCAGTTACCGTTTGGGCGGTTATGAAGGGGCGTTGCTGGGTCt  <  1:824566/148‑1 (MQ=255)

GCACTTCATCAGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAATGCGTTATCGCTGATTGTCGGCAGCTTGATTGGCGTTCTCGCGTACTACGTCAGTTACCGTTTGGGCGGTTATGAAGGGGCGTTGCTGGGTCT  >  CP009273/3969879‑3970173

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