Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	803,952	C→T	6.3%	S150S (AGC→AGT)	ybhJ →	aconitase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	803,952	0	C	T	6.3%	195.3 / 2.3	65	S150S (AGC→AGT)	ybhJ	aconitase family protein
Reads supporting (aligned to +/- strand):  ref base C (30/30);  new base T (2/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

CTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTG  >  NC_000913/803807‑804096                                                                                                                                                  |                                                                                                                                                 cTATTAACGGTGATGTCCTTGTTTTTGGTTTATCGGCGGCCCAGCGTTTTGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCTTGAGATGATGGCAGGCTGCGGCAAAATGATCCTCGGGTCAGACAGCcaca                                                                                                                                              <  1:1521010/150‑1 (MQ=255)   aTTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACAcc                                                                                                                                            >  1:1007820/1‑150 (MQ=255)      aaCGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGt                                                                                                                                         >  1:1675905/1‑150 (MQ=255)       aCGGAGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGCGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGtt                                                                                                                                        <  2:1140249/150‑1 (MQ=255)           tgatgaCCATGTTTTTTGTTTATCGTCGGCCCAGCTTTATTGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGAGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACgg                                                                                                                                    <  1:608588/150‑1 (MQ=255)           tgatgaCCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACgg                                                                                                                                    <  2:223925/150‑1 (MQ=255)              tgaCCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGc                                                                                                                                 <  2:1619891/150‑1 (MQ=255)               gaCCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCa                                                                                                                                <  2:1256115/150‑1 (MQ=255)                aCCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCAt                                                                                                                               <  2:930449/150‑1 (MQ=255)                  cATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGTCACACCCGTTACGGTGCATTa                                                                                                                             >  1:694318/1‑150 (MQ=255)                  cATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGTCACACCCGTTACGGTGCATTa                                                                                                                             >  1:225527/1‑150 (MQ=255)                     gTTTTTGGTTTATCGGCGGCCCAGCGTTTTGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAggg                                                                                                                          <  1:519898/150‑1 (MQ=255)                        tttGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACa                                                                                                                       >  1:1585738/1‑150 (MQ=255)                           ggTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGCCAGCCACCCCCGTTACGGTGCATTAGGGACAATg                                                                                                                    >  2:1002504/1‑150 (MQ=255)                            gTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATgg                                                                                                                   <  2:1168406/150‑1 (MQ=255)                             tgtATCGGAGGCCCAGCGTTATGGAGGTATTTTTGTGCCTCCGCATAGTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGAGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGc                                                                                                                  <  1:137175/148‑1 (MQ=255)                                       gcCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGAGGCAAAATGATCCTCGGGTCTGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAg                                                                                                        <  2:1036802/150‑1 (MQ=255)                                         ccAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAggg                                                                                                      <  1:1794420/150‑1 (MQ=255)                                         ccAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAggg                                                                                                      >  2:1459903/1‑150 (MQ=255)                                         ccAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAggg                                                                                                      >  2:1824600/1‑150 (MQ=255)                                                tATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGt                                                                                               <  1:1207381/150‑1 (MQ=255)                                                 aTGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGACAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTg                                                                                              <  2:1355397/150‑1 (MQ=255)                                                  tGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGa                                                                                             >  2:1130180/1‑150 (MQ=255)                                                   ggCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGTTGATGGAAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACCCCCGTTACGGTGCATTAGGGCCAATGGCAGGCGGTGAGGGCGGCGGGGAg                                                                                            >  2:293079/1‑150 (MQ=255)                                                   ggCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAg                                                                                            <  1:1250295/150‑1 (MQ=255)                                                       gTATTTTTGTGCCGCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTgg                                                                                        <  2:906672/150‑1 (MQ=255)                                                             ttGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAAc                                                                                  <  1:190073/150‑1 (MQ=255)                                                              tgcgCCTCCGCATATAGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACa                                                                                 <  1:1824261/147‑1 (MQ=255)                                                               gtgCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCTGGGCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAg                                                                                <  1:237995/150‑1 (MQ=255)                                                               gtgCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAg                                                                                >  2:1692258/1‑150 (MQ=255)                                                                  ccTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAgctg                                                                             >  2:734993/1‑150 (MQ=255)                                                                   cTCCGCATATTGCGGTCATCCATCAATATATGCGTCAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGTCAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAgctgc                                                                            <  2:1320292/150‑1 (MQ=255)                                                                    tCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACTGTGCATTAGGGACAAGGGCAGTCGGTGAGGGCGGCGGTGAGTGGGTAACACAgaggct                                                                           >  2:170791/1‑145 (MQ=255)                                                                    tCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAgctgct                                                                           <  1:680469/150‑1 (MQ=255)                                                                    tCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGGGCATTATGGACAATGGCAGTCGGTGAGGGCGGCGTGGAGTTGGTAAAACAgctgaa                                                                           >  2:165205/1‑148 (MQ=255)                                                                       gCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTGACGGTGCATTAGGGACAAAGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTaa                                                                        >  1:136050/1‑150 (MQ=255)                                                                         atatTGCGGTCATCCATCAATATATGCGGGCGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATg                                                                      <  1:1292690/150‑1 (MQ=255)                                                                              gCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACAcc                                                                 >  2:419228/1‑150 (MQ=255)                                                                              gCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACAcc                                                                 >  2:1369701/1‑150 (MQ=255)                                                                               cGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCt                                                                >  2:1511607/1‑150 (MQ=255)                                                                                     tccatcAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGata                                                          <  2:830765/150‑1 (MQ=255)                                                                                       catcAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATc                                                        >  2:1705928/1‑150 (MQ=255)                                                                                          cAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGTCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGAc                                                     <  2:225527/150‑1 (MQ=255)                                                                                          cAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGTCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGAc                                                     <  2:694318/150‑1 (MQ=255)                                                                                               tatGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATcc                                                >  1:5730/1‑150 (MQ=255)                                                                                                 tGCGTGAGATGATGGCAGGCGGCGGCAAACTGATCCTCGGGTCAGACAGCCACACCCGGTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCgg                                              <  2:1416376/150‑1 (MQ=255)                                                                                                   cGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAACACAGCTGCTTAATGCCACCTGGGATATCGACTATCCGGGc                                            >  2:1221005/1‑150 (MQ=255)                                                                                                   cGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGGCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGc                                            <  2:1406761/150‑1 (MQ=255)                                                                                                      gagaTGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTg                                         >  2:785380/1‑150 (MQ=255)                                                                                                        gatgatGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGt                                       <  1:1726458/150‑1 (MQ=255)                                                                                                            atgGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCg                                   >  2:1745282/1‑150 (MQ=255)                                                                                                                 aGGCGGCGTCAAAATGATCCTCGGGGCAGACAGCCACACGCGTTACGGTGCATTATGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCa                              <  1:192966/150‑1 (MQ=255)                                                                                                                       ctgcAAAATGATCCTCGGGTCAGACAGCCACACACGTTAAGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTAGGTAAAACAGCTGCTTAATGACACCTGGTATATCGACTATCCTTGCGTGGTTGCGGTGCATCTGAc                        >  1:410302/3‑150 (MQ=255)                                                                                                                        ggcAAAATGATCCTCGGTTCAGACATGCACACTCGTTACGGTGCATTAGGGACAATGGCAGTCGGGGAGGGGGGCGGTGTGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGAcc                       <  1:582782/150‑1 (MQ=255)                                                                                                                         gcAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCg                      >  2:1366008/1‑150 (MQ=255)                                                                                                                           aaaaTGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGCCACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGCCCGGa                    >  2:99516/1‑150 (MQ=255)                                                                                                                             aaTGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGaaa                  >  2:1693267/1‑150 (MQ=255)                                                                                                                              aTGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACTGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGaaaa                 >  1:95998/1‑150 (MQ=255)                                                                                                                              aTGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATACCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGaaaa                 >  2:1366871/1‑150 (MQ=255)                                                                                                                               tGATCCTCGGGTCAGACAGCCTCACCCGTAACGGTGCAGTAGGGACAATGGCAGGCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAAc                <  1:476470/150‑1 (MQ=255)                                                                                                                                  tCCTCGGGTCAGCCAGCCACCCCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTg             <  2:629202/150‑1 (MQ=255)                                                                                                                                  tCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCTGGGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTg             <  2:1384480/150‑1 (MQ=255)                                                                                                                                   ccTCGGGTCAGACAGCCACACCCGTTACGGAGCATTAGGGACAATGGCAGTCGGAGAGGGCGGCGCTGAGTAGGTAAAACAGCTGCTAAATGACACCTCCCAAATCGACTATCAGGGACTGGATGCGGTCCATCTGACAGGAAAACCTgc            >  2:508496/1‑150 (MQ=255)                                                                                                                                           cAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGGGAGGGCGGCGGAGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTAtg    >  2:1286273/1‑150 (MQ=255)                                                                                                                                             gACAGCCACACCCGTTACGGTGCACTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTAtgtg  >  1:1201294/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                 CTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTG  >  NC_000913/803807‑804096

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