Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,444,415	A→C	14.5%	*159G (TGA→GGA)	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,415	0	A	C	14.5%	175.5 / 20.0	69	*159G (TGA→GGA)	sulA	cell division inhibitor SulA
Reads supporting (aligned to +/- strand):  ref base A (30/29);  new base C (4/6);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 7.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

GTGGCCAAGGCCATGCTCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCG  >  NC_002947/2444267‑2444563                                                                                                                                                      |                                                                                                                                                       gTGGCCAAGGCCATGCTCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCa                                                                                                                                                        <  1:118949/149‑1 (MQ=255)     ccAAGGCCATGCTCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGcc                                                                                                                                                     >  3:154502/1‑148 (MQ=255)       aaGGCCATGCTCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCCGcctcg                                                                                                                                                  >  6:122476/1‑146 (MQ=255)           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                                                                                                                                        <  6:213112/142‑1 (MQ=255)                        gggCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCCGCCTCGGCTGATGTTc                                                                                                                                        >  5:213112/1‑142 (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                                                                                                                                       >  6:109995/1‑107 (MQ=255)                                                            gCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCa                                                                                                                                       <  5:109995/107‑1 (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                                                              >  2:103448/1‑149 (MQ=255)                                                                                         tgcACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCccagc                                                              >  3:208378/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                                      >  6:308069/1‑149 (MQ=255)                                                                                                                 cAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAgg                                      >  2:256680/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                                    >  2:234462/1‑149 (MQ=255)                                                                                                                   ggCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGcc                                    >  4:147943/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                                          <  4:124811/125‑1 (MQ=255)                                                                                                                                      gCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCt                                          >  3:124811/1‑125 (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                          >  4:178191/1‑138 (MQ=255)                                                                                                                                         aaGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCaaa                          <  3:178191/138‑1 (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)                                                                                                                                                      |                                                                                                                                                       GTGGCCAAGGCCATGCTCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTG‑AAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCG  >  NC_002947/2444267‑2444563

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