Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,056,257	C→A	17.1%	intergenic (‑1758/+1909)	sspA ← / ← USA300HOU_0997	S01 family glutamyl endopeptidase/bifunctional N‑acetylmuramoyl‑L‑alanine amidase/mannosyl‑glycoprotein endo‑beta‑N‑acetylglucosaminidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,056,257	0	C	A	17.1%	117.5 / 6.1	63	intergenic (‑1758/+1909)	sspA/USA300HOU_0997	S01 family glutamyl endopeptidase/bifunctional N‑acetylmuramoyl‑L‑alanine amidase/mannosyl‑glycoprotein endo‑beta‑N‑acetylglucosaminidase
Reads supporting (aligned to +/- strand):  ref base C (34/17);  new base A (3/9);  total (37/26)
Fisher's exact test for biased strand distribution p-value = 1.96e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGAAAATTGGCGAATACTTGGTGCTCTTAAATATTTAGAATTAGAATTTAAAGAAAGTTTCATTTGTGACACCTCAAATATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTAATAATATCATAACATACTTTTAAAGCAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTAAATACACTGTAGTAACTACTTTTCAAAATGAATAAAGTGGTTACAGTTAATGTACTTAGACAGTATAAAATTATGAATCTTTAAACTGTG  >  CP000730/1056118‑1056396                                                                                                                                             |                                                                                                                                            ttGAAAATTGGCGAATACTTGGCGCTTCTAACTCTTTAGGAATAGAATTTAAAGAAAGTTTCATGTGTACCACCGTGATTTTTGATCAAATATTGTGTTTAAAT‑TTGGAAGGTGTTAATATCATAACATACTTTTAAAGc                                                                                                                                             <  2:364876/140‑1 (MQ=255)    aaaaTTGGCGAATACTTGGTGCTCTTAAATATTTAGAATTAGAATTTAAAGAAAGTTTCATTTGTGACACCTCAAATATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTAATAATATCATAACATACTTTTAAAGCaaa                                                                                                                                          <  1:186533/140‑1 (MQ=255)     aaaTTGGCGAATACTTGGTGCTCTTAAATATTTAGAATTAGAATTTAAAGAAAGTTTCATTTGTGACACCTCAAATATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTAATAATATCATAACATACTTTTAAAGCaaaa                                                                                                                                         >  2:655559/1‑140 (MQ=255)     aaaTTGGCGAATACTTGGTGCTCTTAAATATTTAGAATTAGAATTTAAAGAAAGTTTCATTTGTGACACCTCAAATATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTAATAATATCATAACATACTTTTAAAGCaaaa                                                                                                                                         >  1:480479/1‑140 (MQ=255)     aaaTTGGCGAATACTTGGTGCTCTTAAATATTTAGAATTAGAATTTAAAGAAAGTTTCATTTGTGACACCTCAAATATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTAATAATATCATAACATACTTTTAAAGCaaaa                                                                                                                                         >  2:640038/1‑140 (MQ=255)     aaaTTGGCGAATACTTGGTGCTCTTAAATATTTAGAATTAGAATTTAAAGAAAGTTTCATTTGTGACACCTCAAATATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTAATAATATCATAAAATACTTTTAAAGCaaaa                                                                                                                                         >  2:175023/1‑140 (MQ=255)            cGAATACTTGGTGCTCTTAAATATTTAGAATTAGAATTTAAAGAAAGTTTCATTTGTGACACCTCAAATATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTAATAATATCATAACATACTTTTAAAGCAAAAAAGGGTa                                                                                                                                  >  1:661974/1‑140 (MQ=255)                  cTTGGTGCTCTTAAATATTTAGAATTAGAATTTCTTCTGAAACACATTTGTGTATACTCATAGATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTAATAATATCTTAACATAATTTTAAATAAAAAAAGAGTAAAATAc                                                                                                                            >  1:364604/1‑140 (MQ=255)                    tGGTGCTCTTAAATATTTAGAATTAGAATTTAAAGAAAGTTTCATTTGTGACACCTCAAATATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTAATAATATCATAACATACTTTTAAAGCAAAAAAGAGTAAATTACtt                                                                                                                          >  1:268402/1‑140 (MQ=255)                    tGGTGCTCTTAAATATTTAGAATTAGAATTTAAAGAAAGTTTCATTTGTGACACCTCAAATATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTAATAATATCATAACATACTTTTAAAGAAAAAAAGGGTAAATTACtt                                                                                                                          >  1:161091/1‑140 (MQ=255)                     ggTGCTCTTAAATATTTAGAATTAGAATTTAAAGAAAGTTTCATTTGTGACACCTCAAATATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTAATAATATCATAACATACTTTTAAAGCAAAAAAGGGTAAATTACTTa                                                                                                                         >  2:634293/1‑140 (MQ=255)                      gTGCTATTAAATATTTTGAATTAGTTTTTAAAGAAAGTTTCATTTGTGACACCTCAAATATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTAATTCTACCATAACATACTTTTAAAGCAAAAAAGGGTAAATTACTTaa                                                                                                                        <  1:334133/140‑1 (MQ=255)                           ctTAAATATTTAGAATTAGAATTTAAAGAAAGTTTCATTTGTGACACCTCAAATATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTCGGGACCGGAAAACATACTTTTTAAGCAAAAAAGGGTAAATTACTTAAGCttt                                                                                                                   >  1:549731/1‑140 (MQ=255)                           ctTAAATATTTAGAATTAGAATTTAAAGAAAGTTTCATTTGTGACACCTCAAATATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTAATAATATCATAACATACTTTTAAAGCAAAAAAGGGTAAATTACTTAAGCttt                                                                                                                   >  2:386019/1‑140 (MQ=255)                           ctTAAATATTTAGAATTAGAATTTAAAGAAAGTTTCATTTGTGACACCTCAAATATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTAATAATATCATAACATACTTTTAAAGAAAAAAAGGGTAAATGACTTAAGAttt                                                                                                                   >  1:336458/1‑140 (MQ=255)                              aaaTATTTAGAATTAGAATTTAAAGAAAGTTTCATTTGTGACACCTCAAATATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTAATAATATCATAACATACTTTTAAAGCAAAAAAAGGTAAATTACTTAAGCTTTAtt                                                                                                                >  1:579120/1‑140 (MQ=255)                              aaaTATTTAGAATTAGAATTTAAAGAAAGTTTCATTTGTGACACATAAAATATAAATCGAA‑AAAGAAAAGCCATATAAGATATTAATAATATCATAACATACTTTTAAAGCAAAAAAGGGTAAATTACTTAAGCTTTAtt                                                                                                                >  2:271811/1‑140 (MQ=255)                               aaTATTTAGAATTAGAATTTAAAGAAAGTTTCATTTGTGACACCTCAAATATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTAATAATATCATAACATACTTTTAAAGCAAAAAAGGGTAAATTACTTAAGCTTTAttt                                                                                                               >  2:294982/1‑140 (MQ=255)                                 tatTTAGAATTAGAATTTAAAGAAAGTTTCATTTGTGACACCTCAAATATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTAATAATATCATAACATACTTTTAAAGCAAAAAAGGGTAAATTACTTAAGCTTTATTTGa                                                                                                             >  1:431443/1‑140 (MQ=255)                                         atttgaatgCAAAGAATGTTTCACTTGTGAAACCGCAAATATCACTCAAATATTTACTAAAAAT‑TTACAAAATAAGAATATCATAACAAACTTTTAAAGCAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGaa                                                                                                     <  2:63650/140‑1 (MQ=255)                                          ttagaattTAAAGAAAGTTTCATTTGTGACACCTCAAATATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTAATAATATCATAACATACTTTTAAAGCAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGaaa                                                                                                    <  1:131724/140‑1 (MQ=255)                                            agaattTAAAGAAAGTTTCATTTGTGACACCTCAAATATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTAATAATATCATAACATACTTTTAAAGCAAAAAAGGGTAAATTACTTAAAATTAATTTGAAAAACGaaaaa                                                                                                  >  2:727085/1‑140 (MQ=255)                                            agaattTAAAGAAAGTTTCATTTGAGACACCTCAAATATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTAATAATATCATAACATACTTTTAAAGCAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGaaaaa                                                                                                  >  1:487146/1‑140 (MQ=255)                                            agaaatTAAAGAAAGTTTCATTTGTGACACCTCAAATATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTAATAATATCATAACATACTTTTAAAGCAAAAAAGGATAAATTACTTAAGCTTTATTTGAAAAACGaaaaa                                                                                                  >  1:517438/1‑140 (MQ=255)                                             gaattTAAAGAAAGTTTCATTTGTGACACCTCAAATATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTAATAATATCATAACATACTTTTAAAGCAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAAt                                                                                                 >  1:487132/1‑140 (MQ=255)                                                    aagaaagTTTCATTTGTGACACCTCAAATATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTAATAATATCATAACATACTTTTAAAGCAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTaaa                                                                                          <  1:730850/140‑1 (MQ=255)                                                      gaaagTTTCATTTGTGACACCTCAAATATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTAATAATATCATAACATACTTTTAAAGCAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTAAATa                                                                                        >  1:732813/1‑140 (MQ=255)                                                         ggTTTCATTTGTGACACCTCAAATATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTAATAATATCATAACATACTTTTAAAGCAAAAAAGGGTAAATTACTTAAGCTTTATTTAAAAAACGAAAAATGTCTAAATacac                                                                                     >  2:316230/2‑140 (MQ=255)                                                            ttCATTTGTGACACCTCAAATATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTAATAATATCATAACATACTTTTAAAGCAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTAAATACACTgt                                                                                  >  1:497861/1‑140 (MQ=255)                                                             tCATTTGTGACACCTCAAATATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTAATACTATCATAAAATACTTTTAAAGCAAAAAAGGGTAAATTACTTAAGATTTATTTGAAAAACGAAATATGTCGAAATACACTgaa                                                                                 >  2:657608/1‑138 (MQ=255)                                                               aTTTGTGACACCTCAAATATATAAGCTGTAGAGGCTAAAAAT‑TTAGAAATTAATAATATCATAACATACTTTTAAAGCAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTAAATACACagtagt                                                                               >  1:754952/1‑140 (MQ=255)                                                               aTTTGTGACACCTCAAATATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTAATAATATCATAACATACTTTTAAAGCAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTAAATACACTgtagt                                                                               >  2:244484/1‑140 (MQ=255)                                                                    tgACACCTCAAATATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTACTAATATCATAACATACTTTTAAAGCAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATATCTAAATACACTGTAGTAacta                                                                          >  2:493080/1‑140 (MQ=255)                                                                    tgACACCTCAAATATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTAATAATATCATAACATACTTTTAAAGCAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTAAATACACTGTAGTAacta                                                                          >  2:579847/1‑140 (MQ=255)                                                                      aaacatCAACTATAAATCAAACATTATCCAAAAAT‑TTATAAATTAATCATATCATAACATAAAAAAAAAAAAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTAAATACACTGTAGTAactact                                                                        <  2:747689/135‑1 (MQ=255)                                                                        acCTCAAATATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTAATAATATCATAACATACTTTTAAAGCAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTAAATACACTGTAGTAACTACttt                                                                      >  1:227173/1‑140 (MQ=255)                                                                        acCTCAAATATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTAATAATATCATAACATACTTTTAAAGCAAAAAAGGGTAAATTACTAAAGCTTTATTTGAAAAACGAAAAATGTCTAAATACACTGTAGTAACTACttt                                                                      >  2:22/1‑140 (MQ=255)                                                                                    aaTCAATTATTGTACAAAAAT‑TTAGAAATTAATACTATCATAACATAAATAAAAAAAAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTAAATACACTGTAGTAACTACTTTTCAAAATGAATa                                                          <  2:449703/140‑1 (MQ=255)                                                                                    aaTCAAATATTGTCTAAAAAT‑TTAGAAATTAATAATATCATAACATACTTTTAAAGCAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATTTCTAAATACACTGTAGTAACTACTTTTCAAAATGAATa                                                          <  1:487840/140‑1 (MQ=255)                                                                                    aaTCAAATATTGTCTAAAAAT‑TTAGAAATTAATAATATCATAACATACTTATAAAGCAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTAAATACACTGTAGTAACTACTTTTCAAAATGAATa                                                          <  1:60904/140‑1 (MQ=255)                                                                                    aaTAAAATATTGTCTAAAAAT‑TTAGAAATTAATAATATCATAACATACTATTAAAGCAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTAAATACACTGTAGTAACTACTTTTCAAAATGAATa                                                          <  2:635565/140‑1 (MQ=255)                                                                                          atatTGTCTAAAAAT‑TTAGAAATTAATAATATCTTAACATAATTTTAAAGCAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTAAATACACTGTAGTAACTACTTTTCAAAATGAATAAAGTgg                                                    <  2:456887/140‑1 (MQ=255)                                                                                          ataatttCTAAAAAT‑TTATAAATTAATAATATCATAACATAAAAAAAAAGAAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTAAATACACTGTAGTAACTACTTTTCAAAATGAATAAAGTgg                                                    <  2:480479/134‑1 (MQ=255)                                                                                               gTCTAAAAAT‑TTAGAACTTAATAATATCATAACATACTTTTAAAGCAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTAAATACACTGTAGTAACTACTTTTCAAAATGAATAAAGTGGTTACa                                               >  1:452414/1‑140 (MQ=255)                                                                                                 cTAAAAAT‑TTAGAAATTAATAATATCATAACATACTTTTAAAGCAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTAAATACACTGTAGTAACTACTTTTCAAAATGAATAAAGTGGTTACAGt                                             <  1:294982/140‑1 (MQ=255)                                                                                                 cTAAAAAT‑ATAGAAAGTTATAATATCATAACATACTTTTAAAGCAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTACCTCCCATGTCATATTAAGGTTTCAAAATGAATAAAGTGGTTACAGt                                             >  2:533291/1‑140 (MQ=255)                                                                                                  tAAAAAT‑TTAGAAATTAATAATATCATAACATAAAATAAAAAAAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTAAATACACTGTAGTAACTACTTTTCAAAATGAATAAAGTGGTTACAGtt                                            <  2:431443/140‑1 (MQ=255)                                                                                                       aT‑TTAGAAATTAATAATATCATAACATACTTTTAAAGCAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTAAATACACTGTAGTAACTACTTTTCAAAATGAATAAAGTGGTTACAGTTAATGt                                       <  2:497861/140‑1 (MQ=255)                                                                                                        acttACAAATTAATAATATCCTAACATACAAAAAAAGAAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTAAATACACTGTAGTAACTACTTTTCAAAATGAATAAAGTGGTTACAGTTAATGt                                       <  2:713/138‑1 (MQ=255)                                                                                                           taaaaaTTAATAATATCATAACATAAAAAAAAAAAAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTACATACACTGTAGTAACTACTTTTCAAAATGAATAAAGTGGTTACAGTTAATGTACt                                    <  2:654304/137‑1 (MQ=255)                                                                                                            aGAAATTAATAATATTATAACATACTTTTAAAGCAAAAGACGGTAACATACGAAATTTTTATTTGAAAAACACAATATGCCTCAATACACCGTAGATACTACTTTAAAAAAGATATAACGGGGTTACAGGTGATCTAcat                                   >  2:123023/1‑138 (MQ=255)                                                                                                            aGAAATTAATAATATCATAACATACTTTTAAAGCAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTAAATACACTGTAGTAACTACTTTTCAAAATGAATAAAGTGGTTACAGTTAATGTACtt                                   >  1:600259/1‑140 (MQ=255)                                                                                                            aGAAATTAATAATATCATAACATACTTTTAAAGCAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTAAATACACTGTAGTAACTACTTTTCAAAATGAATAAAGTGGTTACAGTTAATGTACtt                                   >  2:424511/1‑140 (MQ=255)                                                                                                            aGAAATTAATAATATCATAACATACTTTTAAAGCAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTAAATACACTGTAGTAACTACTTTTCAAAATGAATAAAGTGGTTACAGTTAATATACtt                                   >  1:219521/1‑140 (MQ=255)                                                                                                                 ttAATAATATTATAACATACTTTTAAAGCAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTAAATACACTGTAGTAACTACTTTTCAAAATGAATAAAGTGGTTACAGTTAATGTACTTAGACa                              <  1:296426/140‑1 (MQ=255)                                                                                                                 ttAATAATATTATAACATAATAAAAAAAAAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTAAATACACTGTAGTAACTACTTTTCAAAATGAATAAAGTGGTTACAGTTAATGTACTTAGACa                              <  2:504174/140‑1 (MQ=255)                                                                                                                 ttAATAATATCATAACATACTTTTAAAGCAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTAAATACACTGTAGTAACTACTTTTCAAAATGAATAAAGTGGTTACAGTTAATGTACTTAGACa                              <  1:725643/140‑1 (MQ=255)                                                                                                                           cataacataaaaaaaaaGCAAAATAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTAAATACACTGTAGTAACTACTTTTCGAAATGAATAAAGTGGTTACAGTTAATGTACTTAGACAGTATAAAAtt                    <  2:475190/126‑1 (MQ=255)                                                                                                                            aTAACATACTTTTAAAGAAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTAAATACACTGTAGTAACTACTTTTCAAAATGAATAAAGTGGTTACAGTTAATGTACTTAGACAGTATAAAATTa                   <  1:193008/140‑1 (MQ=255)                                                                                                                               aCATACTTTTAAAGCAAAAAAGTGTAATTTACTTAATTTTTATTTGAAAAACGAAAAATGTCTAAATACACTGTAGTAACTACTTTTCAAAATGAATAAAGTGGTTACAGTTAATGTACTTAGACAGTATAAAATTATGa                <  1:453958/140‑1 (MQ=255)                                                                                                                                       ttAAAGAAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTAAATACACTGTAGTAACTACTTTTCAAAATGAATAAAGTGGTTACAGTTAATGTACTTAGACAGTATAAAATTATGAATCTTTaa        <  2:511531/140‑1 (MQ=255)                                                                                                                                       ggaaaGCAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTAAATACACTGTAGTAACTACTTTTCAAAATGAATAAAGTGGTTACAGTTAATGTACTTAGACAGTATAAAATTATGAATCTTTaa        <  2:185610/138‑1 (MQ=255)                                                                                                                                       aaaaaaCAAAAAAAGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTAAATACACTGTAGTAACTACTTTTCAAAATGAATAAAGTGGTTACAGTTAATGTACTTAGACAGTATAAAATTATGAATCTTTaa        <  2:78055/138‑1 (MQ=255)                                                                                                                                          aaGCAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTAAATACACTGTATTAACTACTTTTCAAAATGAATAAAGTGGTTACAGTTAATGTACTTAGACAGTATAAAATTATGAATCTTTAAACt     >  1:300175/1‑140 (MQ=255)                                                                                                                                             cAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTAAATACACTGTAGTAACTACTTTTCAAAATGAATAAAGTGGTTACAGTTAATGTACTTAGACAGTATAAAATTATGAATATTTAAACagtg  >  1:685857/1‑140 (MQ=255)                                                                                                                                             |                                                                                                                                            TTGAAAATTGGCGAATACTTGGTGCTCTTAAATATTTAGAATTAGAATTTAAAGAAAGTTTCATTTGTGACACCTCAAATATAAATCAAATATTGTCTAAAAAT‑TTAGAAATTAATAATATCATAACATACTTTTAAAGCAAAAAAGGGTAAATTACTTAAGCTTTATTTGAAAAACGAAAAATGTCTAAATACACTGTAGTAACTACTTTTCAAAATGAATAAAGTGGTTACAGTTAATGTACTTAGACAGTATAAAATTATGAATCTTTAAACTGTG  >  CP000730/1056118‑1056396

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