Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,250,324	G→A	19.7%	A26T (GCA→ACA)	yjbL →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,250,324	0	G	A	19.7%	188.5 / 40.3	81	A26T (GCA→ACA)	yjbL	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (26/39);  new base A (6/10);  total (32/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.78e-01

TTTGATAGTGACCAGATGGCAAAAGAATTATAAGGAGGACACGGTACGACAAAAAAGGTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCGGGATATAATTATTGCAAAAATAGATTTGTTTATTCAACGAATCACGCA  >  CP009273/4250181‑4250466                                                                                                                                                |                                                                                                                                               tttGATAGTGACCAGATGGCAAAAGAATTATAAGGAGGACACGGTACGACAAAAAAGGTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCa                                                                                                                                          <  2:3177994/150‑1 (MQ=255)  ttGATAGTGACCAGATGGCAAAAGAATTATAAGGAGGACACGGTACGACTAAAAAGGTAGATGAATCATGTTGAAGATAAGACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCaa                                                                                                                                         <  2:2203661/150‑1 (MQ=255)      tAGTGACCAGATTTCAAAAGAATTATAAGGATGACACTGTACGACAAAAAAGGTATATGAATCATGTTTAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGAt                                                                                                                                     <  2:2552058/150‑1 (MQ=255)                aTGGCAAAAGAATTATAAGGAGGACACGGTACGACAAAAAAGGTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGata                                                                                                                           <  1:1457977/150‑1 (MQ=255)                 tgtcAAAAGAATTATAAGGAGGACACGGTACGACAAAAAAGGTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGataa                                                                                                                          <  2:456971/147‑1 (MQ=255)                 tGGCAAAAGAATTATAAGGAGGACACGGTACGACAAAAAAGGTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGataa                                                                                                                          <  2:2835656/150‑1 (MQ=255)                 tGGCAAAAGAATTATAAGGAGGACACGGTACGACAAAAAAGGTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGataa                                                                                                                          >  1:3031916/1‑150 (MQ=255)                 tGGCAAAAGAATTATAAGGAGGACACGGTACGACAAAAAAGGTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGataa                                                                                                                          >  1:3030798/1‑150 (MQ=255)                    cAAAAGAATTATAAGGAGGACACGGTACGACAAAAAAGGTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATACAATCAAAGATAAATTAGATAATag                                                                                                                       >  1:3018042/1‑150 (MQ=255)                    cAAAAGAATTATAAGGAGGACACGGTACGACAAAAAAGGTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATACAATCAAAGATAAATTAGATAATag                                                                                                                       >  1:3018034/1‑150 (MQ=255)                     aaaaGAATTATAAGGAGGACACGGTACGACAAAAAAGGTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATaga                                                                                                                      <  1:4716887/150‑1 (MQ=255)                     aaaaGAATTATAAGGAGGACACGGTACGACAAAAAAGGTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATaga                                                                                                                      <  1:1596224/150‑1 (MQ=255)                     aaaaGAATTATAAGGAGGACACGGTACGACAAAAAAGGTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATACAATCAAAGATAAATTAGATAATaga                                                                                                                      <  2:443857/150‑1 (MQ=255)                     aaaaGAATTATAAGGAGGACACGGTACGACAAAAAAGGTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATACAATCAAAGATAAATTAGATAATaga                                                                                                                      <  2:2045142/150‑1 (MQ=255)                     aaaaGAATTATAAGGAGGACACGGTACGACAAAAAAGGTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATACAATCAAAGATAAATTAGATAATaga                                                                                                                      <  2:2044329/150‑1 (MQ=255)                        aGAATTATAAGGAGGACACGGTACGACAAAAAAGGTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATACAATCAAAGATAAATTAGATAATAGAGGt                                                                                                                   >  2:1341158/1‑150 (MQ=255)                          aaTTATAAGGAGGACACGGTACGACAAAAAAGGTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTtc                                                                                                                 <  2:2427507/150‑1 (MQ=255)                          aaTTATAAGGAGGACACGGTACGACAAAAAAGGTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTtc                                                                                                                 <  2:2427495/150‑1 (MQ=255)                             tataAGGAGGACACGGTACGACAAAAAAGGTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTtcaat                                                                                                              <  2:1614961/150‑1 (MQ=255)                              ataAGGAGGACACGGTACGACAAAAAAGGTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTtcaatc                                                                                                             <  1:1252653/150‑1 (MQ=255)                              ataAGGAGGACACGGTACGACAAAAAAGGTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTtcaatc                                                                                                             >  2:2956163/1‑150 (MQ=255)                              ataAGGAGGACACGGTACGACAAAAAAGGTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTtcaatc                                                                                                             >  2:2956780/1‑150 (MQ=255)                              ataAGGAGGACACGGTACGACAAAAAAGGTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTtcaatc                                                                                                             >  2:2957163/1‑150 (MQ=255)                               taAGGAGGACACGGTACGACAAAAAAGGTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATCGAGGTtcaatca                                                                                                            >  2:720108/1‑150 (MQ=255)                                    aggaCACGGTACGACAAAAAAGTTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTTGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGt                                                                                                       <  1:3205329/150‑1 (MQ=255)                                     ggaCACGGTACGACAAAAAAGGTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGta                                                                                                      <  1:4309262/150‑1 (MQ=255)                                       acacGGTACGACAAAAAAGGTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATACAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGtaaa                                                                                                    <  2:725436/150‑1 (MQ=255)                                        cacGGTACGACAAAAAAGGTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGtaaat                                                                                                   <  1:747717/150‑1 (MQ=255)                                            gTACGACAAAAAAGGTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGtaaattaaa                                                                                               <  1:2731689/150‑1 (MQ=255)                                            gTACGACAAAAAAGGTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGtaaattaaa                                                                                               <  1:2731711/150‑1 (MQ=255)                                               cGACAAAAAAGGTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAaataa                                                                                            >  2:333396/1‑150 (MQ=255)                                                 aCAAAAAAGGTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAataata                                                                                          <  1:966790/150‑1 (MQ=255)                                                    aaaaaGGTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATACAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAatatat                                                                                       >  1:3064039/1‑150 (MQ=255)                                                     aaaaGGTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATACAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAatatata                                                                                      <  2:4215464/150‑1 (MQ=255)                                                     aaaaGGTAGATGAATCATGTTGAAGATAATACCCGGAGCAACTTGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATACAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAatatata                                                                                      <  2:4215387/150‑1 (MQ=255)                                                         ggTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGt                                                                                  <  1:2928858/150‑1 (MQ=255)                                                              aTGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTATATAATAGATGTTAAATAAAATTTAAATTAAAAAAAATATATTTTAAAta                                                                             >  2:3733714/1‑149 (MQ=255)                                                              aTGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATc                                                                             <  1:487359/150‑1 (MQ=255)                                                              aTGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATc                                                                             <  1:487388/150‑1 (MQ=255)                                                               tGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCt                                                                            <  1:4745607/150‑1 (MQ=255)                                                                gAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTa                                                                           <  1:4266493/150‑1 (MQ=255)                                                                 aaTCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATACAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTAt                                                                          >  2:621471/1‑150 (MQ=255)                                                                  aTCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTAtt                                                                         <  2:1981380/150‑1 (MQ=255)                                                                   tCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTg                                                                        >  1:3934914/1‑150 (MQ=255)                                                                         tGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATACAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATg                                                                  <  2:2992232/150‑1 (MQ=255)                                                                                         tGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGAGTATGCAGCGCTGCGTCCTCg                                                  >  2:2422004/1‑150 (MQ=255)                                                                                          gCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCgg                                                 <  2:3934914/150‑1 (MQ=255)                                                                                           cAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCggg                                                <  2:3030798/150‑1 (MQ=255)                                                                                           cAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCggg                                                <  2:3031916/150‑1 (MQ=255)                                                                                            aaCTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCGGGa                                               >  1:4093721/1‑150 (MQ=255)                                                                                              cTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCGGGata                                             <  1:901427/150‑1 (MQ=255)                                                                                               tGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCGGGatat                                            >  1:280682/1‑150 (MQ=255)                                                                                               tGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCGGGatat                                            >  1:280364/1‑150 (MQ=255)                                                                                               tGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCGGGatat                                            >  1:2315263/1‑150 (MQ=255)                                                                                               tGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATACAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCGGGatat                                            <  1:1168410/150‑1 (MQ=255)                                                                                               tGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATACAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCGGGatat                                            <  1:1168390/150‑1 (MQ=255)                                                                                               tGGATACTTTAATAAAACACTCAATTAAAATCAATTCGATAATGAAGATACAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCGGGatat                                            <  2:2847171/150‑1 (MQ=255)                                                                                                ggATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCGGGAtata                                           <  1:2484795/150‑1 (MQ=255)                                                                                                 gATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCGGGATATaa                                          >  2:21221/1‑150 (MQ=255)                                                                                                   tACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCGGGATATAatt                                        <  1:4265945/150‑1 (MQ=255)                                                                                                      tttAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCGGGATATAattatt                                     <  1:1855570/150‑1 (MQ=255)                                                                                                        taataaAACACTCAATTCAAATCAATTCGATAATGAAGATACAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCGGGATATAATTATTGc                                   >  2:1755682/1‑150 (MQ=255)                                                                                                            aaaaCACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCGGGATATAATTATTGCaaaa                               >  1:3195324/1‑150 (MQ=255)                                                                                                               acacTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTCTTGATTATGCAGCGCTGCGTCATCGGGATATAATTATTGCAAAAATa                            >  2:886315/1‑150 (MQ=255)                                                                                                               acacTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCGGGATATAATTATTGCAAAAATa                            <  1:3814590/150‑1 (MQ=255)                                                                                                                   tCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCGGGATATAATTATTGCAAAAATAGAtt                        <  1:4188166/150‑1 (MQ=255)                                                                                                                   tCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCGGGATATAATTATTGCAAAAATAGAtt                        <  1:3642945/150‑1 (MQ=255)                                                                                                                   tCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCGGGATATAATTATTGCAAAAATAGAtt                        <  1:3641781/150‑1 (MQ=255)                                                                                                                    cAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCGGGATATAATTATTGCAAAAATAGAttt                       >  2:3554342/1‑150 (MQ=255)                                                                                                                     aaTTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCGGGATATAATTATTGCAAAAATAGATTTg                      >  1:4509120/1‑150 (MQ=255)                                                                                                                     aaTTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCGGGATATAATTATTGCAAAAATAGATTTg                      <  2:4635298/150‑1 (MQ=255)                                                                                                                     aaTTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCGGGATATAATTATTGCAAAAATAGATTTg                      >  1:3301542/1‑150 (MQ=255)                                                                                                                              aaaTTCGATAATGAAGATGCAATCAAAGATAAATTATATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCGGGATATAATTATTGCAAAAATAGATTTGTTTATTaaa             >  2:273538/2‑150 (MQ=255)                                                                                                                               aaTTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCGGGATATAATTATTGCAAAAATAGATTTGTTTATTCAAc            >  1:4221348/1‑150 (MQ=255)                                                                                                                                aTTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCGGGATATAATTATTGCAAAAATAGATTTGTTTATTCAACg           >  1:851833/1‑150 (MQ=255)                                                                                                                                  tCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCGGGATATAATTATTGCAAAAATAGATTTGTTTATTCAACGaa         >  1:1778848/1‑150 (MQ=255)                                                                                                                                      tAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCGGGATATAATTATTGCAAAAATAGATTTGTTTATTCAACGAATCAc     <  1:2066194/150‑1 (MQ=255)                                                                                                                                      tAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCGGGATATAATTATTGCAAAAATAGATTTGTTTATTCAACGAATCAc     <  1:2066175/150‑1 (MQ=255)                                                                                                                                       aaTGAAGATGCAAACAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCGGGATATAATTATTGCAAAAATAGATTTGTTTATTCAACGAATCACg    <  1:3733714/150‑1 (MQ=255)                                                                                                                                         tGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCCAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCGGGATATAATTATTGCAAAAATAGATTTGTTTATTCAACGAATCACGCa  >  2:896263/1‑150 (MQ=255)                                                                                                                                         tGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCGGGATATAATTATTGCAAAAATAGATTTGTTTATTCAACGAATCACGCa  >  2:3268593/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                               TTTGATAGTGACCAGATGGCAAAAGAATTATAAGGAGGACACGGTACGACAAAAAAGGTAGATGAATCATGTTGAAGATAATACCCGGTGCAACTGGATACTTTAATAAAACACTCAATTCAAATCAATTCGATAATGAAGATGCAATCAAAGATAAATTAGATAATAGAGGTTCAATCAAAGGTAAATTAAATAATATATATGGTAAATCTATTGATTATGCAGCGCTGCGTCATCGGGATATAATTATTGCAAAAATAGATTTGTTTATTCAACGAATCACGCA  >  CP009273/4250181‑4250466

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