Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	355,954	A→G	5.0%	M345V (ATG→GTG)	codB →	cytosine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	355,954	0	A	G	5.0%	246.1 / 2.5	79	M345V (ATG→GTG)	codB	cytosine transporter
Reads supporting (aligned to +/- strand):  ref base A (41/34);  new base G (2/2);  total (43/36)
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.49e-01

CCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGA  >  NC_000913/355805‑356101                                                                                                                                                      |                                                                                                                                                    ccGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCtga                                                                                                                                                     >  1:1231916/1‑150 (MQ=255) ccGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCtga                                                                                                                                                     <  2:460087/150‑1 (MQ=255)  cGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCtgat                                                                                                                                                    <  1:1079259/150‑1 (MQ=255)  cGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCtgat                                                                                                                                                    <  1:499219/150‑1 (MQ=255)   gggATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCtga                                                                                                                                                     >  2:1350557/1‑148 (MQ=255)   gggATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCtga                                                                                                                                                     <  1:1350557/148‑1 (MQ=255)        gTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCg                                                                                                                                              >  1:267225/1‑150 (MQ=255)          cGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGGTGAACCGTc                                                                                                                                            >  1:788949/1‑150 (MQ=255)          cGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGGTGAACCGTc                                                                                                                                            >  1:783704/1‑150 (MQ=255)              cAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTcgccg                                                                                                                                        <  1:1097372/150‑1 (MQ=255)                   cccTTTCGGTGATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATg                                                                                                                                   <  2:1518219/150‑1 (MQ=255)                   cccTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCAGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATg                                                                                                                                   >  2:19033/1‑150 (MQ=255)                     cTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAg                                                                                                                                 <  1:1229284/150‑1 (MQ=255)                      tttCGGTAATCAACGGTATTTTCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGc                                                                                                                                <  1:630874/150‑1 (MQ=255)                      tttCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGc                                                                                                                                <  2:613137/150‑1 (MQ=255)                          ggTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGCCTATCTGATGAACCGTCGCCGCTATGAGCACtt                                                                                                                            >  2:832108/1‑150 (MQ=255)                            tAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTg                                                                                                                          >  1:675004/1‑150 (MQ=255)                              aTCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCg                                                                                                                        <  2:767944/150‑1 (MQ=255)                                cAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGAc                                                                                                                      <  2:1114986/150‑1 (MQ=255)                                 aaCGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGAcc                                                                                                                     >  1:541368/1‑150 (MQ=255)                                         tatCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTat                                                                                                             >  2:485716/1‑150 (MQ=255)                                          atCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTatg                                                                                                            <  1:1010983/150‑1 (MQ=255)                                             ggTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTatgatg                                                                                                         <  2:471667/150‑1 (MQ=255)                                             ggTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTatgatg                                                                                                         <  2:1147401/150‑1 (MQ=255)                                                     cTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCaa                                                                                                 >  1:608015/1‑150 (MQ=255)                                                      tGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAAt                                                                                                >  1:1274816/1‑150 (MQ=255)                                                       gcgcATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAAtt                                                                                               >  1:1322739/1‑150 (MQ=255)                                                       gcgcATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCACTGGGTGGCGTTATCATCGCCGACTATCTGATTAAACGTCGCAGATATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAAtt                                                                                               >  2:491079/1‑150 (MQ=255)                                                        cgcATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTg                                                                                              <  2:675004/150‑1 (MQ=255)                                                          cATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTggg                                                                                            >  2:1083530/1‑150 (MQ=255)                                                               tGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCg                                                                                       >  2:341571/1‑150 (MQ=255)                                                                   tGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTc                                                                                   <  1:297086/150‑1 (MQ=255)                                                                   tGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTc                                                                                   >  1:1176520/1‑150 (MQ=255)                                                                        aaCAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCg                                                                              >  1:216072/1‑150 (MQ=255)                                                                        aaCAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTCATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTCAGACCACGCGTATGATGAGAGGCAAGTGCGTGGCGAATCTgcag                                                                              >  2:668694/1‑147 (MQ=255)                                                                            aTTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCg                                                                          >  1:1393010/1‑150 (MQ=255)                                                                                   cGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTggg                                                                   >  2:1248563/1‑150 (MQ=255)                                                                                   cGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTggg                                                                   <  2:1322739/150‑1 (MQ=255)                                                                                    ggCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTgggg                                                                  <  1:1240582/150‑1 (MQ=255)                                                                                    ggCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTgggg                                                                  >  2:952013/1‑150 (MQ=255)                                                                                    ggCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTgggg                                                                  >  2:1402772/1‑150 (MQ=255)                                                                                    ggCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTgggg                                                                  <  2:216072/150‑1 (MQ=255)                                                                                    ggCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTgggg                                                                  <  1:341571/150‑1 (MQ=255)                                                                                     gCTGGTTGACCTTCCTTTCGGCAGCTTTTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCGATGAGCACTATGCGCCCACGCGTAGGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCGTGGGGa                                                                 >  2:1487944/1‑150 (MQ=255)                                                                                      cTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGAt                                                                >  2:838457/1‑150 (MQ=255)                                                                                            gACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGGTGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATtgctgc                                                          <  2:788949/150‑1 (MQ=255)                                                                                            gACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGGTGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATtgctgc                                                          <  2:783704/150‑1 (MQ=255)                                                                                            gACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCgg                                                                             <  2:546209/131‑1 (MQ=255)                                                                                            gACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCgg                                                                             >  1:546209/1‑131 (MQ=255)                                                                                            gACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATtgctgc                                                          <  1:1248563/150‑1 (MQ=255)                                                                                            gACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATtgctgc                                                          <  2:267225/150‑1 (MQ=255)                                                                                            gACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATtgctgc                                                          <  2:1375262/150‑1 (MQ=255)                                                                                              ccttccttTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAg                                                        >  1:1044608/1‑150 (MQ=255)                                                                                               cttccttTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAgg                                                       <  1:176240/150‑1 (MQ=255)                                                                                               cttccttTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAgg                                                       >  1:356839/1‑150 (MQ=255)                                                                                                 tccttTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGcc                                                     <  1:824284/150‑1 (MQ=255)                                                                                                 tccttTCGGCAGCTATTCCTCCAGTGGGTGGAGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACTCGTATGATGAGTGACAATTGGGTGGCGATTCTGGAGGTCGCCTTGGGGATTGCTGCAGGcc                                                     >  1:1438103/1‑150 (MQ=255)                                                                                                      tCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTgg                                                <  1:759299/150‑1 (MQ=255)                                                                                                        ggCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGtt                                              >  2:615982/1‑150 (MQ=255)                                                                                                          cAGGTATTCCTCCAGTGGGGGGCGTGATCATCGCCGACTATCTGATGATCCGTCGCCGCTATGAGGACTTTGCGACCACGCGTATGAGGAGTGTCAATTGGGTGGCTATTCTGGCGGTCGCCTTGGTGATTGCTGCAGGCCACTGGTTAc                                            <  1:1300106/150‑1 (MQ=255)                                                                                                          cAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTAc                                            >  1:1002118/1‑150 (MQ=255)                                                                                                                   ctccAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTg                                   >  1:1391122/1‑150 (MQ=255)                                                                                                                     ccAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGtt                                 >  1:1432892/1‑150 (MQ=255)                                                                                                                      cAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGcc                                                     >  2:536186/1‑129 (MQ=255)                                                                                                                      cAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGcc                                                     >  2:541452/1‑129 (MQ=255)                                                                                                                      cAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGcc                                                     <  1:541452/129‑1 (MQ=255)                                                                                                                      cAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGcc                                                     <  1:536186/129‑1 (MQ=255)                                                                                                                       aGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTcc                               >  2:1267009/1‑150 (MQ=255)                                                                                                                       aGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTcc                               <  1:485716/150‑1 (MQ=255)                                                                                                                         tgggtggCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCgg                             <  2:645486/150‑1 (MQ=255)                                                                                                                          gggtggCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGt                            <  1:156781/150‑1 (MQ=255)                                                                                                                           ggtggCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTc                           >  1:900625/1‑150 (MQ=255)                                                                                                                               ccctGATCATCAACGAATATCAGAAGAACCGAAGCCAATATGAGCACTTTGCGACCAAGAGTATGATGAGTGTCAATTGGGTGGCGATTCTCGCGGTAGCCAAGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTGCCGGTCAAcc                       >  2:1448903/4‑149 (MQ=255)                                                                                                                                 gTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAAcgcg                     >  2:1347108/1‑150 (MQ=255)                                                                                                                                  tGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGAGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCgg                    <  2:1057970/150‑1 (MQ=255)                                                                                                                                   gATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGt                   <  2:541368/150‑1 (MQ=255)                                                                                                                                   gATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACAACGCGTATGATCAGAGTCAATTGGGTGCCGATTCTGGCGGTCGCCTTCGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGt                   >  2:640976/1‑150 (MQ=255)                                                                                                                                    atcatcGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTa                  >  1:871930/1‑150 (MQ=255)                                                                                                                                      catcGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTAtt                >  2:202449/1‑150 (MQ=255)                                                                                                                                       atcGCCGACTATCTGATGACCCGTCGCCGCTATGAGCACTTTGAGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCAGCAGGCCCCTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTa               >  2:468825/1‑150 (MQ=255)                                                                                                                                         cGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAgg             <  2:1002118/150‑1 (MQ=255)                                                                                                                                          gCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGt            >  2:638824/1‑150 (MQ=255)                                                                                                                                                   cTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGaa                                      <  2:1009181/115‑1 (MQ=255)                                                                                                                                                   cTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGaa                                      >  1:1009181/1‑115 (MQ=255)                                                                                                                                                   cTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTg   <  2:1044608/150‑1 (MQ=255)                                                                                                                                                    tgatgaACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGa  >  1:978003/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                    CCGGGATGTCGAGCAAAACCCTTTCGGTAATCAACGGTATTATCGGTACGGTCTGCGCATTATGGCTGTATAACAATTTTGTCGGCTGGTTGACCTTCCTTTCGGCAGCTATTCCTCCAGTGGGTGGCGTGATCATCGCCGACTATCTGATGAACCGTCGCCGCTATGAGCACTTTGCGACCACGCGTATGATGAGTGTCAATTGGGTGGCGATTCTGGCGGTCGCCTTGGGGATTGCTGCAGGCCACTGGTTACCGGGAATTGTTCCGGTCAACGCGGTATTAGGTGGCGCGCTGA  >  NC_000913/355805‑356101

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