Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	2,286,977	A→G	6.0%	I248T (ATC→ACC)	glmM ←	phosphoglucosamine mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,286,977	0	A	G	6.0%	199.4 / 3.3	67	I248T (ATC→ACC)	glmM	phosphoglucosamine mutase
Reads supporting (aligned to +/- strand):  ref base A (25/38);  new base G (2/2);  total (27/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.34e-01

CTTTGTAAAAACCTAAATTACTCATAACAGTAGAAACAATCATGTCATTATTCAATTCTTGATTTTTATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAATTTTTCAGGATGTGTAGAGCCACATTTCTCATTGATATTATATCCATCAGGACTACATCCAATTGTTTCAGTA  >  CP000730/2286837‑2287117                                                                                                                                             |                                                                                                                                             cTTTGTAAAAACCTAAATTACTCATAACAGTAGAAACAATCATGTCATTATTCAATTCTTGATTTTTATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGa                                                                                                                                              >  1:102959/1‑141 (MQ=255)   ttGTAAAAACCTAAATTACTCATAACAGTAGAAACAATCATGTCATTATTCAATTCTTGATTTTTATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGAtt                                                                                                                                            >  1:483945/1‑141 (MQ=255)        aaaaCCTACATGACTCATAGCAGTAGAAACAATCATGTCATTATTCAATTCTTGATTTTTATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGtc                                                                                                                                       <  2:570333/141‑1 (MQ=255)           aCCTAAATTACTCATAACAGTAGAAACAATCATGTCATTATTCAATTCTTGATTTTTATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTcc                                                                                                                                    <  1:144586/141‑1 (MQ=255)               aaaTTACTCATAACAGTAGAAACAATCATGTCATTATTCAATTCTTGATTTTTATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGGCCATTCTCATCTACTGCTATGATTCTGTCTCCATCg                                                                                                                                <  2:9850/141‑1 (MQ=255)                aaTTACTCATAACAGTAGAAACAATCATGTCATTATTCAATTCTTGATTTTTATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGc                                                                                                                               >  2:512441/1‑141 (MQ=255)                aaTTACTCATAACAGTAGAAACAATCATGTCATTATTCAATTCTTGATTTTTATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGc                                                                                                                               >  2:534890/1‑141 (MQ=255)                     cTCATAACAGTAGAAACAATCATGTCATTATTCAATTCTTGATTTTTATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCa                                                                                                                          <  1:234415/141‑1 (MQ=255)                      tCATAACAGTAGAAACAATCATGTCATTATTCAATTCTTGATTTTTATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCaa                                                                                                                         <  2:483945/141‑1 (MQ=255)                      tCATAACAGTAGAAACAATCATGTCATTATTCAATTCTTGATTTTTATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCaa                                                                                                                         >  2:254275/1‑141 (MQ=255)                               tAGAAACAATCATGTCATTATTCAATTCTTGATTTTTATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAAcc                                                                                                                <  2:297047/141‑1 (MQ=255)                               tAGAAACAATCATGTCATTATTCAATTCTTGATTTTTATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAAcc                                                                                                                >  1:265398/1‑141 (MQ=255)                                 gAAACAATCATGTCATTATTCAATTCTTGATTTTTATTTTTTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCa                                                                                                              >  2:429320/1‑141 (MQ=255)                                 gAAACAATCATGTCATTATTCAATTCTTGATTTTTATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCa                                                                                                              >  1:222297/1‑141 (MQ=255)                                  aaaCAATCATGTCATTATTCAATTCTTGATTTTTATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCaa                                                                                                             <  1:268600/141‑1 (MQ=255)                                  aaaCAATCATGTCATTATTCAATTCTTGATTTTTATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCaa                                                                                                             <  1:648374/141‑1 (MQ=255)                                  aaaCAATCATGTCATTATTCAATTCTTGATTTTTATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCaa                                                                                                             <  1:509779/141‑1 (MQ=255)                                      aaTCATGTCATTATTCAATTCTTGATTTTTATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATc                                                                                                         <  2:449008/141‑1 (MQ=255)                                       aTCATGTCATTATTCAATTCTTGATTTTTATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCttattccagtaat                                                                                                        >  2:559684/1‑129 (MQ=255)                                        tCATGTCATTATTCAATTCTTGATTTTTATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCAc                                                                                                       >  2:653071/1‑141 (MQ=255)                                         cATGTCATTATTCAATTCTTGATTTTTATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACt                                                                                                      <  1:162267/141‑1 (MQ=255)                                          aTGTCATTATTCAATTCTTGATTTTTATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCACAATCACtt                                                                                                     >  1:185328/1‑141 (MQ=255)                                            gTCATTATTCAATTCTTGATTTTTATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTc                                                                                                   >  1:38904/1‑141 (MQ=255)                                                 tattCAATTCTTGATTTTTATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGAGTTTACCATTCGCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGttt                                                                                              >  1:172061/1‑141 (MQ=255)                                                  attCAATTCTTGATTTTTATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTc                                                                                             <  1:142445/141‑1 (MQ=255)                                                   ttCAATTCTTGATTTTTATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCa                                                                                            <  1:677512/141‑1 (MQ=255)                                                   ttCAATTCTTGATTTTTATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCa                                                                                            >  1:286873/1‑141 (MQ=255)                                                     cAATTCTTGATTTTTATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAac                                                                                          <  2:256929/141‑1 (MQ=255)                                                     cAATTCTTGATTTTTATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAac                                                                                          >  1:520681/1‑141 (MQ=255)                                                      aaTTCTTGATTTTTATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAact                                                                                         <  2:267087/141‑1 (MQ=255)                                                            tGATTTTTATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACtttt                                                                                   >  1:117755/1‑141 (MQ=255)                                                               tttttATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCa                                                                                >  2:403341/1‑141 (MQ=255)                                                                ttttATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAg                                                                               >  1:362322/1‑141 (MQ=255)                                                                 tttATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGc                                                                              <  2:363295/141‑1 (MQ=255)                                                                      gCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAAtt                                                                         >  2:353322/1‑141 (MQ=255)                                                                       cATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAAttt                                                                        <  2:264093/141‑1 (MQ=255)                                                                        aTTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAAtttt                                                                       <  1:533054/141‑1 (MQ=255)                                                                               gACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCGTCTGCTATGATTCGGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAATTTTTCAGGAt                                                                <  1:599204/141‑1 (MQ=255)                                                                                 ccAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAATTTTTCAGGAtgt                                                              <  2:639186/141‑1 (MQ=255)                                                                                 ccAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAATTTTTCAGGAtgt                                                              >  1:123097/1‑141 (MQ=255)                                                                                         aaaCATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAATTTTTCAGGATGTGTAGAGcc                                                      <  1:382447/141‑1 (MQ=255)                                                                                         aaaCATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTATGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAATTTTTCAGGATGTGTAGAGcc                                                      <  1:219722/141‑1 (MQ=255)                                                                                          aaCATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAATTTTTCAGGATGTGTAGAGCca                                                     <  1:166931/141‑1 (MQ=255)                                                                                            cATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAATTTTTCAGGATGTGTAGAGCcaca                                                   <  1:660248/141‑1 (MQ=255)                                                                                            cATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAATTTTTCAGGATGTGTAGAGCcaca                                                   <  1:653071/141‑1 (MQ=255)                                                                                               aaTTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAATTTTTCAGGATGTGTAGAGCCACAttt                                                >  2:549637/1‑141 (MQ=255)                                                                                               aaTTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAATTTTTCAGGATGTGTAGAGCCACAttt                                                <  1:353322/141‑1 (MQ=255)                                                                                               aaTTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAATTTTTCAGGATGTGTAGAGCCACAttt                                                <  1:338108/141‑1 (MQ=255)                                                                                                 tttGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAATTTTTCAGGATGTGTAGAGCCACATTtct                                              <  1:443339/141‑1 (MQ=255)                                                                                                    ggTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAATTTTTCAGGATGTGTAGAGCCACATTTCTCAt                                           <  1:2172/141‑1 (MQ=255)                                                                                                     gTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAATTTTTCAGGATGTGTAGAGCCACATTTCTCAtt                                          >  1:210175/1‑141 (MQ=255)                                                                                                        aCCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAATTTTTCAGGATGTGTAGAGCCACATTTCTCATTGat                                       <  2:321275/141‑1 (MQ=255)                                                                                                         ccGTCAACGATTTGACCATTCTCATCTACTGCTATGGTTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAATTTTTCAGGATGTGTAGAGCCACATTTCTCATTGata                                      >  1:357994/1‑141 (MQ=255)                                                                                                         ccGTCAACGATTTGACCATTCTCATCTACTGCTATGGTTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAATTTTTCAGGATGTGTAGAGCCACATTTCTCATTGata                                      >  1:357988/1‑141 (MQ=255)                                                                                                         ccGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAATTTTTCAGGATGTGTAGAGCCACATTTCTCATTGata                                      <  2:362322/141‑1 (MQ=255)                                                                                                          cGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAATTTTTCAGGATGTGTAGAGCCACATTTCTCATTGatat                                     >  2:352140/1‑141 (MQ=255)                                                                                                           gTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGTTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAATTTTTCAGGATGTGTAGAGCCACATTTCTCATTGAtatt                                    >  2:393955/1‑141 (MQ=255)                                                                                                             cAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCGGCTAATTTTTCAGGATGTGTAGAGCCACATTTCTCATTGAtgttat                                  <  2:293866/141‑1 (MQ=255)                                                                                                                   tttGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAATTTTTCAGGATGTGTAGAGCCACATTTCTCATTGATATTATATCCAt                            <  1:52687/141‑1 (MQ=255)                                                                                                                               cATCTACTGCTATGGTTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAATTTTTCAGGATGTGTAGAGCCACATTTCTCATTGATATTATATCCATCAGGACTACATc                <  2:357994/141‑1 (MQ=255)                                                                                                                               cATCTACTGCTATGGTTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAATTTTTCAGGATGTGTAGAGCCACATTTCTCATTGATATTATATCCATCAGGACTACATc                <  2:357988/141‑1 (MQ=255)                                                                                                                               cATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAATTTTTCAGGATGTGTAGAGCCACATTTCTCATTGATATTATATCCATCAGGACTACATc                >  1:365709/1‑141 (MQ=255)                                                                                                                               cATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAATTTTTCAGGATGTGTAGAGCCACATTTCTCATTGATATTATATCCATCAGGACTACATc                <  2:222233/141‑1 (MQ=255)                                                                                                                               cATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCCAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAATTTTTCAGGATGTGTAGAGCCACATTTCTCATTGATATTATATCCATCAGGACTACATc                <  2:342075/141‑1 (MQ=255)                                                                                                                                 tCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAATTTTTCAGGATGTGTAGAGCCACATTTCTCATTGATATTATATCCATCAGGACTACATCCa              <  1:536032/141‑1 (MQ=255)                                                                                                                                 tCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAATTTTTCAGGATGTGTAGAGCCACATTTCTCATTGATATTATATCCATCAGGACTACATCCa              <  1:577318/141‑1 (MQ=255)                                                                                                                                     cTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAATTTTTCAGGATGTGTAGAGCCACATTTCTCATTGATATTATATCCATCAGGACTACATCCAATTg          <  2:12978/141‑1 (MQ=255)                                                                                                                                          aTGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAATTTTTCAGGATGTGTAGAGCCACATTTCTCATTGATATTATATCCATCAGGACTACATCCAATTGTTTCa     >  1:334145/1‑141 (MQ=255)                                                                                                                                             aTTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAATTTTTCAGGATGTGTAGAGCCACATTTCTCATTGATATTATATCCATCAGGACTACATCCAATTGTTTCAGTa  >  2:247021/1‑141 (MQ=255)                                                                                                                                             |                                                                                                                                             CTTTGTAAAAACCTAAATTACTCATAACAGTAGAAACAATCATGTCATTATTCAATTCTTGATTTTTATGCATTTCTTGACCAATAATAAACATAATTTGGTCACCGTCAACGATTTGACCATTCTCATCTACTGCTATGATTCTGTCTCCATCGCCGTCAAATGCTAACCCAAAATCACTTTCAGTTTCAACTACTTTTTCAGCTAATTTTTCAGGATGTGTAGAGCCACATTTCTCATTGATATTATATCCATCAGGACTACATCCAATTGTTTCAGTA  >  CP000730/2286837‑2287117

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