Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	656,852	T→C	5.8%	T31T (ACT→ACC)	USA300HOU_0595 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	656,852	0	T	C	5.8%	177.6 / 3.1	70	T31T (ACT→ACC)	USA300HOU_0595	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (45/20);  new base C (2/2);  total (48/22)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.99e-01

GTTCTTCAAATAAAAAATAATTTGAATAAGGTGGAAGTGATAAACAGTGCTTTGCGAATCTAGAGTCATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAAATAATACTTTCTGGCCTGTTGCAGGATTTGCAG  >  CP000730/656714‑656991                                                                                                                                           |                                                                                                                                            gTTCTTCAAATAAAAAATAATTTGAATAAGGTGGAAGTGATAAACAGTGCTTTGCGAATCTAGAGTCATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTg                                                                                                                                           >  1:393065/1‑141 (MQ=255) gTTCTTCAAATAAAAAATAATTTGAATAAGGGGGAAGTGATAAACAGTGCTTTGGGAATCTCGAGTCATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACtag                                                                                                                                           >  2:297713/1‑139 (MQ=255)       cAAATAAAAAATAATTTGAATAAGGTGGAAGTGATAAACAGTGCTTTGCGAATCTAGAGTCATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACCTGGATTGc                                                                                                                                     >  1:792758/1‑141 (MQ=255)       cAAATAAAAAATAATTTGAATAAGGTGGAAGTGATAAACAGTGCTTTGCGAATCTAGAGTCATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATCGTTGATTTAGTAAGTACCTGGATTGc                                                                                                                                     >  1:959004/1‑141 (MQ=255)        aaaTAAAAAATAATTTGAATAAGGTGGAAGTGATAAACAGTGCTTTTCGAATCGAGAGTCATTAATCAAAACCCTAAATATAGAAGTATTAAATACAATAGTGAGTATTTCATGGTTGAGTTACTGAGTATTAGGATTGCt                                                                                                                                    >  2:215315/1‑141 (MQ=255)        aaaTAAAAAATAATTTGAATAAGGTGGAAGTGATAAACAGTGCTTTGCGAATCTAGAGTCATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCt                                                                                                                                    >  2:40524/1‑141 (MQ=255)        aaaTAAAAAATAATTTGAATAAGGTGGAAGTGATAAACAGTGCTTTGCGAATCTAGAGTCATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCt                                                                                                                                    >  1:358697/1‑141 (MQ=255)        aaaTAAAAAATAATTTGAATAAGGTGGAAGTGATAAACAGTGCTTTGCGAATCTAGAGTCATTAATCAAAACCCTAAATATAGAATTATTAAATAATATAATGAATATTTAATGGTTGATTTAGTAAGTACTTGGATTgat                                                                                                                                    >  2:565168/1‑139 (MQ=255)         aaTAAAAAATAATTTGAATAAGGTGGAGGTGATAAACAGTGCTTTGCGAATCTAGAGTCATTAATAAAAACCCTAAATATAGAATTATTAAATATAATAATGATTATTTAATTCGTGATTTAGTAAGAAATTGgattgatt                                                                                                                                   >  2:803712/1‑141 (MQ=18)          aTAAAAAATAATTTGAATAAGGTGGAAGTGATAAACAGTGCTTTGCGAATCTAGAGTCATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTa                                                                                                                                  <  1:458375/141‑1 (MQ=255)           tAAAAAATAATTTGAATAAGGTGGAAGTGATAAACAGTGCTTTGCGAATCTAGAGTCATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTAt                                                                                                                                 <  2:357007/141‑1 (MQ=255)              aaaaTAATTTGAATAAGGTGGAAGTGATAAACAGTGCTTTGCGAATCTAGAGTCATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTGAGTAAGTACTTGGATTGCTTAtttt                                                                                                                              >  1:417907/1‑141 (MQ=255)              aaaaTAATTTGAATAAGGTGGAAGTGATAAACAGTGCTTTGCGAATCTAGAGTCATTAATCAAAACACTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTAtttt                                                                                                                              >  2:969845/1‑141 (MQ=255)              aaaaCAATTTGAATGGGGTGGAAGTTGTAAACAGTGCTTTGCGAATCTTGAGACATTAAGCAAATCCCTAAATACAAAATTATTAAATACAATAATGAAGATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTAtttt                                                                                                                              <  2:166867/141‑1 (MQ=21)                 ataatTTGAATAAGGTGGAAGTGATAAACAGTGCTTTGCGAATCTAGAGTCATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTg                                                                                                                           >  1:714573/1‑141 (MQ=255)                      ttGAATAAGGTGGAAGTGATAAACAGTGCTTTGCGAATCTAGAGTCATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCAt                                                                                                                      <  1:198372/141‑1 (MQ=255)                        taaTAAGGTGGAAGAGATAAACAGTGCTTTGCGAATCTAGAGTCATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTA‑TTTTTGCCCATGAt                                                                                                                   <  1:837341/140‑1 (MQ=255)                        gAATAAGGTGGAAGTGATAAACAGTGCTTTGCGAATCTAGAGTCATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGa                                                                                                                    <  2:406221/141‑1 (MQ=255)                        gAATAAGGTGGAAGTGATAAACAGTGCTTTGCGAATCTAGAGTCATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTA‑TTTTTGCCCATGAt                                                                                                                   <  1:837023/141‑1 (MQ=255)                         aagaaGGTGCAAGTGATAAACAGTGCTTTGCGAATCTAGATTCATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGAt                                                                                                                   <  1:901085/138‑1 (MQ=255)                                      tGATAAACAGTGCTTTGCGAATCTAGAGTCATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCAATGATTAATTGGTTTAtt                                                                                                      >  2:428959/1‑141 (MQ=255)                                        aTAAACAGTGCTTTGCGAATCTAGAGTCATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTcc                                                                                                    >  1:160960/1‑141 (MQ=255)                                         tAAACAGTGCTTTGCGAATCTAGAGTCATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTccc                                                                                                   <  1:425534/141‑1 (MQ=255)                                               gTGCTTTGCGAATCTAGAGTCATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAg                                                                                             >  2:944279/1‑141 (MQ=255)                                                   tttGCGAATCTAGAGTCATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTTTTTTTTTTTAATTGGTTTATTCCCATAAAGTAcg                                                                                         >  2:333299/1‑141 (MQ=25)                                                   tttGCGAATCTAGAGTCATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCTAGAAGTAcg                                                                                         >  2:1075880/1‑141 (MQ=255)                                                   tttGCGAATCTAGAGTCATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTAcg                                                                                         <  1:432461/141‑1 (MQ=255)                                                   tttGCGAATCTAGAGTCATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTAcg                                                                                         >  1:75008/1‑141 (MQ=255)                                                   tttGCGAATCTAGAGTCATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTAcg                                                                                         >  1:947649/1‑141 (MQ=255)                                                   tttGCGAATCTAGAGTCATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTAcg                                                                                         >  1:881608/1‑141 (MQ=255)                                                           tCTAGAGTCATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAAtt                                                                                 <  1:1021848/141‑1 (MQ=255)                                                           tCTAGAGTCATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAAtt                                                                                 <  2:34961/141‑1 (MQ=255)                                                                agTCATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATtagtag                                                                            >  1:881407/1‑141 (MQ=255)                                                                agTCATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATtagtag                                                                            >  1:837204/1‑141 (MQ=255)                                                                agTCATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACGTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCGAAAAGTACGCGCACATtggtag                                                                            >  1:178527/1‑141 (MQ=255)                                                                 gTCTTTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAGTATAGGCTCGTGGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACCCGAAAATTAGTaga                                                                           <  1:577624/141‑1 (MQ=21)                                                                   cATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAgaga                                                                         >  1:1060297/1‑141 (MQ=255)                                                                   cATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAgaga                                                                         >  2:891314/1‑141 (MQ=255)                                                                             aaCCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTTTTTTTTTTTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGaa                                                               >  2:113055/1‑141 (MQ=21)                                                                             aaCCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTTCTTTTGTTTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGaa                                                               >  2:366664/1‑141 (MQ=21)                                                                                  tAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAGTTGGTTTATGCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGttt                                                          <  1:428959/141‑1 (MQ=255)                                                                                  tAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGttt                                                          >  1:198932/1‑141 (MQ=255)                                                                                          gAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTg                                                  >  2:145190/1‑141 (MQ=255)                                                                                                 taAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGCTTGCTTATTTTTTGCCAATGATTCATTGGTTTATTCCATGAAAGTACGCGGAAATTAGTATAGGCGAGTTTGAAAGTGTAAatgatataacatc                                           >  2:857417/1‑130 (MQ=255)                                                                                                      aCAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCACACCCGCTa                                      >  2:679168/1‑141 (MQ=255)                                                                                                      aCAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTa                                      >  1:532343/1‑141 (MQ=255)                                                                                                        aataatGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTaaa                                    >  1:1030714/1‑141 (MQ=255)                                                                                                         ataatGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTaaaa                                   >  1:181786/1‑141 (MQ=255)                                                                                                            atGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAaata                                >  1:1026531/1‑141 (MQ=255)                                                                                                              gAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCGTGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAaataat                              >  1:1000446/1‑141 (MQ=255)                                                                                                              gAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAaataat                              >  1:926029/1‑141 (MQ=255)                                                                                                                atatTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAAATAATAc                            >  2:248114/1‑141 (MQ=255)                                                                                                                atatTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAAATAATAc                            >  2:796771/1‑141 (MQ=255)                                                                                                                     tCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAAATAATACTTTCt                       >  1:274588/1‑141 (MQ=255)                                                                                                                         ggTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAAATAATACTTTCTGGcc                   >  1:178939/1‑141 (MQ=255)                                                                                                                          gTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAAATAATACTTTCTGGCCt                  >  2:250417/1‑141 (MQ=255)                                                                                                                            tGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAAATAATACTTTCTAGCCTGt                >  1:654682/1‑141 (MQ=255)                                                                                                                             gATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGCGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAAATAATACTTTCTGGCCTGtt               >  2:664705/1‑141 (MQ=255)                                                                                                                              aTTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAAATAATACTTTCTGGCCTGTTg              <  1:796771/141‑1 (MQ=255)                                                                                                                              aTTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAAATAATACTTTCTGGCCTGTTg              <  1:248114/141‑1 (MQ=255)                                                                                                                              aTTCAGTAAGTACTTGGATTGTTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAATAATTTGAAAGCTTAAATAGTGTCAAAGCCGCTAAAAATAATACTTTCTGGCCTGTTg              <  2:519772/141‑1 (MQ=255)                                                                                                                               tttAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAAATAATACTTTCTGGCCTGTTGc             >  2:745169/1‑141 (MQ=255)                                                                                                                               tttAGGAAGTACTTGAGTTGCTTATTTTTTGACCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAAATAATACTTTCTGGCCTGTTGc             <  2:111335/141‑1 (MQ=255)                                                                                                                                  agtaagtaCCTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAAATAATACTTTCTGGCCTGTTGCAgg          <  2:792758/141‑1 (MQ=255)                                                                                                                                  agtaagtaCCTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAAATAATACTTTCTGGCCTGTTGCAgg          <  2:959004/141‑1 (MQ=255)                                                                                                                                   gtaagtaCTTGGATTGCTTATTTTTTGCCTATTATTAATTGGTTTATTCCCAAAAAGTACGCGACAATTAGTAGAGAAGAATTTGATAGTTTAAATATTGTCAAACCCGGTAAAATTAATACTTTGTGCACTGTTgtctgt         >  2:1068915/1‑136 (MQ=255)                                                                                                                                   gtaagtaCTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAAATAATACTTTCTGGCCTGTTGCAGGa         >  1:915230/1‑141 (MQ=255)                                                                                                                                    taagtaCTTGGATTGCTTATTTTTTGCCCATGATTAATGGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAAATAATACTTTCTGGCCTGTTGCAGGAt        <  1:420676/141‑1 (MQ=255)                                                                                                                                    gaagtaCTTGGATTGCTTATTTTTTGCCCATGATCAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAAATAATACTTTCTGGCCTGTTGCAGGAt        <  1:319468/140‑1 (MQ=255)                                                                                                                                          cTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAAATAATACTTTCTGGCCTGTTGCAGGATTTGCAg  >  1:621150/1‑141 (MQ=255)                                                                                                                                          cTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAAATAATACTTTCTGGCCTGTTGCAGAATTTGCAg  >  2:457673/1‑141 (MQ=255)                                                                                                                                           |                                                                                                                                            GTTCTTCAAATAAAAAATAATTTGAATAAGGTGGAAGTGATAAACAGTGCTTTGCGAATCTAGAGTCATTAATCAAAACCCTAAATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAAATAATACTTTCTGGCCTGTTGCAGGATTTGCAG  >  CP000730/656714‑656991

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: ‑