Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,542,093	G→A	31.4%	V563M (GTG→ATG)	feoB →	ferrous iron transporter protein B and GTP‑binding protein; membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,542,093	0	G	A	31.4%	111.1 / 63.2	84	V563M (GTG→ATG)	feoB	ferrous iron transporter protein B and GTP‑binding protein; membrane protein
Reads supporting (aligned to +/- strand):  ref base G (23/34);  new base A (11/15);  total (34/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

AACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGA  >  NC_000913/3541994‑3542190                                                                                                    |                                                                                                  aaCAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCg                                                                                                   <  1:1655042/100‑1 (MQ=255)  aCAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGt                                                                                                  >  1:1972618/1‑100 (MQ=255)   cAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCTCCCCGGTCTTAAAGCCAATTGGCGTg                                                                                                 <  2:780359/100‑1 (MQ=255)    aGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGc                                                                                                <  1:1434185/100‑1 (MQ=255)      cTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCAt                                                                                              <  1:955255/100‑1 (MQ=255)      cTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCAt                                                                                              >  2:257448/1‑100 (MQ=255)      cTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCa                                                                                               <  2:1099334/99‑1 (MQ=255)      cTTCGCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCATCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATg                                                                                                 <  2:1202843/97‑1 (MQ=255)       ttCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATg                                                                                             >  2:1058391/1‑100 (MQ=255)          tcGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAg                                                                                          <  2:2026398/100‑1 (MQ=255)           cgttGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGa                                                                                         <  1:1829624/97‑1 (MQ=255)           cGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGAGTGATCACCCCGGTCTTCAAGCCAATTCGCGTGCATGAAGa                                                                                         >  2:1293655/1‑100 (MQ=255)             cTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATa                                                                                       >  1:1945552/1‑100 (MQ=255)                 gCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACt                                                                                    >  1:595632/1‑99 (MQ=255)                 gCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACt                                                                                    >  2:895686/1‑99 (MQ=255)                 gCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTg                                                                                   >  2:1619998/1‑100 (MQ=255)                 gCGGGAAAATCGTCGATAACATCAACGACTCGGCGCGGGCGTCCGTCGTCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTg                                                                                   >  2:743918/1‑100 (MQ=255)                   gggAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCCCCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTggc                                                                                 >  1:1031772/1‑100 (MQ=255)                   gggAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATaaat                                                                                    >  1:147068/1‑95 (MQ=255)                   gggAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCA‑CCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTggca                                                                                <  2:1576380/100‑1 (MQ=255)                     gAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTggcag                                                                               >  2:428740/1‑100 (MQ=255)                      aaaaTCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTggcagg                                                                              <  1:1160283/100‑1 (MQ=255)                             tcgATAACATCAACGACTCGGCGCTGGAGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGt                                                                       <  2:705581/100‑1 (MQ=255)                              cgATAACATCAACTATTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGtt                                                                      <  1:187230/100‑1 (MQ=255)                               gATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGTCAGGCAACGGtcg                                                                     <  1:1521751/100‑3 (MQ=255)                               gATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTg                                                                     <  2:542775/100‑1 (MQ=255)                                agaaCATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTgg                                                                    >  1:1516043/3‑100 (MQ=255)                                aTAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTgg                                                                    <  1:1438846/100‑1 (MQ=255)                                 tAACATCAACGGGTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAACTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGc                                                                   <  2:1192509/100‑1 (MQ=255)                                 tAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGc                                                                   <  2:201514/100‑1 (MQ=255)                                  aaCATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGcc                                                                  <  1:1295769/100‑1 (MQ=255)                                  aaCATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGcc                                                                  >  1:1643108/1‑100 (MQ=255)                                    cATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATAACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCt                                                                 >  2:2061602/1‑99 (MQ=255)                                       cAACGACGCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCCTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGttt                                                             <  2:1382579/100‑1 (MQ=255)                                        aaCGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTa                                                            >  2:1440238/1‑100 (MQ=255)                                        aaCGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGttt                                                             <  2:1783927/99‑1 (MQ=255)                                        aaCGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTa                                                            <  1:301166/100‑1 (MQ=255)                                         aCGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTAc                                                           <  1:429413/100‑1 (MQ=255)                                         aCGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTAc                                                           >  2:415388/1‑100 (MQ=255)                                          cGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACa                                                          <  2:1661963/100‑1 (MQ=255)                                             cTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGt                                                       >  2:153367/1‑100 (MQ=255)                                             cTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACATGt                                                       >  1:458336/1‑100 (MQ=255)                                                ggCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGcc                                                    <  2:919134/100‑1 (MQ=255)                                                     tGGCGTCCGTCAGTCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATgg                                                <  1:615321/99‑1 (MQ=255)                                                     tGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGc                                               <  1:1058393/100‑1 (MQ=255)                                                     tGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGc                                               <  2:1182541/100‑1 (MQ=255)                                                        cgtccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGaa                                            >  1:1283479/1‑100 (MQ=255)                                                        cgtccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGCTTGACCTGTTTACAGGTGCCATGGCGaa                                            >  2:2087148/1‑100 (MQ=255)                                                         gtccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGaaa                                           <  2:940995/100‑1 (MQ=255)                                                          tccgtcAGCAGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACAGTTGGCCTGTTTACAGGTGCCATGGCGAaag                                          <  2:1259858/100‑1 (MQ=255)                                                           ccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGaaaaa                                         >  1:1770889/1‑98 (MQ=255)                                                            cgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAaagaa                                        >  2:1312302/1‑100 (MQ=255)                                                             gtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAGGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAaagaag                                       >  2:1074647/1‑100 (MQ=255)                                                               cAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTg                                     <  1:563928/100‑1 (MQ=255)                                                                aGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGCTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTgg                                    <  1:1579081/100‑1 (MQ=255)                                                                      tttatcatcCCTGTCTTCAATCCAATTGTCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTg                               <  1:43566/96‑1 (MQ=255)                                                                      gTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTg                               <  1:1117825/99‑1 (MQ=255)                                                                      gTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTgg                              <  1:1259588/100‑1 (MQ=255)                                                                       tGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTggg                             >  1:1669665/1‑100 (MQ=255)                                                                       tGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTggg                             >  2:26748/1‑100 (MQ=255)                                                                        gATCACCCCGGTCTTCAAGCCAATTGGCGTCCATGAAGATACCTGGCAGTCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGt                                <  2:962045/96‑1 (MQ=255)                                                                         aTCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTa                           <  2:1606210/100‑1 (MQ=255)                                                                          tCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTAc                          <  2:670524/100‑1 (MQ=255)                                                                           cACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACg                         >  2:1477098/1‑100 (MQ=255)                                                                            aCCCCGGTCTTCAAGCCAATTGGCGTGCATGAAAATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGTCATGGCGAAAGAAGTGGTAGTGGGTACGc                        >  2:1349417/1‑100 (MQ=255)                                                                               ccGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCa                     <  1:1312308/100‑1 (MQ=255)                                                                               ccGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCa                     <  2:1283473/100‑1 (MQ=255)                                                                               ccGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCa                     <  2:2249465/100‑1 (MQ=255)                                                                                 ggTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAac                   <  1:740982/100‑1 (MQ=255)                                                                                 ggTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAac                   <  1:2028078/100‑1 (MQ=255)                                                                                     ttCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACAccc               <  1:596470/100‑1 (MQ=255)                                                                                      tCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACAccc               <  1:1487779/99‑1 (MQ=255)                                                                                      tCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCct              <  1:873938/100‑1 (MQ=255)                                                                                       cAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctc             <  1:1525588/100‑1 (MQ=255)                                                                                       cAAGCCAATTGGCGTGCATGAAGAGAACTGGCAGGCACCGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctc             <  2:1514264/100‑1 (MQ=255)                                                                                        aaGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctct            >  1:1137522/1‑100 (MQ=255)                                                                                         aGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTa           >  2:1373188/1‑100 (MQ=255)                                                                                         aGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTa           >  1:1164758/1‑100 (MQ=255)                                                                                         aGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGACGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTa           >  2:1506801/1‑100 (MQ=255)                                                                                          gtcAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTac          <  1:984806/98‑1 (MQ=255)                                                                                            cAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTacac        <  2:1421880/100‑1 (MQ=255)                                                                                               ttGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGc     <  1:428739/100‑1 (MQ=255)                                                                                                tGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCa    <  1:1705369/100‑1 (MQ=255)                                                                                                 ggCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAg   >  1:1237939/1‑100 (MQ=255)                                                                                                  gCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGa  >  1:194955/1‑100 (MQ=255)                                                                                                    |                                                                                                  AACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGA  >  NC_000913/3541994‑3542190

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