Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	2,120,639	A→G	9.0%	L46L (TTA→CTA)	USA300HOU_1994 ←	hypothetical bacteriophage protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,120,639	0	A	G	9.0%	242.7 / 13.1	89	L46L (TTA→CTA)	USA300HOU_1994	hypothetical bacteriophage protein
Reads supporting (aligned to +/- strand):  ref base A (45/36);  new base G (5/3);  total (50/39)
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.47e-01

TCTCATGCTTATACCTACCATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTT  >  CP000730/2120500‑2120778                                                                                                                                            |                                                                                                                                             tctcATGCTTATACCTACCATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTACAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAAc                                                                                                                                             <  1:784846/141‑1 (MQ=255) tctcATGCTTATACCTACCATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAAc                                                                                                                                             <  1:989805/141‑1 (MQ=255) tctcATGCTTATACCTACCATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAAc                                                                                                                                             >  2:700825/1‑141 (MQ=255) tctcATGCTTATACCTACCATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAAc                                                                                                                                             <  1:61745/141‑1 (MQ=255)  ctcATGCTTATACCTACCATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAAcc                                                                                                                                            >  1:748148/1‑141 (MQ=255)    cATGCTTATACCTACCATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCg                                                                                                                                          >  2:160075/1‑141 (MQ=255)     aTGCTTATACCTACCATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGt                                                                                                                                         >  2:909957/1‑141 (MQ=255)           ataCCTACCATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTg                                                                                                                                   >  1:522359/1‑141 (MQ=255)            tacctaccATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGa                                                                                                                                  <  1:238647/141‑1 (MQ=255)            tacctaccATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACCAGTAACCCGTCAATTGa                                                                                                                                  >  1:603800/1‑141 (MQ=255)             acctaccATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCNCATCAGAAACTANTAACCCGTCAATTGaa                                                                                                                                 <  1:809794/141‑1 (MQ=255)               ctaccATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAAtt                                                                                                                               <  1:449886/141‑1 (MQ=255)                    aTTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCtt                                                                                                                          <  2:1030522/141‑1 (MQ=255)                     tttCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTc                                                                                                                         <  2:890044/141‑1 (MQ=255)                       tCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCAt                                                                                                                       <  2:385348/141‑1 (MQ=255)                        cATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATg                                                                                                                      <  1:886707/141‑1 (MQ=255)                        cATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATg                                                                                                                      >  2:723607/1‑141 (MQ=255)                         aTGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATgg                                                                                                                     >  2:587370/1‑141 (MQ=255)                            aCTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCt                                                                                                                  <  1:592804/141‑1 (MQ=255)                            aCTAAGTTAATTAGTCTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAATACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCGTCAGAAACTAGTAGCCCGTCAATTGAATTTCCTTCATGATCt                                                                                                                  <  2:1076496/141‑1 (MQ=38)                                aGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGaa                                                                                                              <  1:503076/141‑1 (MQ=255)                                 gttaaTTAGTCTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAATACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCGTCAGAAACTAGTAGCCCGTCAATTGAATTTCCTTCATGATCTAAaac                                                                                                             >  2:1080817/5‑141 (MQ=38)                                   tCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCGTAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGaacaa                                                                                                           <  1:852966/141‑1 (MQ=255)                                    cATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGaacaac                                                                                                          >  1:714974/1‑141 (MQ=255)                                        aGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCGGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCAtt                                                                                                      <  1:281820/141‑1 (MQ=255)                                         gTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCAttt                                                                                                     >  1:872486/1‑141 (MQ=255)                                              tCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACg                                                                                                >  2:1065026/1‑141 (MQ=255)                                               ccTGTTCATCTGTGTTATTTTCAATCCATTCATAAATACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCGTCAGAAACTAGTAGCCCGTCAATTGAATTTCCTTCATGATCTAAAACGACTATTTCGAcac                                                                                               >  2:705909/1‑139 (MQ=255)                                               ccTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGc                                                                                               >  1:663286/1‑141 (MQ=255)                                               ccTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGc                                                                                               >  2:760113/1‑141 (MQ=255)                                                   ttCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCtata                                                                                           >  2:538008/1‑141 (MQ=255)                                                   ttCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCtata                                                                                           >  2:923004/1‑141 (MQ=255)                                                   ttCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCtata                                                                                           >  1:644608/1‑141 (MQ=255)                                                    tCATCTGTGTTATTTTCAATCCATTCATAAATACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCGTCAGAAACTAGTAGCCCGTCAATTGAATTTCCTTCATGATCTAAAACGACTATTTCGACACtatat                                                                                          >  1:526314/1‑140 (MQ=255)                                                    tCATCTGTGTTATTTTCAATCCATTCATAAATACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCGTCAGAAACTAGTAGCCCGTCAATTGAATTTCCTTCATGATCTAAAACGACTATTTCGACACtatat                                                                                          >  1:734959/1‑140 (MQ=255)                                                    tCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATAc                                                                                          <  2:1003342/141‑1 (MQ=255)                                                        cTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTc                                                                                      <  1:429292/141‑1 (MQ=255)                                                          gtgtTATTTTCAATCCATTCATAAATACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCGTCAGAAACTAGTAGCCCGTCAATTGAATTTCCTTCATGATCTAAAACGACTATTTCGACACTATATGCTCgc                                                                                    <  1:455632/141‑2 (MQ=255)                                                          gtgtTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGt                                                                                    >  2:635553/1‑141 (MQ=255)                                                             ttATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTc                                                                                 <  2:629719/141‑1 (MQ=255)                                                               aTTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCtt                                                                               <  1:299598/141‑1 (MQ=255)                                                                    cAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAActc                                                                          <  1:309741/141‑1 (MQ=255)                                                                    cAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAActc                                                                          >  1:372840/1‑141 (MQ=255)                                                                     aaTCCATTCATAAATACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCGTCAGAAACTAGTAGCCCGTCAATTGAATTTCCTTCATGATCTAAAACGACTATTTCGACACTATATGCTCGCTTCTTAActct                                                                         <  2:82870/141‑1 (MQ=255)                                                                      aTCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCtt                                                                        <  1:846708/141‑1 (MQ=255)                                                                       tCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTa                                                                       >  1:992887/1‑141 (MQ=255)                                                                       tCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTa                                                                       >  1:475000/1‑141 (MQ=255)                                                                           ttCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTg                                                                   <  1:641314/141‑1 (MQ=255)                                                                               tAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑aaaa                                                              >  2:1046637/1‑141 (MQ=255)                                                                                 aaCACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATc                                                            <  2:730749/141‑1 (MQ=255)                                                                                   cacTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGGTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCaa                                                          >  1:432525/1‑141 (MQ=255)                                                                                    acTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCaaa                                                         >  2:188678/1‑141 (MQ=255)                                                                                       ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGcc                                                      >  1:940092/1‑141 (MQ=255)                                                                                       ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGcc                                                      >  2:101574/1‑141 (MQ=255)                                                                                       ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGcc                                                      >  1:1050983/1‑141 (MQ=255)                                                                                       ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGcc                                                      >  2:1054135/1‑141 (MQ=255)                                                                                       ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGcc                                                      >  1:219361/1‑141 (MQ=255)                                                                                       ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGcc                                                      >  2:944847/1‑141 (MQ=255)                                                                                         gTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCAt                                                    <  1:929619/141‑1 (MQ=255)                                                                                               aaaTATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAAATCAAAGCCATCTACa                                               >  1:897473/1‑141 (MQ=255)                                                                                                 atatCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACatta                                            <  1:460266/141‑1 (MQ=255)                                                                                                 atatCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACatta                                            <  2:669992/141‑1 (MQ=255)                                                                                                  tatCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACattaa                                           <  1:797050/141‑1 (MQ=255)                                                                                                    tCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACattaatt                                         <  1:631767/141‑1 (MQ=255)                                                                                                    tCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACattaatt                                         <  2:46529/141‑1 (MQ=255)                                                                                                       aaaGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAattatt                                      <  1:522513/141‑1 (MQ=255)                                                                                                        aaGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTAttt                                     >  2:1058932/1‑141 (MQ=255)                                                                                                         aGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTAtttt                                    <  1:576778/141‑1 (MQ=255)                                                                                                          gCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTAttttt                                   >  1:712845/1‑141 (MQ=255)                                                                                                             gtgtATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGa                                >  2:409349/1‑141 (MQ=255)                                                                                                                tataGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACAt                             <  1:439432/141‑1 (MQ=255)                                                                                                                        gTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCAcc                     <  1:617052/141‑1 (MQ=255)                                                                                                                         ttCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCAccc                    >  2:60905/1‑141 (MQ=255)                                                                                                                          tctcATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCg                   >  2:625438/1‑141 (MQ=255)                                                                                                                          tctcATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCg                   >  2:883004/1‑141 (MQ=255)                                                                                                                          tctcATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCg                   >  2:927714/1‑141 (MQ=255)                                                                                                                           ctcATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGt                  >  2:723508/1‑141 (MQ=255)                                                                                                                           ctcATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGt                  >  2:723512/1‑141 (MQ=255)                                                                                                                            tcgtcaGAAACTAGTAGCCCGTCAATTGAATTTCCTTCATGATCTAAAACGACTATTTCGACACTATATGCTCGCTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAAATATTTTTCGACATACGTCACCCGtt                 >  1:22892/4‑141 (MQ=255)                                                                                                                              atcaGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGtttt               >  1:29322/1‑141 (MQ=255)                                                                                                                                caGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGt             <  1:34641/141‑1 (MQ=255)                                                                                                                                 aGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCTTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTa            >  1:340122/1‑141 (MQ=255)                                                                                                                                  gAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTaa           <  2:379579/141‑1 (MQ=255)                                                                                                                                       tagtagCCCGTCAATTGAATTTCCTTCATGATCTAAAACGACTATTTCGACACTATATGCTCGCTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAAATATTTTTCGACATACGTCACCCGTTTTGTAATACAt      >  1:715115/1‑141 (MQ=255)                                                                                                                                       tAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCAcc                     >  2:1074934/1‑126 (MQ=255)                                                                                                                                       tAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCAcc                     <  1:1074934/126‑1 (MQ=255)                                                                                                                                        aGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACAtt     >  1:1014721/1‑141 (MQ=255)                                                                                                                                        aGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACAtt     <  1:923004/141‑1 (MQ=255)                                                                                                                                        aGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACAtt     >  2:993461/1‑141 (MQ=255)                                                                                                                                         gTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTg    >  2:176986/1‑141 (MQ=255)                                                                                                                                           aaCCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGtt  >  1:26639/1‑141 (MQ=255)                                                                                                                                            |                                                                                                                                             TCTCATGCTTATACCTACCATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTT  >  CP000730/2120500‑2120778

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