Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	479,096	T→C	5.5%	F41F (TTT→TTC)	USA300HOU_0458 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	479,096	0	T	C	5.5%	226.5 / 3.5	87	F41F (TTT→TTC)	USA300HOU_0458	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (38/41);  new base C (2/3);  total (41/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.02e-01

TATAGGAGGTGAAAGGTATGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGTTACTGATATGTTACTTGAAC  >  CP000730/478957‑479235                                                                                                                                            |                                                                                                                                            tataGGATGTGAAAGGTATGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTg                                                                                                                                            <  2:456326/141‑1 (MQ=255) tataGGAGGTGAAAGGTCTGACAACACAGTTAAATATCCATTCAGTCGTTGAACATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTg                                                                                                                                            <  2:742715/141‑1 (MQ=255) tataGGAGGTGAAAGGTATGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCGGCACGCAATCCATGGGAAGGATTCGAAGCGGATACGTTTg                                                                                                                                            >  1:336877/1‑141 (MQ=255) tataGGAGGTGAAAGGTATGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTg                                                                                                                                            >  1:541615/1‑141 (MQ=255) tagaggaggTGAAAGGTATGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTg                                                                                                                                            <  1:186537/138‑1 (MQ=255)  ataGGAGGTGAAAGGTATGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGa                                                                                                                                           <  1:334142/141‑1 (MQ=255)    aggaggTGAAAGGTATGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAg                                                                                                                                         <  1:106886/141‑1 (MQ=255)     ggaggTGAAAGGTATGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTAAACCATTATCACTAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGa                                                                                                                                        >  2:663102/1‑141 (MQ=255)      gaggTGAAAGGTATGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTGATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGAt                                                                                                                                       >  1:283623/1‑141 (MQ=255)           gAAAGGTATGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGc                                                                                                                                  >  1:607933/1‑141 (MQ=255)             tatgTATGTCAATACAGTTAATGATCAATTCAGTGATTGAAAATGCGAAGCGGGTCATTACACCATTATCACCAATTTCGATTTTTGCAGCACCAAATCCATGTGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCaa                                                                                                                                <  1:339361/138‑1 (MQ=255)             aaGGTATGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCaa                                                                                                                                >  1:682083/1‑141 (MQ=255)             aaGGTATGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATGAGAAGAGGATACGTTTGCAGATGTCGCaa                                                                                                                                >  2:421143/1‑141 (MQ=255)              aGGTATGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGTATACATTTGAAGATGTCGCaaa                                                                                                                               <  2:564199/141‑1 (MQ=255)              aGGTATGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCaaa                                                                                                                               >  1:162523/1‑141 (MQ=255)              aGGTATGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGGGTTATTACACCATTATCACCAATTTCGATGTTTGCAGCACGAAATCCATGGGAAGGAATAGAAGCGGATACGTTTGAAGAGGTAGCaaa                                                                                                                               >  2:104455/1‑141 (MQ=255)                  aTGACAACACGGTTAAATATCACTTGAGTCATTGAAAATGCGAAACGTGTTATTAGACCATTATCACCAAGTTCGTGTTTTGCAGCACGAAATCGATGGTAAGGATTAGAAGCGGATCGGTTTAAAGATGTCGGAAGATgg                                                                                                                           >  2:572119/1‑141 (MQ=255)                  aTGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAGATgg                                                                                                                           >  2:572154/1‑141 (MQ=255)                       aacaCAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACg                                                                                                                      <  1:684174/141‑1 (MQ=255)                          ccccttAAACATCAATTCAGTCATCGAAAATGCGAAACGCGTTACTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGCCTGAAGATGTCGCAAAATGGTTACGTGa                                                                                                                   <  2:573358/137‑1 (MQ=255)                          agaGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACTATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGa                                                                                                                   <  2:686430/139‑1 (MQ=255)                           caGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATTGTTACGTGAt                                                                                                                  >  2:698809/1‑141 (MQ=255)                             gTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGGTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACATTTGAAGATGTCGCAAAAGAGTTACGTGATGt                                                                                                                >  2:482486/1‑141 (MQ=255)                               tAAACATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTc                                                                                                              <  1:89743/141‑1 (MQ=255)                                aaaTATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGttgg                                                                                                             >  2:292631/1‑139 (MQ=255)                                 aaTATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCga                                                                                                            >  2:304540/1‑141 (MQ=255)                                   tatCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCgaga                                                                                                          <  1:402837/141‑1 (MQ=255)                                     tCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACAAATTTCGATTTTTGCAGCACGAAATCCATGGGCAGGATTAGTAGCGGGTACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGAtg                                                                                                        >  1:335010/1‑141 (MQ=255)                                       aaTTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGGtgtg                                                                                                      >  2:498622/1‑141 (MQ=255)                                        gttCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTgg                                                                                                     <  2:378528/140‑1 (MQ=255)                                        aTTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTgg                                                                                                     <  2:495168/141‑1 (MQ=255)                                        aTTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTgg                                                                                                     <  2:263736/141‑1 (MQ=255)                                        aTTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTgg                                                                                                     <  2:246379/141‑1 (MQ=255)                                             gtaaTTGAAAATGGGACACATGTTATTACACCATTATCACCAAGTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTGGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATAtt                                                                                                <  2:147165/138‑1 (MQ=255)                                              tCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATAttt                                                                                               >  1:31908/1‑141 (MQ=255)                                                aTTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGTATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATCtttcc                                                                                             >  2:448408/1‑139 (MQ=255)                                                aTTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTc                                                                                             <  1:545254/141‑1 (MQ=255)                                                 ttGAAAATGCGAAACGTGTTATTCCACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGGTGTTCGAGAGGTGGATATTTTcc                                                                                            >  1:272813/1‑141 (MQ=255)                                                 ttGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTcc                                                                                            <  1:317150/141‑1 (MQ=255)                                                  tGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTccc                                                                                           <  2:676959/141‑1 (MQ=255)                                                   gAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCTATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCa                                                                                          <  1:158257/141‑1 (MQ=255)                                                   gAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCa                                                                                          <  1:34859/141‑1 (MQ=255)                                                     aaaTGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAGATCCATGGGAAGGATTAGAAGCGGATACGTTCGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCaaa                                                                                        >  2:248335/1‑141 (MQ=255)                                                     aaaTGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAGATCCATGGGAAGGATTAGAAGCGGATACGTTCGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCaaa                                                                                        >  2:247916/1‑141 (MQ=255)                                                          cGAAACGTGTTATTACACCATTATCACCAATGTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAGTGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAg                                                                                   >  2:433505/1‑141 (MQ=255)                                                             aaCGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCAt                                                                                <  1:691927/141‑1 (MQ=255)                                                               cGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTa                                                                              >  1:117968/1‑141 (MQ=255)                                                                 tgtTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAAt                                                                            <  1:458673/141‑1 (MQ=255)                                                                     attaCACCGTTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAGTGGTTACGTGATGTTCGAGATGTAGATATTTTCCCAAATAAAGCATTAATAGaa                                                                        >  2:569724/1‑141 (MQ=255)                                                                     attaCACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAAAAGTCACAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAgca                                                                        <  2:389939/141‑3 (MQ=255)                                                                         cacCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAgag                                                                    >  2:315801/1‑139 (MQ=255)                                                                          acCAGTATCACCAACTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAgcgc                                                                   <  1:317556/141‑1 (MQ=255)                                                                              ttATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCtgtg                                                               >  1:108261/1‑141 (MQ=255)                                                                              ttATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAAAAGAAAGCGCtgtg                                                               >  2:300947/1‑141 (MQ=255)                                                                              ttATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGAGGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCtgtg                                                               >  1:607465/1‑141 (MQ=255)                                                                               tATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTgg                                                              >  1:316354/1‑141 (MQ=255)                                                                               tATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTgg                                                              >  2:494917/1‑141 (MQ=255)                                                                               tATCACCAATTTCGATTTTTGCAGCACGAAATCCATGAGAAGGGTTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTgg                                                              >  1:29416/1‑141 (MQ=255)                                                                                       aTTTCGATTTTTGCAGCACGAAATCCCTGGGAAGGGTTAGGAGCGGCTCTATTGGAAGATGTCGCAAAATGGTTACGTGATGTTCGGGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACgtggt                                                      <  1:758235/141‑1 (MQ=255)                                                                                       aTTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATGAGAAGCGGGTTTTCGTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACgtggt                                                      <  1:758267/141‑1 (MQ=255)                                                                                          tCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGaa                                                   >  1:279715/1‑141 (MQ=255)                                                                                            gATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGAGGTGAAtt                                                 >  1:228351/1‑141 (MQ=255)                                                                                             aTTTTTGCAGAACGAAATCCATGGGAAGGATTAGAAGCGGATAGGTTTGAAGATGTGGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCAGGTGATGAATTa                                                >  2:300504/1‑141 (MQ=255)                                                                                                tttGCAGCACGAGATCCATGGGAAGGATTAGAAGCGGATACGTTCGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGAt                                             <  1:248335/141‑1 (MQ=255)                                                                                                tttGCAGCACGAGATCCATGGGAAGGATTAGAAGCGGATACGTTCGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGAt                                             <  1:247916/141‑1 (MQ=255)                                                                                                tttGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGAt                                             <  1:485593/141‑1 (MQ=255)                                                                                                  tGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGa                                           <  2:232873/141‑1 (MQ=255)                                                                                                  tGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGa                                           >  1:755246/1‑141 (MQ=255)                                                                                                  tGCAGCACGAAATCCATGAGAAGGCTTAAAAACGGATCCAATAGCAGATGTCACAAAATGATTATGAGATACTCGAGAAATGTATAGTTGCCCAAGCAAAGCCTTAAGAGAAAGCGCGTTGGCACGTGGTGAATTAGATGa                                           >  2:181358/1‑141 (MQ=255)                                                                                                     agaaCGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGGTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCTCTGTGGCACGTGGTGAATTAGATGAAAg                                        <  2:335010/138‑1 (MQ=255)                                                                                                       caCGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAgtg                                      <  1:745197/141‑1 (MQ=255)                                                                                                                caggggAAGGATTAGGAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATc                             <  2:301155/138‑1 (MQ=255)                                                                                                                cATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGCAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATc                             <  1:744336/141‑1 (MQ=255)                                                                                                                       agtaTTAGAAGCGGATACGTTTGAAGATGTCGTAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGt                      <  1:330938/138‑1 (MQ=255)                                                                                                                        ggATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGtt                     >  2:719053/1‑141 (MQ=255)                                                                                                                            tAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGTTACTg                 <  2:633197/141‑1 (MQ=255)                                                                                                                            tAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGTTACTg                 <  1:114250/141‑1 (MQ=255)                                                                                                                               gaGCGGATACGTTTGAAGATGTCGAAAAATGGTTACGTGTTGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGTTACTGata              <  1:131148/140‑1 (MQ=255)                                                                                                                                  gggATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGCTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGTACGGTGTGAATTAGATGAAAGTGCATTTAATCAACTTGTTACTGATATGt           <  2:23068/140‑1 (MQ=255)                                                                                                                                  cGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGTTACTGATATGt           <  1:609660/141‑1 (MQ=255)                                                                                                                                  cGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGTTACTGATATGt           <  1:421143/141‑1 (MQ=255)                                                                                                                                   ggATACGTTTGAAGGTGTCGCAAAATGGTGACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGTTACTGATATGtt          <  2:39068/141‑1 (MQ=255)                                                                                                                                    gATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGTTACTGATATGTTa         <  2:206159/141‑1 (MQ=255)                                                                                                                                      tACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGTTACTGATATGTTaat       >  2:725948/1‑139 (MQ=255)                                                                                                                                      tACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGTTACTGATATGTTACt       >  1:304581/1‑141 (MQ=255)                                                                                                                                        cGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGTTACTGATATGTTACTTg     >  2:736363/1‑141 (MQ=255)                                                                                                                                         gTTCGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGTTACTGATATGTTACTTGa    <  1:45913/141‑1 (MQ=255)                                                                                                                                           ttGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGTTACTGATATGTTACTTGAAc  <  2:607465/141‑1 (MQ=255)                                                                                                                                           ttGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGTTACTGATATGTTACTTGAAc  <  2:108261/141‑1 (MQ=255)                                                                                                                                           ttGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTACCAAAGAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGCTACTGATATGTTGCTTGaaa  >  2:600294/1‑140 (MQ=255)                                                                                                                                            |                                                                                                                                            TATAGGAGGTGAAAGGTATGACAACACAGTTAAATATCAATTCAGTCATTGAAAATGCGAAACGTGTTATTACACCATTATCACCAATTTCGATTTTTGCAGCACGAAATCCATGGGAAGGATTAGAAGCGGATACGTTTGAAGATGTCGCAAAATGGTTACGTGATGTTCGAGATGTGGATATTTTCCCAAATAAAGCATTAATAGAAAGCGCTGTGGCACGTGGTGAATTAGATGAAAGTGTCTTTAATCAACTTGTTACTGATATGTTACTTGAAC  >  CP000730/478957‑479235

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