Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	2,018,834	G→C	11.2%	P29A (CCA→GCA)	USA300HOU_1869 ←	transglycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,018,834	0	G	C	11.2%	166.8 / 10.8	64	P29A (CCA→GCA)	USA300HOU_1869	transglycosylase
Reads supporting (aligned to +/- strand):  ref base G (26/28);  new base C (3/4);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

TTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTAAAACTCACTTAGTATATACCTTGAGTTTACCAGTAC  >  CP000730/2018696‑2018973                                                                                                                                            |                                                                                                                                            ttCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGtt                                                                                                                                           <  1:292279/141‑1 (MQ=255)  tCTTAGTTCATCCACATTATCGCTTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTg                                                                                                                                          <  1:240691/141‑1 (MQ=255)  tCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGttt                                                                                                                                          >  2:536535/1‑140 (MQ=255)  tCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTg                                                                                                                                          <  2:12332/141‑1 (MQ=255)   cTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGCTAATGATTAGAATGGTTAATAATATTTTTAATAGTCT‑TCGTTTAATTTTTTTCTTTTTTTGTGGTTTGCAAATTGGGTGa                                                                                                                                         >  2:371671/1‑141 (MQ=255)   cTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAATTGGTTgt                                                                                                                                         >  1:503001/1‑140 (MQ=255)   cTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACGGGTTGa                                                                                                                                         >  1:91052/1‑141 (MQ=255)      aGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGataa                                                                                                                                      <  2:669544/141‑1 (MQ=255)      aGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGAGATAATGATTAGAATGGTTAATAATATTTTTAATAGTCT‑TCGTTAACTTTTTTTCTTTTTGGGCGTTTTGGCACGTGGTTGataa                                                                                                                                      >  1:230204/1‑141 (MQ=255)          cATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATCATGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGTTGGTTTGCCAACTGGTTGTTAATACg                                                                                                                                  >  1:203313/1‑141 (MQ=255)          cATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACg                                                                                                                                  >  2:132082/1‑141 (MQ=255)             ccACATTATCGCGTGTAGCTAAAAAATACATGATACCAATAAACAATGCGATAACGATTAGAGTGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGTCGGTTTGTCACTTGGTGGCTAATACGtat                                                                                                                               >  1:74138/1‑141 (MQ=255)             ccACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGtat                                                                                                                               >  2:23707/1‑141 (MQ=255)                    aTCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATAAGTATTATAGTg                                                                                                                        >  1:332824/1‑141 (MQ=255)                        cgTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGt                                                                                                                    <  1:357872/141‑1 (MQ=255)                        cgTGTAGATAAAAAATACATGATACAAATTAACAATGCGATGCTGCTTTGAATGCTTAATAATA‑TTTTAATCGTTTGTCGCTTAACCACGCGCTGTGTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGt                                                                                                                    <  1:389519/141‑1 (MQ=255)                                tAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGttt                                                                                                            <  1:196314/141‑1 (MQ=255)                                    aaaTACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTTTTTTTCTTCTATTTGCCAACTGCTTTATAATACGTATTATAGTGAGGTTGGTGTTTCatt                                                                                                        >  2:405223/1‑140 (MQ=255)                                    aaaTACATGATACCAATAAACAATGCGATAATGATTACAATGGTTAATAATATTTTTAATAGAAT‑TCGTTTACTTTTTTTGTTATTCGGCGGTTTTCCAACAGCTTGCTAATACGTCTTATCGTTCGGTTCGTGTGTGata                                                                                                        >  2:223650/1‑141 (MQ=255)                                       taaaTGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGc                                                                                                     <  1:709109/138‑1 (MQ=255)                                         cATGATACCAATATACAATGCGATAATGATTCGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGGTGATAATACGTATTATAGTGAGGTTCGTGTTTCATTTGCTc                                                                                                   <  1:310326/141‑1 (MQ=255)                                         cATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTc                                                                                                   >  1:43962/1‑141 (MQ=255)                                         cATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTc                                                                                                   <  1:151107/141‑1 (MQ=255)                                         cATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGGTTGCCAACTGGTTGATAATACGCATTATAGTGAGGTTCGTGTTTCATTTGCTc                                                                                                   <  1:310653/141‑1 (MQ=255)                                         cATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTGCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTc                                                                                                   <  1:384248/141‑1 (MQ=255)                                          aTGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATGTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCa                                                                                                  >  1:426881/1‑141 (MQ=255)                                            taTACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCaaa                                                                                                <  1:472037/140‑1 (MQ=255)                                             aTACCAAGAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCaaaa                                                                                               <  1:452340/141‑1 (MQ=255)                                                   aTAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCa                                                                                         >  1:354711/1‑141 (MQ=255)                                                   aTAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCATTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCa                                                                                         <  1:550128/141‑1 (MQ=255)                                                   aTAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGGAT‑TCGTTTACTTTTTTTCTTTTTAGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGAGCTTCAGAGGGTCAAAATGTTCa                                                                                         >  2:524106/1‑141 (MQ=255)                                                     aaaCAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTTTTTTTCGTCGGTTTGCCAACTGTTTGATAATATGTATTATTGTGAGGTTCGTGTTTCATTTGCTCAATATGTTCAtt                                                                                       >  2:333627/1‑141 (MQ=255)                                                     aaaCAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGCTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCAtt                                                                                       >  1:404803/1‑141 (MQ=255)                                                         attgagATAATGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCAtttgag                                                                                   <  1:78227/136‑1 (MQ=255)                                                           tGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGCTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCAtttgagtt                                                                                 <  1:331454/141‑1 (MQ=255)                                                           cgCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGCTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCAtttgagtt                                                                                 <  1:41752/140‑1 (MQ=255)                                                              gATAATGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTGGTGTTTCATATGCTGAGACTGTGCAtttgagtttgg                                                                              >  1:69607/1‑140 (MQ=255)                                                              gATAATGATTAGAATGGTTAATAATATTATTAATAGTAA‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTACATATGCTCAAAATGTTCAtttgagtttga                                                                              >  2:739820/1‑141 (MQ=255)                                                               aTAATGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGGGGTTCGGGTTTCATATGCTCAAAATGTTCAtttgagtttgag                                                                             >  1:29401/1‑141 (MQ=255)                                                               aTAATGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCAtttgagtttgag                                                                             >  1:733734/1‑141 (MQ=255)                                                                 aaTGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTCATGCGTCGGTTTTCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTa                                                                           <  1:587446/141‑1 (MQ=255)                                                                        aGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCg                                                                    <  1:401177/141‑1 (MQ=255)                                                                          aaTGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGct                                                                  <  1:112142/141‑1 (MQ=255)                                                                                   taataTTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCTTCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCCAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCAt                                                         >  1:266854/1‑141 (MQ=255)                                                                                     atatTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGc                                                       <  1:302441/141‑1 (MQ=255)                                                                                      tatTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCg                                                      >  1:371559/1‑141 (MQ=255)                                                                                      tatTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCg                                                      >  1:163451/1‑141 (MQ=255)                                                                                       atTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGt                                                     >  2:684295/1‑141 (MQ=255)                                                                                       atTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGGTTGAGTACCTATCGCTTCTTTTCATGCGt                                                     <  1:392214/141‑1 (MQ=255)                                                                                        tttttAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTCAGTTTGAGTACCTATCGCTTCTTTTCATGCAtt                                                    >  2:239009/1‑141 (MQ=255)                                                                                                        tttttCTTTTTGTCTTTTTCGGCGGTTTGTCAACTTGTTGATAATACGTATTATAGTGCGGTTCGTGTTTCATTTGCTCATCATGTTCATTTGAGTTTGAGTATCTTCCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTaa                                     <  2:303432/140‑1 (MQ=255)                                                                                                            aCTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTAAAACt                                 <  1:172816/141‑1 (MQ=255)                                                                                                            aCTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAAAACGTGTTATAGTAAGATTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTAACTATCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTAAAAca                                 >  1:277609/1‑140 (MQ=255)                                                                                                            aCTTGTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTAAAACt                                 >  2:560451/1‑141 (MQ=255)                                                                                                                 ttttctttttcGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTAAAACTCACtt                            <  1:634353/141‑1 (MQ=255)                                                                                                                     ctttttcGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGATCCTTCTTTTAAAACTCACTTAGta                        >  1:660088/1‑141 (MQ=255)                                                                                                                            gggggTTTGGCAACTGTTTGATAATACGTATTATAGTGGGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTTCCTATCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTAAAACTCACTTAGTATATACCt                 <  2:481920/138‑1 (MQ=255)                                                                                                                             gctgTTTGCCAATTGGTTAATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTAAAACTCACTTAGTATATACCtt                <  2:203313/138‑1 (MQ=255)                                                                                                                             gcggTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTAAAACTCACTTAGTATATACCtt                <  1:132082/141‑1 (MQ=255)                                                                                                                              cggTTTGCCAATTGTTTGATAATACGTATTATAGGGAGGTTTGTTTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTAAACCTCACTTAGTATATACCTTg               <  2:282603/141‑1 (MQ=255)                                                                                                                              cggTTTGCCAACTGGTTGATAAGACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTAAAACTCACTTAGTATATACCTTg               <  1:303866/141‑1 (MQ=255)                                                                                                                                 tttGCCAACTGCTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTAAAACTCACTTAGTATATACCTTGAGt            <  2:404803/141‑1 (MQ=255)                                                                                                                                 tctGCCAACGGCTTGATAAAACGGATTAAAGTGAGGTGCGTGTTTCATATGCTCAAAATTTTCATTTGAGTTTGAGCACCTATCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTAAAACTCACTTAGTATATACCTTGAGt            <  1:333627/139‑1 (MQ=255)                                                                                                                                    gCCAACTGGTTTATAATACGTATTATAGTGAGGTTCTTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTAACGCTTCTTTTCATGCGTTTGCGTCTTCTTTTAAAACTCACTTAGTATATACCTTGAGTTTa         <  2:332824/141‑1 (MQ=255)                                                                                                                                           ggTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAACATGTTCGTTTCAATCTGAGTACCTATCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTAAAACTCACTTAGTATATACCTTGAGTTTACCCGTAc  <  1:551433/141‑1 (MQ=255)                                                                                                                                            |                                                                                                                                            TTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTAT‑TCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTAAAACTCACTTAGTATATACCTTGAGTTTACCAGTAC  >  CP000730/2018696‑2018973

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