Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,082,048	C→A	5.3%	pseudogene (589/726 nt)	efeU →	ferrous iron permease (pseudogene)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,082,048	0	C	A	5.3%	288.4 / 3.7	85	pseudogene (589/726 nt)	efeU	ferrous iron permease (pseudogene)
Reads supporting (aligned to +/- strand):  ref base C (44/36);  new base A (3/2);  total (47/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 = 1.00e+00

ATGGACCAGCCTGTTTATTCTCTTCGTCGCCGCAGGGCTGGCAGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCCGCGTAACAAATACGAC  >  NC_000913/1081900‑1082195                                                                                                                                                     |                                                                                                                                                    aTGGACCAGCCTGTTTATTCTCTTCGTCGCCGCAGGGCTGGCAGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGc                                                                                                                                                     >  1:379281/1‑149 (MQ=255)     aCCAGCCTGTTTATTCTCTTCGTCGCCGCAGGGCTGGCAGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGAt                                                                                                                                                 <  4:284713/149‑1 (MQ=255)         gCCTGTTTATTCTCTTCGTCGCCGCAGGGCTGGCAGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGaa                                                                                                                                             <  8:15858/149‑1 (MQ=255)          ccTGTTTATTCTCTTCGTCGCCGCAGGGCTGGCAGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAg                                                                                                                                            <  1:287606/149‑1 (MQ=255)           cTGTTTATTCTCTTCGTCGCCGCAGGGCTGGCAGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAgg                                                                                                                                           <  4:443211/149‑1 (MQ=255)           cTGTTTATTCTCTTCGTCGCCGCACGGCTGGCAGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTCCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAgg                                                                                                                                           >  4:264621/1‑149 (MQ=255)                 aTTCTCTTCGTCGCCGCAGGGCTGGCAGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGATGATGGAAGGGAtttt                                                                                                                                     >  3:105682/1‑149 (MQ=255)                 aTTCTCTTCGTCGCCGCAGGGCTGGCAGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGATGATGGAAGGGAtttt                                                                                                                                     >  3:111971/1‑149 (MQ=255)                  ttCTCTTCGTCGCCGCAGGGCTGTCAGCTGGTGCAATTCGCGAATTTCATGAAGCCGGATTGTGGAACCACTTTCAGTAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGATGATGGAAGGGAttttt                                                                                                                                    <  4:105682/149‑1 (MQ=255)                  ttCTCTTCGTCGCCGCAGGGCTGGCAGCTGGTGCCTTTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGAttttt                                                                                                                                    <  2:402272/149‑1 (MQ=255)                   tctcTTCGTCGCCGCAGGGCTGGCAGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGATGATGGAAGGGAttttt                                                                                                                                    <  4:111971/148‑1 (MQ=255)                     tctTCGTCGCCGCAGGGCTGGCAGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGACCCACTGTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGc                                                                                                                                 >  5:490175/1‑149 (MQ=255)                       ttCGTCGCCGCAGGGCTGGCAGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCGACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTa                                                                                                                               >  7:51012/1‑149 (MQ=255)                       ttCGTCGCCGCAGGGCTGGCAGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTgg                                                                                                                                  <  8:268728/146‑1 (MQ=255)                          gtcgCCGCAGGGCTGGCAGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCTATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCa                                                                                                                            <  2:246181/149‑1 (MQ=255)                          gtcgCCGCAGGGCTGGCAGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCGACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCa                                                                                                                            <  8:51012/149‑1 (MQ=255)                          gtcgCCGCAGGGCTGGCAGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTAt                                                                                                                              >  1:246181/1‑147 (MQ=255)                              ccgcAGGGCTGGCAGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGaa                                                                                                                        <  7:330511/149‑1 (MQ=255)                               cgcAGGGCTGGCAGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAg                                                                                                                       >  2:318287/1‑149 (MQ=255)                                       tGGCAGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGc                                                                                                               >  6:123462/1‑149 (MQ=255)                                           aGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGATGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTcc                                                                                                           >  6:115127/1‑148 (MQ=255)                                             cTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGc                                                                                                         <  5:424965/149‑1 (MQ=255)                                              tGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGc                                                                                                         >  8:216271/1‑148 (MQ=255)                                                 tGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGt                                                                                                             <  1:478065/141‑1 (MQ=255)                                                 tGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGt                                                                                                             >  2:478065/1‑141 (MQ=255)                                                     aTTAGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGc                                                                                                 >  6:488536/1‑149 (MQ=255)                                                      ttCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGcc                                                                                                >  4:192333/1‑149 (MQ=255)                                                         gcgcATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCACCTCACTCGCTGTTTGGCACGCTGATGGAAGGGAtttttca                                                                                                                                  >  6:227212/1‑110 (MQ=255)                                                                   tGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGtttatt                                                                                    >  5:390682/1‑148 (MQ=255)                                                                      aGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTAtc                                                                                >  3:129880/1‑149 (MQ=255)                                                                       gCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTAtct                                                                               >  2:30508/1‑149 (MQ=255)                                                                        ccGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTAtct                                                                               >  8:165277/1‑148 (MQ=255)                                                                         cGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGTGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCa                                                                             <  1:318287/149‑1 (MQ=255)                                                                          ggATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCAt                                                                            >  2:236648/1‑149 (MQ=255)                                                                          ggATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCAt                                                                            >  4:208941/1‑149 (MQ=255)                                                                                  gAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGgcgc                                                                    >  8:288151/1‑149 (MQ=255)                                                                                   aaCCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCt                                                                   >  7:294316/1‑149 (MQ=255)                                                                                    aCCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGGCGCCGTCTGGTTTATTTATCTCATCCCGGCGCtg                                                                  <  3:264621/149‑1 (MQ=255)                                                                                            cAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCAtt                                                          >  4:424649/1‑149 (MQ=255)                                                                                              ggAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTg                                                        <  3:192333/149‑1 (MQ=255)                                                                                               gAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGc                                                       <  6:490175/149‑1 (MQ=255)                                                                                               gAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCAACCCGGCGCTGGTGGCATTTGc                                                       >  7:139860/1‑149 (MQ=255)                                                                                                 aaTCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGctc                                                     <  1:30508/149‑1 (MQ=255)                                                                                                  aTCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGctct                                                    <  8:294316/149‑1 (MQ=255)                                                                                                     gCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGcc                                                 >  8:253695/1‑149 (MQ=255)                                                                                                        ttCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGccacc                                              >  4:387863/1‑149 (MQ=255)                                                                                                           gATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCAcg                                           >  1:60262/1‑149 (MQ=255)                                                                                                            atatGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCAcgc                                          <  5:123462/149‑1 (MQ=255)                                                                                                              atGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCAcgcg                                         >  8:406317/1‑148 (MQ=255)                                                                                                               tGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCa                                       <  7:216271/149‑1 (MQ=255)                                                                                                               tGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCa                                       <  3:143065/149‑1 (MQ=255)                                                                                                                gAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAg                                      >  3:426156/1‑149 (MQ=255)                                                                                                                 aGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAg                                      <  7:288151/148‑1 (MQ=255)                                                                                                                  gTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAggg                                    <  5:211433/149‑1 (MQ=255)                                                                                                                      ggTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATATCATCCCGgc                                                                      >  1:157883/1‑111 (MQ=255)                                                                                                                      ggTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGATATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCAGTCTGGTTTATTTATCTCATCCCGgc                                                                      <  2:157883/111‑1 (MQ=255)                                                                                                                       gTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCt                                                                                            >  2:302656/1‑88 (MQ=255)                                                                                                                       gTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCt                                                                                            <  1:302656/88‑1 (MQ=255)                                                                                                                           tctcAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCg                           >  5:412948/1‑149 (MQ=255)                                                                                                                             tcAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGtc                         >  3:174320/1‑149 (MQ=255)                                                                                                                             tcAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGtc                         >  2:435582/1‑149 (MQ=255)                                                                                                                             tcAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGtc                         >  3:221132/1‑149 (MQ=255)                                                                                                                             tcAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGtc                         >  8:319499/1‑149 (MQ=255)                                                                                                                               aaCTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCAttt                                                         <  4:428951/115‑1 (MQ=255)                                                                                                                               aaCTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCAttt                                                         >  3:428951/1‑115 (MQ=255)                                                                                                                               aaCTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGa                                <  3:323880/140‑1 (MQ=255)                                                                                                                               aaCTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGa                                >  4:323880/1‑140 (MQ=255)                                                                                                                               aaCTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGtctc                       >  5:11599/1‑149 (MQ=255)                                                                                                                               aaCTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGtc                         <  3:208941/147‑1 (MQ=255)                                                                                                                                actcactcGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCg                      <  6:11599/149‑1 (MQ=255)                                                                                                                                 ctcactcGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAgggg                                   >  5:92709/1‑135 (MQ=255)                                                                                                                                 ctcactcGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAgggg                                   <  6:92709/135‑1 (MQ=255)                                                                                                                                  tcactcGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCt                    >  1:458470/1‑149 (MQ=255)                                                                                                                                    actcGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTGGATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTcc                  >  2:244279/1‑149 (MQ=255)                                                                                                                                     ctcGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCcg                 >  2:103679/1‑149 (MQ=255)                                                                                                                                     ctcGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCcg                 >  6:240427/1‑149 (MQ=255)                                                                                                                                      tcGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGCCAGCGTCTCGCTCcgc                >  6:119499/1‑149 (MQ=255)                                                                                                                                      tcGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCcgc                >  5:225587/1‑149 (MQ=255)                                                                                                                                      tcGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCcgc                >  2:329010/1‑149 (MQ=255)                                                                                                                                        gCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCcgc                <  8:257384/147‑1 (MQ=255)                                                                                                                                        gCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCcgc                <  6:390682/147‑1 (MQ=255)                                                                                                                                        gCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCACCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCCGCGt              <  8:139860/149‑1 (MQ=255)                                                                                                                                            tttGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGGCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCCGCGTAAc           <  2:379281/148‑1 (MQ=255)                                                                                                                                                gCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCCGCGTAACAAATa      <  7:253695/149‑1 (MQ=255)                                                                                                                                                 cACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCCGCGTAACAAATac     <  4:129880/149‑1 (MQ=255)                                                                                                                                                    gCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCCGCGTAACAAATacgac  <  3:387863/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    ATGGACCAGCCTGTTTATTCTCTTCGTCGCCGCAGGGCTGGCAGCTGGTGCCATTCGCGCATTTCATGAAGCCGGATTGTGGAACCACTTTCAGGAAATCGCCTTCGATATGAGTGCGGTGCTCTCAACTCACTCGCTGTTTGGCACGCTGATGGAAGGGATTTTTGGCTATCAGGAAGCGCCGAGCGTCAGCGAAGTCGCCGTCTGGTTTATTTATCTCATCCCGGCGCTGGTGGCATTTGCTCTGCCACCACGCGCAGGGGCGACAGCGTCTCGCTCCGCGTAACAAATACGAC  >  NC_000913/1081900‑1082195

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: ‑