Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,194,516	A→G	100%	K551K (AAA→AAG)	yqiK →	PHB family membrane protein, function unknown

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,194,516	0	A	G	100.0%	305.1 / NA	83	K551K (AAA→AAG)	yqiK	PHB family membrane protein, function unknown
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (45/38);  total (45/38)

CGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTG  >  NC_000913/3194419‑3194611                                                                                                  |                                                                                                cGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGc                                                                                               <  1:227715/100‑1 (MQ=255) cGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGc                                                                                               <  2:31000/100‑1 (MQ=255)   cTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCag                                                                                             <  1:2133810/100‑1 (MQ=255)   cTGATTGACGCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCag                                                                                             <  2:1606208/100‑1 (MQ=255)    tGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCaga                                                                                            <  2:968592/100‑1 (MQ=255)     gATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCagag                                                                                           >  2:1330265/1‑100 (MQ=255)     gATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCagag                                                                                           >  1:1310690/1‑100 (MQ=255)     gATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGACTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCagag                                                                                           >  1:418082/1‑100 (MQ=255)         gACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTaaa                                                                                       >  2:938117/1‑100 (MQ=255)         gACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTaaa                                                                                       <  1:757018/100‑1 (MQ=255)          aCTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAAc                                                                                      <  2:126068/100‑1 (MQ=255)           cTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTaaa                                                                                       <  2:2276886/98‑1 (MQ=255)           cTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAAcc                                                                                     >  1:291249/1‑100 (MQ=255)              cTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCctc                                                                                  <  2:1527192/100‑1 (MQ=255)                    cAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAg                                                                            >  1:1559709/1‑100 (MQ=255)                    cAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAg                                                                            <  1:642624/100‑1 (MQ=255)                     aaTGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAgg                                                                           <  1:1302707/100‑1 (MQ=255)                     aaTGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAgg                                                                           <  2:508103/100‑1 (MQ=255)                      aTGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGa                                                                          >  1:1961864/1‑100 (MQ=255)                       tGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACATCTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGaa                                                                         <  1:868762/100‑1 (MQ=255)                         aaaTTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGaata                                                                       <  1:1334629/100‑1 (MQ=255)                           aTTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGaataat                                                                     >  1:1277889/1‑100 (MQ=255)                             tGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATcc                                                                   <  1:1040669/100‑1 (MQ=255)                             tGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATc                                                                    >  2:1020301/1‑99 (MQ=255)                              ggCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATccc                                                                  >  1:1554353/1‑100 (MQ=255)                              ggCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATccc                                                                  >  2:442739/1‑100 (MQ=255)                                cGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTg                                                                >  2:1015326/1‑100 (MQ=255)                                 gTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTgg                                                               >  1:1565790/1‑100 (MQ=255)                                  tttCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGc                                                              <  2:891746/100‑1 (MQ=255)                                   ttCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGcc                                                             <  2:625782/100‑1 (MQ=255)                                    tCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGccc                                                            >  1:1336022/1‑100 (MQ=255)                                     cAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCa                                                           >  2:2116403/1‑100 (MQ=255)                                     cAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCa                                                           >  1:1154295/1‑100 (MQ=255)                                      aGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGGATTATTCCTggaagag                                                          <  1:165191/100‑6 (MQ=255)                                       ggcggcTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCaaa                                                         >  2:1823770/1‑100 (MQ=255)                                          ggcTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCaaaaa                                                       >  1:1064171/1‑99 (MQ=255)                                           gcTCACTGGCGGCCTTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGc                                                     >  1:243958/1‑100 (MQ=255)                                           gcTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGc                                                     >  2:708845/1‑100 (MQ=255)                                            cTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGcc                                                    >  1:933649/1‑100 (MQ=255)                                              cACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCgg                                                  >  1:1223903/1‑100 (MQ=255)                                               aCTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCggg                                                 <  1:127815/100‑1 (MQ=255)                                                 tGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGa                                               <  1:648052/100‑1 (MQ=255)                                                  ggcggcAATGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGa                                               >  1:1236788/1‑99 (MQ=255)                                                     ggcATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAtttt                                           <  1:1513950/100‑1 (MQ=255)                                                      gcATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAttttt                                          <  1:434344/100‑1 (MQ=255)                                                       cATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAtttttt                                         <  1:797184/100‑1 (MQ=255)                                                       cATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAtttttn                                         <  1:43361/100‑2 (MQ=255)                                                        aTTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTta                                        <  1:832923/100‑1 (MQ=255)                                                        aTTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTta                                        <  1:85873/100‑1 (MQ=255)                                                        aTTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTta                                        <  1:938280/100‑1 (MQ=255)                                                           gACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTatat                                     >  1:1148204/1‑100 (MQ=255)                                                            aCTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATc                                    >  2:2187049/1‑100 (MQ=255)                                                             cTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCt                                   >  1:1463001/1‑100 (MQ=255)                                                             cTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCt                                   >  1:439945/1‑100 (MQ=255)                                                             cTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCt                                   >  1:1749919/1‑100 (MQ=255)                                                              ttCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTg                                  <  1:2116421/100‑1 (MQ=255)                                                                  cccTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGtt                              >  2:1773084/1‑100 (MQ=255)                                                                   ccTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTc                             <  1:1625983/100‑1 (MQ=255)                                                                   ccTTAACCTCAACAACTCCCGTAGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGc                                 <  2:225605/96‑1 (MQ=255)                                                                    cTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGCATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTcc                            <  2:1310684/100‑1 (MQ=255)                                                                         ccTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTaa                       >  1:1028544/1‑100 (MQ=255)                                                                         ccTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTaa                       >  2:2202653/1‑100 (MQ=255)                                                                           tCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTaaaa                     >  1:727950/1‑100 (MQ=255)                                                                            caacaaCTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTaaaaa                    >  2:1590791/1‑100 (MQ=255)                                                                            caacaaCTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTaaaaa                    >  2:1305309/1‑100 (MQ=255)                                                                            caacaaCTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTaaaaa                    >  1:34963/1‑100 (MQ=255)                                                                             aacaacTCCCGTCGAAGAAAAGGCAGAGTATACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTaaaaa                    <  2:1582919/99‑1 (MQ=255)                                                                                  cTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCa               >  2:1485216/1‑99 (MQ=255)                                                                                    cccGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCCGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAAt            >  2:468385/1‑100 (MQ=255)                                                                                    cccGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAAt            <  1:819766/100‑1 (MQ=255)                                                                                     ccGTAGAAGAAAAGGCAGAGTAAACCATCTTCCAGGAAGAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATg           <  2:425242/100‑1 (MQ=21)                                                                                       gTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGAt         >  1:1541225/1‑100 (MQ=255)                                                                                       gTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGAt         >  1:1541290/1‑100 (MQ=255)                                                                                       gTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGAt         >  1:1291224/1‑100 (MQ=255)                                                                                         cGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCa       >  2:2253205/1‑100 (MQ=255)                                                                                         cGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCa       <  1:1445457/100‑1 (MQ=255)                                                                                         cGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCa       >  2:1897948/1‑100 (MQ=255)                                                                                         cGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCa       <  1:1705214/100‑1 (MQ=255)                                                                                          gaagaaAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAg      <  2:418079/100‑1 (MQ=255)                                                                                          gaagaaAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAg      >  1:1623929/1‑100 (MQ=255)                                                                                            agaaAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCACATGATAAgg     >  1:2049921/1‑99 (MQ=255)                                                                                            agaaAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgt    >  1:1062880/1‑100 (MQ=255)                                                                                             gaaAAGGCAGAGTAAGCCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgt    <  2:1288334/99‑1 (MQ=255)                                                                                              aaaagGCAGAGTAAACCCTCTTCCCGGAATAACCCCAGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAACGTGCAAATGATCAGgttg  <  1:2239218/100‑1 (MQ=21)                                                                                                  |                                                                                                CGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTG  >  NC_000913/3194419‑3194611

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