Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	11,127	C→G	4.8%	S77T (AGT→ACT)	yaaW ←	UPF0174 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	11,127	0	C	G	4.8%	254.9 / 5.8	83	S77T (AGT→ACT)	yaaW	UPF0174 family protein
Reads supporting (aligned to +/- strand):  ref base C (45/34);  new base G (2/2);  total (47/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-02

CAGGTATTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGC  >  CP009273/10980‑11273                                                                                                                                                    |                                                                                                                                                   cAGGTATTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCTAGCAAAATGACCCGAGACAATTTAATTTGGCCACGCCGTCTGTTCGCGATACt                                                                                                                                                   >  2:619523/1‑149 (MQ=255)    gTATTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATc                                                                                                                                                <  2:466463/149‑1 (MQ=255)     tATTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCg                                                                                                                                               >  2:850883/1‑149 (MQ=255)        tACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGCCGTCGCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGGGTCGGCGCAGTTTGTTGGCGATACTATCGcca                                                                                                                                            >  2:817253/1‑149 (MQ=255)        tACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGACCAGTTTGTTGGCGATCCTATCGcca                                                                                                                                            >  2:1038730/1‑149 (MQ=255)         aCGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATAGTATCGccaa                                                                                                                                           <  1:811638/149‑2 (MQ=255)           gCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGGAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTGGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATAGTATCGccaaca                                                                                                                                         <  2:586389/149‑1 (MQ=255)                aaTTGTTCCAGTAACTGTTGCCCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATg                                                                                                                                    <  1:1029487/149‑1 (MQ=255)                   tGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGAGATCGCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGGGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTg                                                                                                                                 >  2:385832/1‑149 (MQ=255)                      cccAGTAACCGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTaa                                                                                                                              <  1:428732/148‑1 (MQ=255)                          gTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATtct                                                                                                                          >  2:1204475/1‑149 (MQ=255)                          gTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATtct                                                                                                                          <  1:988098/149‑1 (MQ=255)                            aaccgctGCTCAACTTCAAACGTAGACATCTTTTTGTCGGCTTTCAGGTCCAATCGCTTTGAATTCTCGCGCAAAATGCCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTcc                                                                                                                        <  2:345828/143‑1 (MQ=255)                               tGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGACACTTCTGTCTAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTTTTGGCGATACTATCGACACCTAAATGCTTTGATTTTCCGGc                                                                                                                     >  1:722243/1‑149 (MQ=255)                               tGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGc                                                                                                                     >  1:763590/1‑149 (MQ=255)                               tGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGCTTGAAGAAAATACTAGCGCCACCAAAATGCTGTAATTCTCCGGc                                                                                                                     >  1:955877/1‑149 (MQ=255)                                gTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGCAACATCGAGCAAAATGGCCCGATACAATTTACCGCGTCCGCGTAGTTTGTTGGCGATACTATCGCCCCCAAACTGCTGTAATTCTCCGGCa                                                                                                                    >  2:994158/1‑149 (MQ=255)                                     tCAGTTTCAAACGTAGACATCTCTTTGTCGGCTTCCAGCTTCAATCGCTTTGAAACCTCGGGCAAAATGGCCCGATACAATTTACCGTGCGCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAg                                                                                                               <  2:932656/149‑1 (MQ=255)                                     tCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGCTACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAg                                                                                                               >  1:659865/1‑149 (MQ=255)                                     tCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATccg                                                                                                               >  1:47379/1‑147 (MQ=255)                                     tCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAg                                                                                                               <  1:135540/149‑1 (MQ=255)                                     tCAATTTCAAACGTAGACATCTCTTTGTCGGCGTTCAGCTTCAATCGCTTTGAAACATCGAGCCAAATGGCCCGGTACAATTTACCGGGTGCGCGCAGTTTGTTGTCGATACTATCGCCACCAAACTGCAGTGATTCCCCGGCCAACAg                                                                                                               >  2:10582/1‑149 (MQ=255)                                     tCAATTTCAAACGTAGACATCTCTGTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAGAAGGCTGTAATTCTCCGGCAATCAg                                                                                                               >  1:920696/1‑149 (MQ=255)                                       aaTTTCAAACGTAGACATCTCTTTGTCGGTTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCt                                                                                                             <  1:1190537/149‑1 (MQ=255)                                          tttaaaCGTAGACATCTCTTTGTCTGCCTTCAGCTTCAATCGTTTTGAAACAGCGGGCAAAATGGCCCGATACAACTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCAACAAAATGCTGTAATTCTCCGGCAATCAGCTGcc                                                                                                          <  2:402650/146‑1 (MQ=255)                                                 gTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTgct                                                                                                   <  1:798093/149‑2 (MQ=255)                                                 gTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCAGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTACTTTTTCGGCAATTAGTTGGCTGtttgg                                                                                                   >  1:579705/1‑146 (MQ=255)                                                 gTAGACATCTCTTTGGCGGCTCTCAGCTTAAATCGCTGTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTgcg                                                                                                   <  1:536834/149‑1 (MQ=255)                                                 gTAGACAGCTCTTTGTCGGCTTTCAGCTTCAATCTCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTgcg                                                                                                   <  1:431037/149‑1 (MQ=255)                                                  tAGACATTTCTTTGTCCGCTTTAAGCTTCAATCGCTTTGAAACATCGAGCAAAAGGGCCCGATACAATTTACCGTGTCCGTGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTgcgg                                                                                                  <  2:353737/149‑1 (MQ=255)                                                  tAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTgcgg                                                                                                  <  1:850883/149‑1 (MQ=255)                                                        tctcTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCACAATGCTGTAATTCTCCGGCAATCAACTGCCAGTTGCGGCGATGt                                                                                            >  1:17188/1‑149 (MQ=255)                                                        tctcTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGt                                                                                            <  2:1217993/149‑1 (MQ=255)                                                        tctcTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGGTGCGGCGATGt                                                                                            >  1:775480/1‑149 (MQ=255)                                                        tctcTTTGTCGGCTTTCAGATTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGt                                                                                            >  1:899230/1‑149 (MQ=255)                                                        tctcTTTGTCGGCGTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGt                                                                                            >  2:1145131/1‑149 (MQ=255)                                                         ctctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTggggcgccgtt                                                                                           >  2:179605/1‑144 (MQ=255)                                                          tctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTg                                                                                          >  1:501904/1‑149 (MQ=255)                                                          tctTTGTCGGCTTTCAGCTTCAATCACTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTg                                                                                          >  1:639632/1‑149 (MQ=255)                                                            ttcgTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTGGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCt                                                                                        <  1:158703/146‑1 (MQ=255)                                                                   ggtttCAGCCTCAAAGGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGGGTCCGCGCAGTTTGTTGGCGAGACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATg                                                                                 <  2:532802/147‑1 (MQ=255)                                                                       acAGCTTCAATCGCTTTGAAACAGCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCt                                                                             >  1:733940/2‑149 (MQ=255)                                                                           cTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGCTGTTGCTAGGGATGCCCTTCCa                                                                         >  1:247190/1‑149 (MQ=255)                                                                               aaTCGGTTTGAAACATAGAGCAAAATGGCCCGATACAATTTACCGTGGCCGCGCAGTTTGTGGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGa                                                                     <  2:74127/149‑1 (MQ=255)                                                                                aTCGCTTTGAAACATCGAGCAAACTGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCTACCAAAATGCTGTAATTCTCCGGCAATCAGCTGGCAATTGCGGCGATGTTGCTCAGGATGCCCTTCCATCgct                                                                    >  1:681578/1‑147 (MQ=255)                                                                                 tctcTTTGAACCACCGCGCAAAATGTCCCGACACAATTTACCGTGTCCGCGCAGTTTGTCGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGAtt                                                                   <  1:395504/146‑1 (MQ=255)                                                                                   ccTTTGAAACATCGAGCAAATTGGACAGATACAATTTACCGTGTCCGCGCAGTTTGTGGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTa                                                                 <  2:1042963/148‑1 (MQ=255)                                                                                    cTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTaa                                                                >  1:889942/1‑149 (MQ=255)                                                                                     tttGAAACCTCGAGCAGAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATCCTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTaaa                                                               <  2:201826/149‑1 (MQ=255)                                                                                     tttGAAACATCGAGCAAAATGGCCCGCTACAATTTACCGTGTCCGCGCAGTGTGTTGGCGATACAATCGCCACCAAAATGCTGTGATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCACGGGATGCCCTTCAATCGATGTaga                                                               >  2:1062738/1‑147 (MQ=255)                                                                                     tttGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGAATTaaa                                                               >  2:401350/1‑149 (MQ=255)                                                                                     tttGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGGTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTaaa                                                               >  1:1166514/1‑149 (MQ=255)                                                                                     tttGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGATTGTTGGCGATACCATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATAGATTTaaa                                                               >  2:321247/1‑149 (MQ=255)                                                                                           aCATCGAGTAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTc                                                         <  2:696651/149‑1 (MQ=255)                                                                                           aCATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCTGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTc                                                         <  2:793308/149‑1 (MQ=255)                                                                                           aCATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGTTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTc                                                         <  2:1161918/149‑1 (MQ=255)                                                                                             aTCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCAATCGATTTAAACAGTTCGt                                                       >  1:1078370/1‑149 (MQ=255)                                                                                                  gCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTgcgc                                                  <  2:163524/149‑1 (MQ=255)                                                                                                    aaaaTGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTCGAGGCGATGTTGCTCGGGATGACCTTCCATCGATTTAAACAGTTCGTTGCGCAt                                                >  2:683828/1‑149 (MQ=255)                                                                                                              gATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCTTTGCGCATCAGTACGCTg                                      <  1:604315/149‑1 (MQ=255)                                                                                                                   aaTTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGagag                                 >  2:827569/1‑149 (MQ=255)                                                                                                                   aaTTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGagag                                 >  2:827247/1‑149 (MQ=255)                                                                                                                    aTTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAgg                                >  2:157652/1‑149 (MQ=255)                                                                                                                    aTTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGGTCGTTGCGCATCCGTACGCTggggggg                                >  1:623852/1‑144 (MQ=255)                                                                                                                         ccGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCGCCGGAAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTGGTTGCGCATCAGTACGCTGgagaggagagg                           >  2:726438/1‑148 (MQ=255)                                                                                                                           gtgtCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGttt                         <  1:1204475/149‑1 (MQ=255)                                                                                                                            tgtCCGCGCCGTTTATTAGCGATACTATCGCAACCAAAATTCAGTAAGTCTCCAGCAATCAGCTGCCCGGTGCGGCGATGTTGCGCGGGATGCCCTTCCATCGGTTTAGACCGTTCGTTGCGAGTCAGTACGCTGGAGAGGCGAGtgtt                        >  2:604587/1‑149 (MQ=255)                                                                                                                                 gcgcAGTTTGTTGGCGCTACTATCGCCACGAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTACCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCtt                   >  1:784863/1‑149 (MQ=255)                                                                                                                                 gcgcAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCtt                   >  1:304663/1‑149 (MQ=255)                                                                                                                                  agcATTTTGGTGGCGATAGTATCGCCAACAAACTGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTGATACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCttt                  >  1:1158309/2‑149 (MQ=38)                                                                                                                                   gcAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCtttt                 <  1:421505/149‑1 (MQ=255)                                                                                                                                   gcAGTTTGTTGGCGATACCATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCtttt                 <  2:1154753/149‑1 (MQ=255)                                                                                                                                      gTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCa              <  2:659865/149‑1 (MQ=255)                                                                                                                                      gTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCa              >  1:685249/1‑149 (MQ=255)                                                                                                                                         tgttgGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCTAGTTTTGCCTTTTTCATTa           <  2:247190/149‑1 (MQ=255)                                                                                                                                          ggtgGCGATAGTATCGCCAACAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGCTTTGCCTTTTTCCTTAt          >  2:753073/3‑149 (MQ=255)                                                                                                                                           ttgGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGCGTTTTGCCTTTTTCATTATg         >  2:633537/1‑149 (MQ=255)                                                                                                                                            tgGCGATACTACCGCCACCAAAATGCTGTAATTCGCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATgg        <  2:85188/149‑1 (MQ=255)                                                                                                                                             ggCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGGTGCCAGTTGCGGAGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGGGTTTTGCCTTTTTCATTATggg       >  1:859335/1‑149 (MQ=255)                                                                                                                                              gCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGATCGGGATGCCCTTCCATCGATTTCAACAGTTCGTTGCGCATCAGTGCGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGt      >  2:416741/1‑149 (MQ=255)                                                                                                                                                 aTACTATCGCCACCAAAATGCTGTAATTATCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGTGGCGCGGTTTGCCTTTTTAATTATGGGTggg   >  1:300413/1‑147 (MQ=255)                                                                                                                                                  tACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGc  <  2:775480/149‑1 (MQ=255)                                                                                                                                                  tACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTACTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGc  <  2:899230/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   CAGGTATTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGC  >  CP009273/10980‑11273

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑