Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,458,110	C→A	T578T (ACC→ACA)	gspD →	general secretory pathway component, cryptic

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,458,110	0	C	A	100.0%	236.7 / NA	74	T578T (ACC→ACA)	gspD	general secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (44/30);  total (44/30)

ACTATTCAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACCGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGT  >  NC_000913/3457964‑3458256                                                                                                                                                   |                                                                                                                                                   aCTATTCTAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACagag                                                                                                                                                 <  1:3850/150‑1 (MQ=255)  cTATTCAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAgc                                                                                                                                                >  2:712778/1‑150 (MQ=255)       cAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCt                                                                                                                                           >  2:981031/1‑150 (MQ=255)        aaaaCGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTa                                                                                                                                          <  1:355689/150‑1 (MQ=255)              ccGTGCTGGTCAAAACCGGTGAAACGGTTGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCa                                                                                                                                    <  2:779137/150‑1 (MQ=255)                     ggTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGAt                                                                                                                             <  1:604091/150‑1 (MQ=255)                     ggTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAATGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGAt                                                                                                                             <  1:1024853/150‑1 (MQ=255)                      gTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATg                                                                                                                            <  2:312731/150‑1 (MQ=255)                       tCAATACCGGTGAATCGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCATAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGGGCGCTAAACGCAACCTGATgg                                                                                                                           <  2:943859/150‑1 (MQ=255)                        cAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGt                                                                                                                          <  1:578540/150‑1 (MQ=255)                           aaCCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTAtt                                                                                                                       <  1:942616/150‑1 (MQ=255)                           aaCCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTAtt                                                                                                                       <  1:946781/150‑1 (MQ=255)                             ccGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTa                                                                                                                     >  2:707252/1‑150 (MQ=255)                             ccGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTa                                                                                                                     >  1:986395/1‑150 (MQ=255)                               ggTGAAACGGTGGTCCTGGGCGGATTGCT‑GATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTAtcc                                                                                                                  <  1:970666/150‑1 (MQ=255)                                   aaaCGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGTCGATATTCCTTTAGTGGGGCTACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTAtccgtc                                                                                                               <  2:990228/150‑1 (MQ=255)                                   aaaCGGTGGTCCGGGGCGGAAGGCTGGCAGATGTTTCTACAGAGCAAGTGTCAACGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTAAtccgtc                                                                                                               <  1:168981/150‑1 (MQ=255)                                    aaCGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTAtccgtcc                                                                                                              >  1:17967/1‑150 (MQ=255)                                       ggtggtCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCcgac                                                                                                           <  1:628689/150‑1 (MQ=255)                                          ggtCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGAt                                                                                                        >  2:416851/1‑150 (MQ=255)                                               tGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATcc                                                                                                   >  2:278237/1‑150 (MQ=255)                                                 ggCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGt                                                                                                 >  2:866818/1‑150 (MQ=255)                                                  gCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTg                                                                                                <  2:517667/150‑1 (MQ=255)                                                       ttGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACgat                                                                                           >  1:195583/1‑150 (MQ=255)                                                        tGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACgatg                                                                                          <  1:611275/150‑1 (MQ=255)                                                         gCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACgatga                                                                                         >  2:898622/1‑150 (MQ=255)                                                           tGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGatgatg                                                                                       <  2:17967/150‑1 (MQ=255)                                                            ggATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGt                                                                                      >  2:734209/1‑150 (MQ=255)                                                            ggATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGt                                                                                      <  1:1036384/150‑1 (MQ=255)                                                            ggATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGt                                                                                      >  1:283746/1‑150 (MQ=255)                                                               tgatTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTa                                                                                   <  1:712778/150‑1 (MQ=255)                                                               tgatTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTa                                                                                   >  1:383370/1‑150 (MQ=255)                                                                gatTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTAt                                                                                  >  1:957191/1‑150 (MQ=255)                                                                   ttttCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGc                                                                               >  1:994055/1‑150 (MQ=255)                                                                   ttttCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATAATCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGc                                                                               >  1:323279/1‑150 (MQ=255)                                                                    tttCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCt                                                                              >  1:112693/1‑150 (MQ=255)                                                                       cTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCAc                                                                           >  1:573672/1‑150 (MQ=255)                                                                        tAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACt                                                                          >  1:837527/1‑150 (MQ=255)                                                                         aaaGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTg                                                                         <  1:230217/150‑1 (MQ=255)                                                                           agagCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTc                                                                       >  1:881380/1‑150 (MQ=255)                                                                            gagCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCa                                                                      <  2:195583/150‑1 (MQ=255)                                                                                 aGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAaga                                                                 <  1:337757/150‑1 (MQ=255)                                                                                 aGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGAGCGCTAAAAGCAACCTGATGGTATTTATCCGTCCGACGATTATCAGTGACGATGATGTTTATAGCTAACTGTCaaaata                                                                 >  1:947960/1‑148 (MQ=255)                                                                                   tgtCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAgaga                                                               >  1:129513/1‑150 (MQ=255)                                                                                   tgtCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAgaga                                                               >  1:137856/1‑150 (MQ=255)                                                                                    gtCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGaa                                                              >  1:397203/1‑150 (MQ=255)                                                                                     tCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGaaa                                                             >  1:563454/1‑150 (MQ=255)                                                                                      cAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACTCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGATCACTGTCAAAAGAGAAAt                                                            >  1:657991/1‑150 (MQ=255)                                                                                      cAAAGGTTCCTCTGCTTGGCGATATTCCATTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAAt                                                            >  1:364870/1‑150 (MQ=255)                                                                                       aaaGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATa                                                           <  2:799910/150‑1 (MQ=255)                                                                                             tCCTCTGCTTGGCGAAATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCg                                                     <  2:116176/150‑1 (MQ=255)                                                                                              ccTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGt                                                    >  2:98371/1‑150 (MQ=255)                                                                                                   gCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCg                                               >  1:274645/1‑150 (MQ=255)                                                                                                     ttGGCGATATTAATTTAGTGTGGCAACTCTTCCGCTATACCTCCACAGAGAGCGCTAAACTCAACCTGATGGGATTTATCCGTCCGACGATTATCCGAGCCGATGATGTATATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTc                                             >  2:534591/1‑150 (MQ=255)                                                                                                      tGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGGCGATGATGTTTATCGCTCACTGTCAAAAGAGAAATAACACCGTTACCGTCa                                            >  2:484306/1‑150 (MQ=255)                                                                                                        gCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAgg                                          <  2:325896/150‑1 (MQ=255)                                                                                                          gATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAg                                        >  1:402383/1‑150 (MQ=255)                                                                                                                 cTTTAGTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCaacaac                                 >  1:541291/1‑150 (MQ=255)                                                                                                                      gTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCgg                            <  1:898622/150‑1 (MQ=255)                                                                                                                      gTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCgg                            <  1:308816/150‑1 (MQ=255)                                                                                                                      gTGGGGCAACTCTTCCGCTATACCTCCACAGAGCGCGATAAACGCAACCAGATGGTATATATCAGTCCGACGATTATCCGTGACGATGAAGTTTATCGCTCACTCACAAAAGAGAAATACACACGTTACCGTCAGGAGCAACAACAGAgg                            >  1:171117/1‑150 (MQ=255)                                                                                                                        ggggCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGAt                          >  2:784597/1‑150 (MQ=255)                                                                                                                           gCAACTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGAGGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGa                       >  2:943130/1‑150 (MQ=255)                                                                                                                             aaCTCTTCCGCTATACCTCGACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACg                     <  1:98371/150‑1 (MQ=255)                                                                                                                              aCTCTTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACgg                    >  2:763892/1‑150 (MQ=255)                                                                                                                               ctctTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTGTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACggg                   >  1:158461/1‑150 (MQ=255)                                                                                                                                 ctTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGaa                 <  2:784685/150‑1 (MQ=255)                                                                                                                                 ctTCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGGATTTATCCGTCCGCCGATTAGCCGTGACGATGATGTTTATCGCTCACTGTCAACAGAGAAATCCACCCGTGACCGTCAGGAGCGACAACACCGGATCGCCGGGaa                 >  1:98030/1‑150 (MQ=255)                                                                                                                                   tCCGCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAAt               >  1:729666/1‑150 (MQ=255)                                                                                                                                      gCTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCaa            >  2:959012/1‑150 (MQ=255)                                                                                                                                       cTATACCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGATATCaaa           >  2:281788/1‑150 (MQ=255)                                                                                                                                         ataCCTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCGCACTGGCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAg         <  2:259486/150‑1 (MQ=255)                                                                                                                                            ccTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATTCACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAgcgc      >  2:576795/1‑150 (MQ=255)                                                                                                                                            ccTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGAGGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGACAGGAGCAACACCCGCGGATCGACGGGAAATCAAAAgcgc      >  2:580255/1‑150 (MQ=255)                                                                                                                                             cTCCACAGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCt     <  1:734209/150‑1 (MQ=255)                                                                                                                                                cacaGAGCGCGCTGAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGt  >  2:907910/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                   ACTATTCAAAACGCCGTGCTGGTCAAAACCGGTGAAACGGTGGTCCTGGGCGGATTGCTGGATGATTTTTCTAAAGAGCAAGTGTCAAAGGTTCCTCTGCTTGGCGATATTCCTTTAGTGGGGCAACTCTTCCGCTATACCTCCACCGAGCGCGCTAAACGCAACCTGATGGTATTTATCCGTCCGACGATTATCCGTGACGATGATGTTTATCGCTCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGT  >  NC_000913/3457964‑3458256

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