BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,931,717	C→A	92.9%	N307K (AAC→AAA)	fucK →	L‑fuculokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,931,717	0	C	A	92.9%	1060.8 / 66.1	338	N307K (AAC→AAA)	fucK	L‑fuculokinase
Reads supporting (aligned to +/- strand): ref base C (10/14); new base A (164/150); total (174/164)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 341 total aligned reads displayed.

TGGCGCTGGTGCTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTC  >  CP009273/2931570‑2931866

 ggCGCTGGTGCTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCa                                                                                                                                                    >  2:1813395/1‑150 (MQ=255)

   cgcTGGTGCTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAgg                                                                                                                                                  >  1:3717299/1‑150 (MQ=255)

     cTGGTGCTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTa                                                                                                                                                <  1:2511972/150‑1 (MQ=255)

         tGCTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCa                                                                                                                                            <  2:1963172/150‑1 (MQ=255)

         tGCTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCa                                                                                                                                            <  2:410449/150‑1 (MQ=255)

         tGCTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCa                                                                                                                                            >  1:199473/1‑150 (MQ=255)

           cTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAAt                                                                                                                                          >  1:4568673/1‑150 (MQ=255)

               aCAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAAtggct                                                                                                                                      >  1:2616414/1‑150 (MQ=255)

               aCAAAATGAACCCGTGCTCTCTGCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAAtggct                                                                                                                                      <  2:3553451/150‑1 (MQ=255)

                  aaaTGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAAtggctggc                                                                                                                                   >  1:1399774/1‑150 (MQ=255)

                  aaaTGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAAtggctggc                                                                                                                                   <  2:4353970/150‑1 (MQ=255)

                     tGAACCCGTGCTCTCTGCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATc                                                                                                                                <  2:1033862/150‑1 (MQ=255)

                      gAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATcc                                                                                                                               >  2:4477390/1‑150 (MQ=255)

                      gAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATcc                                                                                                                               <  1:3077828/150‑1 (MQ=255)

                       aaCCCGTGCTCTCTTCCGGGACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCg                                                                                                                              <  1:3349698/150‑1 (MQ=255)

                         cccGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGc                                                                                                                            >  1:1237106/1‑150 (MQ=255)

                            gTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTg                                                                                                                         >  1:4310122/1‑150 (MQ=255)

                             tGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGc                                                                                                                        >  1:187565/1‑150 (MQ=255)

                              gCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCt                                                                                                                       >  2:725083/1‑150 (MQ=255)

                               ctctctTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTg                                                                                                                      <  2:1163335/150‑1 (MQ=255)

                                 ctctTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGa                                                                                                                    >  1:3923320/1‑150 (MQ=255)

                                   ctTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAAt                                                                                                                  <  2:647058/150‑1 (MQ=255)

                                       cGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGt                                                                                                              <  2:2777650/150‑1 (MQ=255)

                                         gTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTga                                                                                                            >  1:2400603/1‑150 (MQ=255)

                                          tACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTgag                                                                                                           <  2:1481679/150‑1 (MQ=255)

                                                ggAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACt                                                                                                     >  2:3005019/1‑150 (MQ=255)

                                                      tttAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTg                                                                                               <  1:316393/150‑1 (MQ=255)

                                                      tttAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTg                                                                                               >  1:438340/1‑150 (MQ=255)

                                                       ttAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTgg                                                                                              >  1:4670368/1‑150 (MQ=255)

                                                       ttAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTgg                                                                                              >  2:3495989/1‑150 (MQ=255)

                                                         aaTGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGAc                                                                                            >  1:3558934/1‑150 (MQ=255)

                                                          aTGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTgg                                                                                                                     <  2:2166069/124‑1 (MQ=255)

                                                            ggTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGc                                                                                         >  2:538397/1‑150 (MQ=255)

                                                            ggTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGCTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGc                                                                                         >  2:3988336/1‑150 (MQ=255)

                                                              ttCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTg                                                                                       <  2:1592726/150‑1 (MQ=255)

                                                                cGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGaa                                                                                     <  1:3636320/150‑1 (MQ=255)

                                                                cGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGaa                                                                                     <  1:523653/150‑1 (MQ=255)

                                                                   aGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAaca                                                                                  >  2:3350558/1‑150 (MQ=255)

                                                                        ccAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTg                                                                             <  1:881139/150‑1 (MQ=255)

                                                                        ccAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTg                                                                             <  2:4770286/150‑1 (MQ=255)

                                                                        ccAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTg                                                                             <  2:3641285/150‑1 (MQ=255)

                                                                         cAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACg                                                                                           <  1:2218678/135‑1 (MQ=255)

                                                                              tGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAAt                                                                       <  1:952860/150‑1 (MQ=255)

                                                                              tGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAAt                                                                       >  2:2427667/1‑150 (MQ=255)

                                                                                aTACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGt                                                                     <  2:2916483/150‑1 (MQ=255)

                                                                                 tACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGtt                                                                    >  1:1851561/1‑150 (MQ=255)

                                                                                 tACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGtt                                                                    >  1:4252470/1‑150 (MQ=255)

                                                                                  aCTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGttg                                                                   <  2:1743851/150‑1 (MQ=255)

                                                                                  aCTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGttg                                                                   >  2:4183255/1‑150 (MQ=255)

                                                                                  aCTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGttg                                                                   <  1:4477390/150‑1 (MQ=255)

                                                                                     tCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGttgatt                                                                >  1:2431971/1‑150 (MQ=255)

                                                                                      cGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGttgattg                                                               >  2:3762994/1‑150 (MQ=255)

                                                                                       gCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGttgattga                                                              <  2:2849655/150‑1 (MQ=255)

                                                                                       gCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGttgattga                                                              >  1:3979570/1‑150 (MQ=255)

                                                                                         tGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTgaag                                                            <  2:3855518/150‑1 (MQ=255)

                                                                                           ttAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTgaagaa                                                          <  1:1003215/150‑1 (MQ=255)

                                                                                             aaGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGc                                                        >  1:1600653/1‑150 (MQ=255)

                                                                                              aGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCt                                                       <  1:771275/150‑1 (MQ=255)

                                                                                              aGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCt                                                       <  2:2412352/150‑1 (MQ=255)

                                                                                               gTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTc                                                      <  2:2339601/150‑1 (MQ=255)

                                                                                               gTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTc                                                      <  2:187565/150‑1 (MQ=255)

                                                                                               gTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTc                                                      <  2:2006578/150‑1 (MQ=255)

                                                                                                tCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCg                                                     >  1:678236/1‑150 (MQ=255)

                                                                                                  aGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTc                                                   <  1:2950498/150‑1 (MQ=255)

                                                                                                      cGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGAt                                               >  2:3181782/1‑150 (MQ=255)

                                                                                                       gCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATc                                              >  1:4539626/1‑150 (MQ=255)

                                                                                                          ggTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGACTGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATcgcg                                           >  2:1233290/1‑150 (MQ=255)

                                                                                                             tCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCt                                        <  1:2736776/150‑1 (MQ=255)

                                                                                                               cACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTgg                                      >  1:410110/1‑150 (MQ=255)

                                                                                                                aCCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGgc                                     <  1:1778320/150‑1 (MQ=255)

                                                                                                                   tGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGcgcg                                  >  2:3138402/1‑150 (MQ=255)

                                                                                                                     cGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGa                                <  2:4258884/150‑1 (MQ=255)

                                                                                                                     cGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGa                                <  2:3022473/150‑1 (MQ=255)

                                                                                                                      gAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGAt                               <  2:2260414/150‑1 (MQ=255)

                                                                                                                       aaCTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATg                              >  1:4288149/1‑150 (MQ=255)

                                                                                                                          tGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCg                           >  1:3321582/1‑150 (MQ=255)

                                                                                                                           ggATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGt                          <  1:1055720/150‑1 (MQ=255)

                                                                                                                           ggATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGt                          <  1:4653391/150‑1 (MQ=255)

                                                                                                                            gATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTa                         <  1:2919964/150‑1 (MQ=255)

                                                                                                                            gATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTa                         <  2:4219319/150‑1 (MQ=255)

                                                                                                                                 ccAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAAt                    <  2:26356/150‑1 (MQ=255)

                                                                                                                                  caggcaggGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATg                   <  1:1430771/150‑1 (MQ=255)

                                                                                                                                   aggcaggGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGc                  <  2:1705433/150‑1 (MQ=255)

                                                                                                                                   aggcaggGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGc                  <  1:863925/150‑1 (MQ=255)

                                                                                                                                    ggcaggGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCa                 >  1:3859433/1‑150 (MQ=255)

                                                                                                                                       aggGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAgtg              >  2:248013/1‑150 (MQ=255)

                                                                                                                                       aggGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAgtg              >  1:477392/1‑150 (MQ=255)

                                                                                                                                         ggTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGtgtg            <  2:884241/150‑1 (MQ=255)

                                                                                                                                         ggTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGtgtg            <  2:3717299/150‑1 (MQ=255)

                                                                                                                                          gTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGa           >  1:743/1‑150 (MQ=255)

                                                                                                                                          gTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGa           <  2:1753491/150‑1 (MQ=255)

                                                                                                                                          gTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGa           >  2:592092/1‑150 (MQ=255)

                                                                                                                                            tGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGAtt         >  1:4440459/1‑150 (MQ=255)

                                                                                                                                              tataAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTa       >  2:4312741/1‑150 (MQ=255)

                                                                                                                                               ataACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTAt      <  1:3734279/150‑1 (MQ=255)

                                                                                                                                               ataAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTAt      >  2:1630848/1‑150 (MQ=255)

                                                                                                                                               ataAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTAt      >  1:1828211/1‑150 (MQ=255)

                                                                                                                                               ataAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTAt      >  1:1514996/1‑150 (MQ=255)

                                                                                                                                                taAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTAtt     <  2:2701453/150‑1 (MQ=255)

                                                                                                                                                taAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTAtt     <  2:309972/150‑1 (MQ=255)

TGGCGCTGGTGCTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTC  >  CP009273/2931570‑2931866

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