Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	2,608,825	A→T	6.2%	intergenic (‑676/+154)	USA300HOU_2470 ← / ← USA300HOU_2471	dehydrogenase/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,608,825	0	A	T	6.2%	236.5 / 6.2	77	intergenic (‑676/+154)	USA300HOU_2470/USA300HOU_2471	dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (31/40);  new base T (3/2);  total (34/43)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

ATTTTAAACAAGTTGAAAATGAAGGATTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTGCCGCCGACTAGGGGTAGTAAT  >  CP000730/2608688‑2608963                                                                                                                                          |                                                                                                                                           aTTTTAAACAAGTTGCAAATGAAGGATTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTATAGCCACATCCAATATAGGa                                                                                                                                         >  2:78999/1‑141 (MQ=255) aTTTTAAACAAGTTGAAAATGAAGGATTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGa                                                                                                                                         >  2:674796/1‑141 (MQ=255)  ttatAAACAAGTTGAAAATGAAGGATTCTACTATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGAc                                                                                                                                        >  2:270649/4‑141 (MQ=255)   tttAAACNAGTTGAAAATGAAGGATTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATCGGACg                                                                                                                                       <  2:590769/141‑1 (MQ=255)   tttAAACAAGTTGAAAATGAAGGATTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACg                                                                                                                                       <  1:11575/141‑1 (MQ=255)     tAAACAAGTTGAAAATGAAGGATTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACGTCCTATATTGGAtgtg                                                                                                                                     <  1:534484/141‑1 (MQ=255)       aaCAAGTTGAAAATGAAGGATTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGc                                                                                                                                   >  2:522174/1‑141 (MQ=255)       aaCAAGTTGAAAATGAAGGATTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGGGGc                                                                                                                                   >  2:665796/1‑141 (MQ=255)          aaGTTGAAAATGAAGGATTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCttt                                                                                                                                <  2:89239/141‑1 (MQ=255)            gTTGAAAATGAAGGATTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTa                                                                                                                              >  2:66545/1‑141 (MQ=255)            gTTGAAAATGAAGGATTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTa                                                                                                                              <  2:402660/141‑1 (MQ=255)            gTTGAAAATGAAGGATTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTa                                                                                                                              >  2:666895/1‑141 (MQ=255)            gTTGAAAATGAAGGATTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTa                                                                                                                              <  2:408998/141‑1 (MQ=255)             ttGAAAATGAAGGATTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAg                                                                                                                             >  2:138776/1‑141 (MQ=255)               gAAAATGAAGGATTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTa                                                                                                                           >  1:119398/1‑141 (MQ=255)               gAAAATGAAGGATTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTa                                                                                                                           >  2:262781/1‑141 (MQ=255)               gAAAATGAAGGATTCAAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACCTgcaagtaaatg                                                                                                                           >  2:228093/1‑131 (MQ=255)                aaaaTGAAGGATTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTAt                                                                                                                          <  1:465579/141‑1 (MQ=255)                   aTGAAGGATTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGt                                                                                                                       >  1:502786/1‑141 (MQ=255)                        ggATTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAg                                                                                                                  <  2:8736/141‑1 (MQ=255)                          aTTCTAATATCATTTCGTTTATATATATCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGaa                                                                                                                >  2:668327/1‑140 (MQ=255)                             cTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTtg                                                                                                             <  2:465739/141‑1 (MQ=255)                                atatCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGAAATAAAAAAAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTgg                                                                                                          <  2:512677/141‑1 (MQ=255)                                  atCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATACAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTgggg                                                                                                        <  1:768545/141‑1 (MQ=255)                                     aTTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGt                                                                                                     <  2:484834/141‑1 (MQ=255)                                     aTTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGt                                                                                                     <  2:714812/141‑1 (MQ=255)                                              atatataGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTa                                                                                            <  1:775229/141‑1 (MQ=255)                                               tatataGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTGAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTAc                                                                                           <  1:768431/141‑1 (MQ=255)                                                atataGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACt                                                                                          <  2:539045/141‑1 (MQ=255)                                                atataGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACt                                                                                          <  2:546951/141‑1 (MQ=255)                                                  ataGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACtta                                                                                        <  1:502887/141‑1 (MQ=255)                                                  ataGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACtta                                                                                        <  1:25734/141‑1 (MQ=255)                                                  ataGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACtta                                                                                        <  2:531693/141‑1 (MQ=255)                                                    aGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaat                                                                                      >  1:739477/1‑141 (MQ=255)                                                             gATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGt                                                                             >  2:264736/1‑141 (MQ=255)                                                                 gAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGc                                                                         <  2:158571/141‑1 (MQ=255)                                                                    ttttATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCt                                                                      <  2:768873/141‑1 (MQ=255)                                                                    ttttATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCt                                                                      >  2:437223/1‑141 (MQ=255)                                                                      ttATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGTATTGTGGGGAGTTAGGGTTTACTTAATTAAATGGGTCTTCCTTCa                                                                    >  2:180017/1‑141 (MQ=255)                                                                      ttATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCa                                                                    <  2:282798/141‑1 (MQ=255)                                                                       tataTGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCaa                                                                   >  2:267755/1‑141 (MQ=255)                                                                        atatGTAAATCTTGTGGGTAATCGTGTTAAAAATAATATAAGTATGACTAAAGCACCATCCAATATAGTGCGACGCTTGTGGGGTGGTTCTTGATGTGGGGGGAGTTAGGGTGTACTTAATTATAGGGGTCTGCCCGCagt                                                                  >  1:196385/1‑139 (MQ=255)                                                                        atatGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAAt                                                                  >  1:420124/1‑141 (MQ=255)                                                                         tatGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAAtt                                                                 <  1:75620/141‑1 (MQ=255)                                                                           tggaaaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGTGTCTGCCCTCAATTGc                                                               <  2:556675/138‑1 (MQ=255)                                                                             tAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTTGCCTCAATTGCTc                                                             >  2:594178/1‑141 (MQ=255)                                                                              aaaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCa                                                            <  2:169893/141‑1 (MQ=255)                                                                              aaaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCa                                                            <  1:138776/141‑1 (MQ=255)                                                                                  cTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCaa                                                        <  2:518103/141‑1 (MQ=255)                                                                                  cTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCaa                                                        <  2:571201/141‑1 (MQ=255)                                                                                   ttGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCTAATATTGGACGTGGCTTTTAGTATCGTTCCCGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAAc                                                       >  1:377449/1‑141 (MQ=255)                                                                                     gTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAACTGCTCACCAACGa                                                     >  2:365700/1‑141 (MQ=255)                                                                                       aGGTAATCGTTTTAGAAATAATATAAGTATGACTAAAGCCACATCCAATAGAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAgg                                                   <  2:480268/141‑1 (MQ=255)                                                                                          tAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCa                                                >  2:716643/1‑141 (MQ=255)                                                                                          tAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCa                                                <  1:650002/141‑1 (MQ=255)                                                                                              cGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGtt                                            >  2:96095/1‑141 (MQ=255)                                                                                               gTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTa                                           <  1:169860/141‑1 (MQ=255)                                                                                                 tttAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACa                                         <  1:522723/141‑1 (MQ=255)                                                                                                  ttAAAAATAATATAAGTATGACTAACGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACat                                        <  1:416600/141‑1 (MQ=255)                                                                                                  ttAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACat                                        >  2:478556/1‑141 (MQ=255)                                                                                                   tAAAAATAATATAAGTATGACTAAAGCCACATTCAATATAGGACGTGGCTTTTAGTATCGTCCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACatc                                       <  1:137032/141‑1 (MQ=255)                                                                                                        ataataTAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACatcaatca                                  >  1:401811/1‑141 (MQ=255)                                                                                                          aataTAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATcaca                                <  1:29471/141‑1 (MQ=255)                                                                                                             ggaacgatactaaaaGCCACGTCCTATATTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACtt                             <  2:773536/129‑1 (MQ=37)                                                                                                              taAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTc                            >  1:293838/1‑141 (MQ=255)                                                                                                                aGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCaa                          <  1:348816/141‑1 (MQ=255)                                                                                                                aGTATGACTAAAGCCACAACCAGTACACGTAGGGGCGTTGAGTATCTTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCaa                          >  2:626090/1‑141 (MQ=255)                                                                                                                  tATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAAtt                        <  1:463895/141‑1 (MQ=255)                                                                                                                      aCTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCACTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTgccg                    >  1:146718/1‑141 (MQ=255)                                                                                                                        tAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTgccgcc                  >  2:433795/1‑141 (MQ=255)                                                                                                                           aGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTGCCGCCGAc               <  1:342623/141‑1 (MQ=255)                                                                                                                           aGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTGCCGCCGAc               <  1:342555/141‑1 (MQ=255)                                                                                                                                  tCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTGCCGCCGACTAGGGgt        <  2:366809/141‑1 (MQ=255)                                                                                                                                   ccAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTGCCGCCGACTAGGGgta       >  2:752337/1‑141 (MQ=255)                                                                                                                                   ccAAGATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGG‑TTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTGACATCAATCACACTTCAATTGCCGCCGACTAGGGgtag      <  2:10678/141‑1 (MQ=255)                                                                                                                                     tatatTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTGCCGCCGACTAGGGgtagt     >  1:192706/2‑141 (MQ=255)                                                                                                                                     tatatTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTGCCGCCGACTAGGGgtagt     >  1:289622/2‑141 (MQ=255)                                                                                                                                        ataGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTGCCGCCGACTAGGGGTAGTAAt  <  2:543176/141‑1 (MQ=255)                                                                                                                                          |                                                                                                                                           ATTTTAAACAAGTTGAAAATGAAGGATTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTGCCGCCGACTAGGGGTAGTAAT  >  CP000730/2608688‑2608963

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