Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	1,211,770	A→G	93.4%	A22A (GCA→GCG)	U712_06040 →	Spore coat protein O

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	1,211,770	0	A	G	93.4%	228.2 / 5.4	76	A22A (GCA→GCG)	U712_06040	Spore coat protein O
Reads supporting (aligned to +/- strand):  ref base A (2/3);  new base G (37/34);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.50e-01

CTTTTTTATTTCTTCTGATTTTCAGCTTTCTGTCATATAGATAGAATATGACACAATCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCAGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCTCCGCCAAAAGCTGAGGAAAAGCCGGCTTACCATGAGAAGCAGCAGGAGGAAGCCGTGGCACAGGATGTCC  >  CP006881/1211635‑1211900                                                                                                                                        |                                                                                                                                   cTTTTTTATTTCTTCTGATTTTCAGCTTTCTGTCTTATAGATAGAATATGACACAATCGGAGGGATGAAGATGTATAATAAGAAAAGAAACGCGGATAAAAAACGGCTCATGTATATTGTTCAGCCGGATTATGCGGa                                                                                                                                  <  3:151902/138‑1 (MQ=255) cTTTTTTATTTCTTCTGATTTTCAGCTTTCTGTCATATAGATAGAATATGACACAATCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGa                                                                                                                                  <  5:357641/138‑1 (MQ=255) cTTTTTTATTTCTTCTGATTTTCAGCTTTCTGTCATATAGATAGAATATGACACAATCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGa                                                                                                                                  <  7:402865/138‑1 (MQ=255) cTTTTTTATTTCTTCTGATTTTCAGCTTTCTGTCATATAGATAGAATATGACACAATCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGa                                                                                                                                  <  7:17645/138‑1 (MQ=255) cTTTTTTATTTCTTCTGATTTTCAGCTTTCTGTCATATAGATAGAATATGACACAATCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGa                                                                                                                                  <  5:168116/138‑1 (MQ=255) cTTTTTTATTTCTTCTGATTTTCAGCTTTCTGTCATATAGATAGAATATGACACAATCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGa                                                                                                                                  <  5:160227/138‑1 (MQ=255) cTTTTTTATTTCTTCTGATTTTCAGCTTTCTGTCAAAAACATAGAATATGACACAATCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGa                                                                                                                                  <  2:443331/138‑1 (MQ=255) cTTTTTTATTTCTTCTGATTTTCAGCTTTATGTCATATAGATAGAATATGACACAATCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCAGa                                                                                                                                  <  1:187233/138‑1 (MQ=255)           tcttctGATTTTCAGCTTTCTGTCATATAGATAGAATATGACACAATCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGa                                                                                                                                  <  4:57744/128‑1 (MQ=255)           tcttctGATTTTCAGCTTTCTGTCATATAGATAGAATATGACACAATCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGa                                                                                                                                  >  3:57744/1‑128 (MQ=255)           tcttctGATTTTCAGCTTTCTGTCATATAGATAGAATATGACACAATCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGAAACTACAAGGt                                                                                                                       >  7:399681/1‑139 (MQ=255)            cttctGATTTTCAGCTTTCTGTCATATAGATAGAATATGACACAATCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCAGAAACTACAAGGt                                                                                                                       <  7:374291/138‑1 (MQ=255)               ctGATTTTCAGCTTTCTGTCATATAGATAGAATATGACACAATCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGAAACTACAAGGTCAt                                                                                                                    <  6:230349/138‑1 (MQ=255)                  aTTTTCAGCTTTCTGTCATATAGATAGAATATGACACAATCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGAAACTACAAGGTCATCCAt                                                                                                                >  2:37710/1‑139 (MQ=255)                       cAGCTTTCTGTCATATAGATAGAATATGACACAATCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAACACCGCTCATGTATATTGTTCAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAg                                                                                                            >  5:166235/1‑138 (MQ=255)                          cTTTCTGTCATATAGATAGAATATGACACAATCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGgaga                                                                                                        <  1:37710/139‑1 (MQ=255)                                 tCATATAGATAGAATATGACACAATCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCAGAAACTACAAGGTCATCCATGCAGGAGATTGtcat                                                                                                 <  6:31631/139‑1 (MQ=255)                                  cATATAGATAGAATATGACACAATCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGtcaa                                                                                                 <  4:33216/138‑2 (MQ=255)                                         atagaATATGACACAAACGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCAAGTATATTGTTCAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAaaaga                                                                                          >  3:114730/1‑138 (MQ=255)                                           agaATATGACACAATCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGa                                                                                                                                  <  5:89677/96‑1 (MQ=255)                                           agaATATGACACAATCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGa                                                                                                                                  >  6:89677/1‑96 (MQ=255)                                           agaATATGACACAATCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAaaagaacg                                                                                       >  6:165092/1‑137 (MQ=255)                                             aaTATGACACAATCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGc                                                                                     >  2:241762/1‑139 (MQ=255)                                              atatGACACAATCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGcc                                                                                    <  2:350723/139‑1 (MQ=255)                                                  gACACAATCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACAGCTCATGTATATTGTTCAGCCGGATTATGAGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGAc                                                                                 >  1:231852/1‑138 (MQ=255)                                                   acacAATCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGATACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACa                                                                                >  7:100191/1‑138 (MQ=255)                                                   acacAATCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACaa                                                                               >  8:384447/1‑139 (MQ=255)                                                   acacAATCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACa                                                                                >  3:30408/1‑138 (MQ=255)                                                       aaTCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGa                                                                                                                                  >  2:37539/1‑84 (MQ=255)                                                       aaTCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGa                                                                                                                                  <  1:37539/84‑1 (MQ=255)                                                       aaTCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCTACAAAccc                                                                           >  3:136286/1‑139 (MQ=255)                                                       aaTCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAAccc                                                                           >  4:153387/1‑139 (MQ=255)                                                       aaTCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAAccc                                                                           >  3:51465/1‑139 (MQ=255)                                                       aaTCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCAGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAAc                                                                             >  4:81059/1‑137 (MQ=255)                                                           ggAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCt                                                                        <  6:320398/138‑1 (MQ=255)                                                           ggAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCt                                                                        <  1:241762/138‑1 (MQ=255)                                                            gAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGa                                                                                                                                  >  8:289753/1‑79 (MQ=255)                                                            gAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGa                                                                                                                                  <  7:289753/79‑1 (MQ=255)                                                             aGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAACGGCCGACAAACCCGCTCCg                                                                     >  2:439781/1‑139 (MQ=255)                                                                      gATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCTCCGCCAAAAGCt                                                            <  3:153387/139‑1 (MQ=255)                                                                          tCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGAAACTACAAGGt                                                                                                                       >  7:318554/1‑76 (MQ=255)                                                                          tCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGAAACTACAAGGt                                                                                                                       <  8:318554/76‑1 (MQ=255)                                                                          tCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCTCCGCCAAAAGCTGAg                                                         >  5:436425/1‑138 (MQ=255)                                                                          tCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGCAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCTCCGCCAAAAGCTGAg                                                         >  5:64179/1‑138 (MQ=255)                                                                               taaCAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCTCCGCCAAAAGCTGAGGAAAAg                                                   <  4:51465/139‑1 (MQ=255)                                                                                 agaaaagaaaCGCGGATTAAATACCGCTCATGTATATTTTTAAGTCGGTTTATTCGGAAACTACAAGGTCATCAATTCAAGAGATAGTCATCAAAAGAAAGGCCGACAAACCCGCTACGCCAAAAGCTCAGGAAAAGcc                                                 >  5:194403/1‑139 (MQ=255)                                                                                 agaaaagaaaCGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGa                                                                                                                                  >  5:209475/1‑58 (MQ=255)                                                                                 agaaaagaaaCGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGa                                                                                                                                  <  6:209475/58‑1 (MQ=255)                                                                                    aaagaaaCGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCTCCGCCAAAAGCTGAGGAAAAGCCGGc                                              <  6:436425/139‑1 (MQ=255)                                                                                    aaagaaaCGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCTCCGCCAAAAGCTGAGGAAAAGCAGGc                                              >  6:212187/1‑139 (MQ=255)                                                                                     aagaaaCGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCTCCGCCAAAAGCTGAGGAAAAGCCGGCt                                             >  6:110251/1‑139 (MQ=255)                                                                                               taTAAAAAACCGCTCATGTATTTTTTTCAGCCGTATTATGCGGAAACTACAAGGTCATACATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCTACGCCAAAATGTGAGGAAAAGCCTGCTTACCATgaga                                   <  4:114730/138‑1 (MQ=255)                                                                                               gATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCTCCGCCAAAAGCTGAGGAAAAGCCGGCTTACCATgaga                                   <  8:399681/139‑1 (MQ=255)                                                                                               gATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGAGGAAACTACAAGGTGATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCTCCGCCAAAAGCTGAGGAAAAGCCGGCTTACCATgaga                                   <  6:228091/139‑1 (MQ=255)                                                                                                    aaaaCCGCTCATGTATATTGTTCAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCTCCGCCAAAAGCTGAGGAAAAGCCGGCTTACCATGAGAagcag                              <  3:442180/139‑1 (MQ=255)                                                                                                          gCTCATGTATATTGTTCAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCTCCGCCAAAAGCTGAGGAAAAGCCGGCTTACCATGAGAAGCAGCaggag                        >  7:16360/1‑139 (MQ=255)                                                                                                          gCTCATGTATATTGTTCAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCTCCGCCAAAAGCTGAGGAAAAGCCGGCTTACCATGAGAAGCAGCaggag                        <  5:110251/139‑1 (MQ=255)                                                                                                          gCTCATGTATATTGTTCAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCTCCGCCAAAAGCTGAGGAAAAGCCGGCTTACCATGAGAAGCAGCaggag                        >  5:102525/1‑139 (MQ=255)                                                                                                          gCTCATGTATATTGTTCAGCCGGATTATGCAGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCTCCGCCAAAAGCTGAGGAAAAGCCGGCTTACCATGAGAAGCAGCaggag                        >  5:221106/1‑139 (MQ=255)                                                                                                              aTGTATATTGTTCAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCTCCGCCAAAAGCTGAGGAAAAGCCGGCTTACCATGAGAAGCAGCAggagga                      <  4:75068/137‑1 (MQ=255)                                                                                                                 tataTTGTTCAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCTCCGCCAAAAGCTGAGGAAAAGCCGGCTTACCATGAGAAGCAGCAGGAGGAAGCCg                 <  1:63959/139‑1 (MQ=255)                                                                                                                 tataTTGTTCAGCCGGATTATGCGGAAACTACAACGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCTCCGCCAAAAGCTGAGGAAAAGCCGGCTTACCATGAGAAGCAGCAGGAGGAAGCCg                 <  6:294491/139‑1 (MQ=255)                                                                                                                      tGTTAAGCAGTATTATGCGGATACTACAAGGTCATCCATGCATGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCTCCGCCAAAAGCTGAGGAAAAGCCGGCTTACCATGAGAAGCAGCAGGAGGAAGCCGTGGc             <  2:231852/138‑1 (MQ=255)                                                                                                                        ttCAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCTCCGCCAAAAGCTGAGGAAAAGCCGGCTTACCATGAGAagcagc                             >  7:142520/1‑120 (MQ=255)                                                                                                                        ttCAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCTCCGCCAAAAGCTGAGGAAAAGCCGGCTTACCATGAGAagcagc                             <  8:142520/120‑1 (MQ=255)                                                                                                                        ttCAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCTCCGCCAAAAGCTGAGGAAAAGCCGGCTTACCATGAGAAGCAGCAGGAGGAAGCCGTGGca            >  5:406558/1‑137 (MQ=255)                                                                                                                        ttCAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCTCCGCCAAAAGCTGAGGAAAAGCCGGCTTACCATGAGAAGCAGCAGGAGGAAGCCGTGGca            >  8:25518/1‑137 (MQ=255)                                                                                                                        ttCAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCTCCGCCAAAAGCTGAGGAAAAGCCGGCTTACCATGAGAAGCAGCAGGAGGAAGCCGTGGca            >  1:301464/1‑137 (MQ=255)                                                                                                                          cAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCTCCGCCAAAAGCTGAGGAAAAGCCGGCTTACCATGAGAAGCAGCAGGAGGAAGCCGTGGc             >  2:278044/1‑134 (MQ=255)                                                                                                                          cAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCTCCGCCAAAAGCTGAGGAAAAGCCGGCTTACCATGAGAAGCAGCAGGAGGAAGCCGTGGc             <  1:278044/134‑1 (MQ=255)                                                                                                                          cAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCTCCGCCAAAAGCTGAGGAAAAGCCGGCTTACCATGAGAAGCAGCAGGAGGAAGCCGTGGCACAg         >  5:127202/1‑138 (MQ=255)                                                                                                                          cAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCTCCGCCAAAAGCTGAGGAAAAGCCGGCTTACCATGAGAAGCAGCAGGAGGAAGCCGTGGCACAg         >  7:155836/1‑138 (MQ=255)                                                                                                                          cAGCCGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCTCCGCCAAAAGCTGAGGAAAAGCCGGCTTACCATGAGAAGCAGCAGGAGGAAGCCGTGGCACAg         >  4:355213/1‑138 (MQ=255)                                                                                                                             ccGGATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCCCTCCGCCAAAAGCTGAGGAAAAGCCGGCTTACCATGAGAAGCAGCAGGAGGACGCCGTGGCACAGGATg     >  4:28353/1‑139 (MQ=255)                                                                                                                             ccGGATTATGCAGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCTCCGCCAAAAGCTGAGGAAAAGCCGGCTTACCATGAGAAGCAGCAGGAGGACGCCGTGGCACAGGATg     >  5:230017/1‑139 (MQ=255)                                                                                                                                gATTATGCGGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCTCCGCCAAAAGCTGAGGAAAAGCCGGCTTACCATGAGAAGCAGCAGGAGGAAGCCGTGGCACAGGATGTcc  <  4:30408/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                   CTTTTTTATTTCTTCTGATTTTCAGCTTTCTGTCATATAGATAGAATATGACACAATCGGAGGGATGAAGATGTCTAATAAGAAAAGAAACGCGGATAAAAAACCGCTCATGTATATTGTTCAGCCGGATTATGCAGAAACTACAAGGTCATCCATGCAGGAGATTGTCATCAAAAGAAAGGCCGACAAACCCGCTCCGCCAAAAGCTGAGGAAAAGCCGGCTTACCATGAGAAGCAGCAGGAGGAAGCCGTGGCACAGGATGTCC  >  CP006881/1211635‑1211900

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