Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,444,424	G→T	13.2%	R156S (CGC→AGC)	sulA ←	cell division inhibitor SulA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,444,424	0	G	T	13.2%	179.0 / 17.1	68	R156S (CGC→AGC)	sulA	cell division inhibitor SulA
Reads supporting (aligned to +/- strand):  ref base G (30/29);  new base T (3/6);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.68e-01

CCATGCTCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTG  >  NC_002947/2444277‑2444568                                                                                                                                                     |                                                                                                                                                   ccATGCTCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCgg                                                                                                                                                   <  1:247714/149‑1 (MQ=255)       tCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGtt                                                                                                                                              >  5:42850/1‑148 (MQ=255)         aGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCCGCCTCGGCTGATGTTCa                                                                                                                                            <  5:122476/148‑1 (MQ=255)         aGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCCGCCTCGGCTGATGTTCAg                                                                                                                                           <  8:157731/149‑1 (MQ=255)         aGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCCGCCTCGGCTGATGTTCAg                                                                                                                                           <  5:202484/149‑1 (MQ=255)         aGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTc                                                                                                                                             <  7:29790/147‑1 (MQ=255)           tgtgGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGc                                                                                                                                         >  3:301603/1‑149 (MQ=255)             tgGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTTAGTGTGCAGCAGCGTATCAGGCTATTTGCAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCt                                                                                                                                        <  1:71884/149‑1 (MQ=255)             tgGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCt                                                                                                                                        <  4:163969/148‑1 (MQ=255)              gggCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCCGCCTCGGCTGATGTTc                                                                                                                                             >  5:213112/1‑142 (MQ=255)              gggCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCCGCCTCGGCTGATGTTc                                                                                                                                             <  6:213112/142‑1 (MQ=255)               ggCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCttt                                                                                                                                      <  8:16692/148‑1 (MQ=255)                 cGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTtgtg                                                                                                                                   <  7:63886/149‑1 (MQ=255)                         tGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTcc                                                                                                                            <  4:137271/148‑1 (MQ=255)                            ccGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGc                                                                                                                         >  1:45929/1‑148 (MQ=255)                             cGGTCTCCTGATGGATGGCGATCTCGAGCGTGATGAAGGGGCGCGTTTGGTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTT‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCa                                                                                                                        <  2:45929/148‑1 (MQ=255)                                tCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTgg                                                                                                                    <  6:122567/149‑1 (MQ=255)                                    cTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGc                                                                                                                <  3:239710/149‑1 (MQ=255)                                     tGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCa                                                                                                               >  8:158151/1‑149 (MQ=255)                                               cGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGAACGTAGCAgtt                                                                                                     >  4:251437/1‑149 (MQ=255)                                                 aGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAgttgt                                                                                                   <  5:233040/149‑1 (MQ=255)                                                  gCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCa                                                                                                                                            <  5:109995/107‑1 (MQ=255)                                                  gCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCa                                                                                                                                            >  6:109995/1‑107 (MQ=255)                                                   cTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAgttgtt                                                                                                  <  6:5147/148‑1 (MQ=255)                                                            gatgaAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGc                                                                                        >  6:144134/1‑149 (MQ=255)                                                                    gCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGa                                                                              >  5:176468/1‑149 (MQ=255)                                                                     ccGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGAc                                                                             >  7:293221/1‑149 (MQ=255)                                                                            tgctgcACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAACGTCA‑‑ATGCCGCTGcac                                                                      <  7:85570/149‑3 (MQ=255)                                                                               tgcACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCccagc                                                                   >  3:208378/1‑149 (MQ=255)                                                                               tgcACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCccagc                                                                   >  2:103448/1‑149 (MQ=255)                                                                                 cACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCccagcc                                                                  <  6:42850/148‑1 (MQ=255)                                                                                    tGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCt                                                              >  7:81224/1‑149 (MQ=255)                                                                                      aGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTca                                                            >  4:249/1‑149 (MQ=255)                                                                                             agcagcGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCATGgt                                                     <  4:301603/149‑1 (MQ=255)                                                                                               cagcGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACTTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGgtgt                                                   >  5:269803/1‑149 (MQ=255)                                                                                                 tcGTATCAGGCTATTTG‑AAGTTACAGGCAAAGCATGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCAAGTAGCAGTTGTTGGAGTGCAGCGGAA‑‑TTGACGTTGCCCAGCCAGCTCAACACGGtgtgt                                                 <  6:269803/148‑2 (MQ=255)                                                                                                    tATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTg                                              <  8:293221/149‑1 (MQ=255)                                                                                                       cAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAgg                                           >  2:256680/1‑149 (MQ=255)                                                                                                       cAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAgg                                           >  6:308069/1‑149 (MQ=255)                                                                                                        aGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑tt                                                                                >  8:258196/1‑111 (MQ=255)                                                                                                        aGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑tt                                                                                <  7:258196/111‑1 (MQ=255)                                                                                                         ggCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGcc                                         >  4:147943/1‑149 (MQ=255)                                                                                                         ggCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGcc                                         >  2:234462/1‑149 (MQ=255)                                                                                                          gCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCcca                                                                     <  4:123760/120‑1 (MQ=255)                                                                                                          gCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCcca                                                                     >  3:123760/1‑120 (MQ=255)                                                                                                          gCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGcc                                         >  2:72602/1‑148 (MQ=255)                                                                                                                     gTGACAGGCAAAGGCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCaaa                               <  7:158151/148‑1 (MQ=255)                                                                                                                     gTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAg                              >  5:5664/1‑149 (MQ=255)                                                                                                                      tGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGc                             >  4:88147/1‑149 (MQ=255)                                                                                                                           ggCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCAGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGACCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTca                         >  5:255675/1‑148 (MQ=255)                                                                                                                            gCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCt                                               >  3:124811/1‑125 (MQ=255)                                                                                                                            gCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCt                                               <  4:124811/125‑1 (MQ=255)                                                                                                                             cAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAg                      <  5:308069/149‑1 (MQ=255)                                                                                                                              aaaGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAgg                     >  7:11466/1‑149 (MQ=255)                                                                                                                               aaGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATCTTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAg                      <  1:256680/149‑1 (MQ=255)                                                                                                                               aaGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCaaa                               <  3:178191/138‑1 (MQ=255)                                                                                                                               aaGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCaaa                               >  4:178191/1‑138 (MQ=255)                                                                                                                               aaGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGc                    <  8:81224/149‑1 (MQ=255)                                                                                                                               aaGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGc                    >  3:286845/1‑149 (MQ=255)                                                                                                                                   cTGAACATCCGCCTCGGCTGATGTTCAGGCTTTGCGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAg                              >  6:146660/1‑135 (MQ=39)                                                                                                                                   cTGAACATCCGCCTCGGCTGATGTTCAGGCTTTGCGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAg                              <  5:146660/135‑1 (MQ=39)                                                                                                                                    tGAACATCCGCCTCGGCTGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCaaa                               <  7:283238/133‑1 (MQ=255)                                                                                                                                    tGAACATCCGCCTCGGCTGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCaaa                               >  8:283238/1‑133 (MQ=255)                                                                                                                                       aCATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCa            <  6:176468/149‑1 (MQ=255)                                                                                                                                         aTCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGt                                    <  2:121824/123‑1 (MQ=255)                                                                                                                                         aTCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGt                                    >  1:121824/1‑123 (MQ=255)                                                                                                                                         aTCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAgc          <  3:251437/149‑1 (MQ=255)                                                                                                                                            aGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAgcgcg       >  8:30142/1‑149 (MQ=255)                                                                                                                                                 aGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGTTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTg  <  1:234462/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   CCATGCTCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTG  >  NC_002947/2444277‑2444568

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