Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,635,808	C→T	5.2%	A14T (GCA→ACA)	stfQ ←	Qin prophage; putative side tail fibre assembly protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,635,808	0	C	T	5.2%	251.0 / 2.5	76	A14T (GCA→ACA)	stfQ	Qin prophage; putative side tail fibre assembly protein
Reads supporting (aligned to +/- strand):  ref base C (36/36);  new base T (2/2);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

AATGCCGGAGCATGAGCACCGCTTGTACCACTCCAGCCGATGAGTAACTCACCTTCGCCAACGGCTGTCATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGG  >  NC_000913/1635663‑1635947                                                                                                                                                  |                                                                                                                                            aaTGCCGGAGCATGAGCACCGCTTGTACCACTCCAGCCGATGAGTAACTCACCTTCGCCAACGGCTGTCATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACctgct                                                                                                                                            <  7:240104/147‑1 (MQ=255)     ccGGAGCATGAGCACCGCTTGTACCACTCCAGCCGATGAGTAACTCTCCTTCGCCAACGGCTGTCATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTAATTCACCTCGTATAACTCAAGACCTGCTGccc                                                                                                                                        <  8:283763/147‑1 (MQ=255)         aGCATGAGCACCGCTTGTACCACTCCAGCCGATGAGTAACTCACCTTCGCCAACGGCTGTCATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGc                                                                                                                                   >  4:343505/1‑148 (MQ=255)             tGAGCACCGCTTGTACCACTCCAGCCGATGAGTAACTCACCTTCGCCAACGGCTGTCATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTAt                                                                                                                               <  8:122172/148‑1 (MQ=255)              gAGCACCGCTTGTACCACTCCAGCCGATGAGTAACTCACCTTCGCCAACGGCTGTCATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTa                                                                                                                                >  5:304829/1‑146 (MQ=255)              gAGCACCGCTTGTACCACTCCAGCCGATGAGTAACTCACCTTCGCCAACGGCTGTCATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTAtt                                                                                                                              >  3:58832/1‑148 (MQ=255)                 cACCGCTTGTACCACTCCAGCCGATGAGTAACTCACCTTCGCCAACGGCTGTCATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATtgtc                                                                                                                           <  7:311693/148‑1 (MQ=255)                   ccGCTTGTACCACTCCAGCCGATGAGTAACTCACCTTCGCCAACGGCTGTCATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATtgtctg                                                                                                                         <  8:69787/148‑1 (MQ=255)                              aCTCCAGCCGATGAGTAACTCACCTTCGCCAACGGCTGTCATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGtt                                                                                                              <  7:16445/148‑1 (MQ=255)                              aCTCCAGCCGATGAGTAACTCACCTTCGCCAACGGCTGTCATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGt                                                                                                               >  5:371291/1‑147 (MQ=255)                                   aGCCGATGAGTAACTCACCTTCGCCAACGGCTGTCATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATAtt                                                                                                         >  2:235902/1‑148 (MQ=255)                                   aGCCGATGAGTAACTCACCTTCGCCAACGGCTGTCATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATAtt                                                                                                         >  4:28722/1‑148 (MQ=255)                                    gCCGATGAGTAACTCACCTTCGCCAACGGCTGTCATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGTTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATAtt                                                                                                         >  4:207121/1‑147 (MQ=255)                                    gCCGATGAGTAACTCACCTTCGCCAACGGCTGTCATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGTTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATAtt                                                                                                         >  4:202020/1‑147 (MQ=255)                                           aGTAACTCACCTTCGCCAACGGCTGTCATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAg                                                                                                  >  3:365217/1‑147 (MQ=255)                                              aaCTCACCTTCGCCAACGGCTGTCATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTc                                                                                                    >  3:196683/1‑142 (MQ=255)                                              aaCTCACCTTCGCCAACGGCTGTCATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTc                                                                                                    <  4:196683/142‑1 (MQ=255)                                              aaCTCACCTTCGCCAACGGCTGTCATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAAc                                                                                              >  4:40347/1‑148 (MQ=255)                                                  cACCTTCGCCAACGGCTGTCATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGccc                                                                                                                                        <  2:378183/102‑1 (MQ=255)                                                  cACCTTCGCCAACGGCTGTCATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGccc                                                                                                                                        >  1:378183/1‑102 (MQ=255)                                                  cACCTTCGCCAACGGCTGTCATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCg                                                                                          >  7:137486/1‑148 (MQ=255)                                                   aCCTTCGCCAACGGCTGTCATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCgg                                                                                         <  4:365217/148‑1 (MQ=255)                                                   aCCTTCGCCAACGGCTGTCATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCgg                                                                                         <  4:58832/148‑1 (MQ=255)                                                    ccTTCTCCAACGGCGGTCATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATAAACCTCGAATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGt                                                                                        <  7:323831/148‑1 (MQ=255)                                                     cTTCGCCAACGGCTGTCATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTa                                                                                       >  3:34458/1‑148 (MQ=255)                                                     cTTCGCCAACGGCTGTCATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTa                                                                                       >  2:122713/1‑148 (MQ=255)                                                       tCGCCAACGGCTGTCATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTAt                                                                                      >  2:380873/1‑147 (MQ=255)                                                        cGCCAACGGCTGTCATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATc                                                                                     <  6:321338/147‑1 (MQ=255)                                                              cGGCTGTCATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGTTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCAt                                                                              <  3:207121/148‑1 (MQ=255)                                                              cGGCTGTCATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGTTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCAt                                                                              <  3:202020/148‑1 (MQ=255)                                                                       tCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATc                                                                                     <  6:73930/132‑1 (MQ=255)                                                                       tCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATc                                                                                     >  5:73930/1‑132 (MQ=255)                                                                       tCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGcgcg                                                                      <  1:347521/147‑1 (MQ=255)                                                                       tCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCgg                                                                     >  2:123419/1‑148 (MQ=255)                                                                       tCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCAGGTATCAGTTCATCCAGcgcg                                                                      >  2:347521/1‑147 (MQ=255)                                                                        cccTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGgc                                                                    >  4:254656/1‑148 (MQ=255)                                                                            ttCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGgctgct                                                                >  8:260736/1‑148 (MQ=255)                                                                             tCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGgctgct                                                                <  7:260736/147‑1 (MQ=255)                                                                                ggTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGt                                                            <  6:304829/148‑1 (MQ=255)                                                                                 gTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGtt                                                           >  8:243595/1‑148 (MQ=255)                                                                                  tGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTc                                                          <  7:243595/148‑1 (MQ=255)                                                                                    aaTGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGccccc                                                                                                                                      <  6:354544/70‑1 (MQ=255)                                                                                    aaTGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGccccc                                                                                                                                      >  5:354544/1‑70 (MQ=255)                                                                                     aTGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCAt                                                        <  8:336255/147‑1 (MQ=255)                                                                                     aTGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCAt                                                        >  7:336255/1‑147 (MQ=255)                                                                                       gATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGc                                                     >  4:44882/1‑148 (MQ=255)                                                                                           tttCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTtg                                                 <  1:235902/148‑1 (MQ=255)                                                                                           tttCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTtg                                                 <  3:28722/148‑1 (MQ=255)                                                                                             tCCATACGCTGTTGGATATCCGTTGTTACACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTtg                                                 <  4:18002/146‑1 (MQ=255)                                                                                             tCCATACGCTGTTGGATATCCGTTGTTACACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTtg                                                 >  3:18002/1‑146 (MQ=255)                                                                                                aTACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTgatgat                                            <  1:380873/148‑1 (MQ=255)                                                                                                  aCGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCGGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTtga                                                <  4:239102/142‑1 (MQ=255)                                                                                                  aCGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCGGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTtga                                                >  3:239102/1‑142 (MQ=255)                                                                                                        ttGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCACCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCg                                      <  8:137486/146‑1 (MQ=255)                                                                                                         tGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGtt                                    >  7:411311/1‑147 (MQ=255)                                                                                                          ggATATCCGTTGGTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCAGTCAGCGCGGGAATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTc                                  <  4:65752/148‑1 (MQ=255)                                                                                                           gATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGttt                                   >  7:112255/1‑146 (MQ=255)                                                                                                               tCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGa                              >  5:173814/1‑147 (MQ=255)                                                                                                                ccGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTtgatg                                              >  2:349004/1‑130 (MQ=255)                                                                                                                ccGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTtgatg                                              <  1:349004/130‑1 (MQ=255)                                                                                                                  gttgttATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGtt                                    >  7:405070/1‑138 (MQ=255)                                                                                                                  gttgttATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGtt                                    <  8:405070/138‑1 (MQ=255)                                                                                                                  gttgttATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAAt                           >  3:150337/1‑147 (MQ=255)                                                                                                                     gttATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCaa                        <  3:44882/147‑1 (MQ=255)                                                                                                                     gtgataCACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCAGTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCaa                        <  6:86737/144‑1 (MQ=255)                                                                                                                      ttATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAAc                       <  6:371291/147‑1 (MQ=255)                                                                                                                      ttATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAAc                       <  1:122713/147‑1 (MQ=255)                                                                                                                        ataCACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATg                    >  8:217600/1‑148 (MQ=255)                                                                                                                           cacCTCGTATAACTCAAGACATGCTGCCCCCTGAGTATTGTCTGTCAGCGAGGTTATATTCACTAATCAAACACGGTATCAGTTCATCCAGCGAGTCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAAAATGTct                 >  1:319894/1‑147 (MQ=18)                                                                                                                            acCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCgg                <  8:411311/148‑1 (MQ=255)                                                                                                                               tCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTtg                                                 >  7:285793/1‑112 (MQ=255)                                                                                                                               tCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTtg                                                 <  8:285793/112‑1 (MQ=255)                                                                                                                               tCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGttt             >  5:193540/1‑148 (MQ=255)                                                                                                                                  tataACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTcc          <  7:217600/148‑1 (MQ=255)                                                                                                                                   ataACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTc           >  4:96465/1‑146 (MQ=255)                                                                                                                                        tCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTcc    <  3:96465/148‑1 (MQ=255)                                                                                                                                          aaGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCgg  <  3:40347/148‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                            AATGCCGGAGCATGAGCACCGCTTGTACCACTCCAGCCGATGAGTAACTCACCTTCGCCAACGGCTGTCATCCCTTTCAGGTGAATGATATTTCCATACGCTGTTGGATATCCGTTGTTATACACCTCGTATAACTCAAGACCTGCTGCCCCCTGCGTATTGTCTGTCAGCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGG  >  NC_000913/1635663‑1635947

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