Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,069,029	C→T	5.6%	A271T (GCC→ACC)	rutG ←	pyrimidine permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,069,029	0	C	T	5.6%	236.0 / 3.5	71	A271T (GCC→ACC)	rutG	pyrimidine permease
Reads supporting (aligned to +/- strand):  ref base C (32/35);  new base T (2/2);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-01

TGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGC  >  NC_000913/1068888‑1069176                                                                                                                                              |                                                                                                                                                    tGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGcc                                                                                                                                              >  4:393049/1‑149 (MQ=255) tGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGcc                                                                                                                                              >  1:435442/1‑149 (MQ=255)   gTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCtt                                                                                                                                            <  3:85774/149‑1 (MQ=255)    tCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCtt                                                                                                                                            <  2:434244/148‑1 (MQ=255)    tCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTg                                                                                                                                           >  4:57857/1‑149 (MQ=255)    tCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTg                                                                                                                                           >  4:286698/1‑149 (MQ=255)    tCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTg                                                                                                                                           <  1:361535/149‑1 (MQ=255)    tCACGGCCATCACACCGATATGTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCtt                                                                                                                                            >  1:434244/1‑148 (MQ=255)        ggCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGAt                                                                                                                                       >  4:470078/1‑149 (MQ=255)          ccATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGa                                                                                                                                     >  2:48017/1‑149 (MQ=255)            aTCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGAcc                                                                                                                                   <  1:295219/149‑1 (MQ=255)                 acCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGt                                                                                                                              >  4:212118/1‑149 (MQ=255)                  ccGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGtt                                                                                                                             >  3:221873/1‑149 (MQ=255)                   cGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCtt                                                                                                                                            >  3:427905/1‑133 (MQ=255)                   cGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCtt                                                                                                                                            <  4:427905/133‑1 (MQ=255)                    gATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGtttt                                                                                                                           >  4:358777/1‑149 (MQ=255)                     atatTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTc                                                                                                                          >  2:368519/1‑149 (MQ=255)                     atatTTCCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCACCGTTGCAAGTCCATCGCCTACGAATGCCCGCCCCATACAAGGGTCCCTGTTGCGTCCGGTCATCGCGGCGGCCGCCTTGAGATGGACTAAGtccgc                                                                                                                          <  2:250480/149‑5 (MQ=255)                                    gtggtCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAAt                                                                                                              >  4:406499/1‑146 (MQ=255)                                    gtggtCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAAt                                                                                                              <  3:406499/146‑1 (MQ=255)                                    gtggtCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCAc                                                                                                           <  3:212118/149‑1 (MQ=255)                                     tggtCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCAcc                                                                                                          <  3:358777/149‑1 (MQ=255)                                                   ccgccgACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAAt                                                                                            <  4:92007/149‑1 (MQ=255)                                                       cgACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGc                                                                                        <  3:204069/149‑1 (MQ=255)                                                           agagCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTgag                                                                                                                                         >  3:84392/1‑96 (MQ=255)                                                           agagCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTgag                                                                                                                                         <  4:84392/96‑1 (MQ=255)                                                           agagCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGcatca                                                                                    >  2:219551/1‑149 (MQ=255)                                                                 ggAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCt                                                                              >  4:291281/1‑149 (MQ=255)                                                                       cATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCAt                                                                        <  3:286698/149‑1 (MQ=255)                                                                       cATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCAt                                                                        <  1:363552/149‑1 (MQ=255)                                                                        aTCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCAtt                                                                       >  4:447176/1‑149 (MQ=255)                                                                              gCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGgc                                                                                                 <  3:423694/117‑1 (MQ=255)                                                                              gCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGgc                                                                                                 >  4:423694/1‑117 (MQ=255)                                                                                cAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGTGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAg                                                               >  4:252602/1‑149 (MQ=255)                                                                                cAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAg                                                               <  4:221873/149‑1 (MQ=255)                                                                                  gACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGa                                                             <  3:394204/149‑1 (MQ=255)                                                                                  gACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGa                                                             <  3:264875/149‑1 (MQ=255)                                                                                    ccATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGTGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGt                                                           <  3:252602/149‑1 (MQ=255)                                                                                     cATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTgaga                                                                                                                                        >  1:137057/1‑71 (MQ=255)                                                                                     cATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTgaga                                                                                                                                        <  2:137057/71‑1 (MQ=255)                                                                                     cATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTg                                                          <  3:291281/149‑1 (MQ=255)                                                                                     cATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTg                                                          >  2:26736/1‑149 (MQ=255)                                                                                               gAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAAtgtg                                                >  2:208579/1‑149 (MQ=255)                                                                                               gAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAAtgtg                                                <  1:219551/149‑1 (MQ=255)                                                                                                aaTGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTgg                                               >  1:180200/1‑149 (MQ=255)                                                                                                aaTGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTgg                                               <  4:280109/149‑1 (MQ=255)                                                                                                     ccgcccCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTgg                                               <  3:162057/144‑1 (MQ=255)                                                                                                     ccgcccCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTgg                                               >  4:162057/1‑144 (MQ=255)                                                                                                         cccATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGaa                                      <  3:393049/149‑1 (MQ=255)                                                                                                           cATATACGGGTCCATATTGCGTCCGGTCATCCCGGTGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGcatc                                                                                     <  3:338028/100‑1 (MQ=255)                                                                                                           cATATACGGGTCCATATTGCGTCCGGTCATCCCGGTGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGcatc                                                                                     >  4:338028/1‑100 (MQ=255)                                                                                                             tataCGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAg                                  <  1:154400/149‑1 (MQ=255)                                                                                                             tataCGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAg                                  >  2:63600/1‑149 (MQ=255)                                                                                                                  gggTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGgcagc                             <  3:470078/149‑1 (MQ=255)                                                                                                                       cATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGc                        <  1:26736/149‑1 (MQ=255)                                                                                                                             gCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACga                  <  3:447176/149‑1 (MQ=255)                                                                                                                             gCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACga                  <  1:48017/149‑1 (MQ=255)                                                                                                                             gCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACAAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACga                  <  3:274416/149‑1 (MQ=255)                                                                                                                               gTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACgaga                >  4:131894/1‑149 (MQ=255)                                                                                                                               gTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACgaga                <  2:180200/149‑1 (MQ=255)                                                                                                                               gTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACgaga                >  1:475826/1‑149 (MQ=255)                                                                                                                                tCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCAcc                                                                                                          >  2:237525/1‑58 (MQ=255)                                                                                                                                tCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCAcc                                                                                                          <  1:237525/58‑1 (MQ=255)                                                                                                                                  cGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTa             <  1:208579/149‑1 (MQ=255)                                                                                                                                    ttCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTaaa           >  4:278168/2‑149 (MQ=255)                                                                                                                                    ttCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTaaa           >  2:75999/2‑149 (MQ=255)                                                                                                                                    gTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGaa                                      >  1:358541/1‑122 (MQ=255)                                                                                                                                    gTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGaa                                      <  2:358541/122‑1 (MQ=255)                                                                                                                                        tCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCg       >  3:36727/1‑149 (MQ=255)                                                                                                                                          ccGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGAGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGAc     >  1:82258/1‑149 (MQ=255)                                                                                                                                             gCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGc  <  1:175635/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                    TGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGC  >  NC_000913/1068888‑1069176

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