Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	727,315	C→T	10.5%	W6* (TGG→TAG) ‡	PP_0621 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	727,315	0	C	T	10.5%	231.2 / 13.3	87	W6* (TGG→TAG) ‡	PP_0621	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (40/37);  new base T (4/5);  total (45/42)
Fisher's exact test for biased strand distribution p-value = 7.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.61e-01

GTGCGACTGCTGGCTGACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGCAGCAGCTCTCGCGGCA  >  NC_002947/727167‑727463                                                                                                                                                      |                                                                                                                                                     tttcGACTGCTGGCTGAATTTTTATTTGCGCCACAGCCAGAACGCGGCGGAGATCAGGGCGATCCTGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGccc                                                                                                                                                      <  3:263770/146‑1 (MQ=255) gTGCGACTGCTGGCTGACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGccc                                                                                                                                                      <  7:187974/149‑1 (MQ=255)   gCGACTGCTGGCTGACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGccc                                                                                                                                                      >  8:230866/1‑147 (MQ=255)      aCTGCTGGCTGACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGc                                                                                                                                                 <  7:261936/149‑1 (MQ=255)          ctggctgACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCTATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTcc                                                                                                                                             <  7:64018/149‑1 (MQ=255)          ctggctgACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTc                                                                                                                                              <  7:230866/148‑1 (MQ=255)            ggctgACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCt                                                                                                                                           >  7:89318/1‑149 (MQ=255)             gctgACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCTTGGGCCTCCCt                                                                                                                                           <  4:149601/148‑1 (MQ=255)             gctgACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCtt                                                                                                                                          <  8:89318/149‑1 (MQ=255)             gctgACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCt                                                                                                                                           <  3:42125/148‑1 (MQ=255)                gACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTttca                                                                                                                                       >  7:359579/1‑149 (MQ=255)                 aCTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCTGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTTGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCTTGGGCCTCCCttttt                                                                                                                                       <  4:314703/148‑3 (MQ=255)                 aCTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCTTGGGCCTCCCtttt                                                                                                                                        <  7:123768/147‑2 (MQ=255)                 aCTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTttcat                                                                                                                                      >  8:218534/1‑149 (MQ=255)                  cTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAACGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCTTGGGCCTCCCtttttta                                                                                                                                     <  7:18745/149‑5 (MQ=255)                   tttGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCGACGGCCGCGGTtgcattc                                                                                                                                    >  4:72604/1‑149 (MQ=255)                   tttGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTttcatt                                                                                                                                     >  7:221179/1‑148 (MQ=255)                   tttGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTttcatt                                                                                                                                     >  6:125249/1‑148 (MQ=255)                      gTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGCCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCAc                                                                                                                                  <  4:195488/148‑1 (MQ=255)                       tATTTGCGCCACAGCCAGACCGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATTGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCg                                                                                                                                <  2:368869/149‑1 (MQ=255)                         tttGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCt                                                                                                                              >  3:422484/1‑149 (MQ=255)                           tGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTg                                                                                                                             <  5:431610/148‑1 (MQ=255)                             cgcCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGTAAAGTAGGCGTAACATTGTGTCCATCTTTTCGATTGAACAATCCAACAGTCAATGACGCCGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAg                                                                                                                          <  1:399811/149‑1 (MQ=255)                             cgcCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGaa                                                                                                                           <  1:225492/148‑1 (MQ=255)                                cacaGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGt                                                                                                                        >  8:209498/1‑148 (MQ=255)                                cacaGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGt                                                                                                                        >  8:198512/1‑148 (MQ=255)                                 acaGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCTACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCGCCGCTGAAGGt                                                                                                                        >  4:87577/1‑147 (MQ=255)                                  caGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAg                                                                                                                     >  1:423514/1‑149 (MQ=255)                                     ccAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGt                                                                                                                        >  3:328098/1‑143 (MQ=255)                                     ccAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGt                                                                                                                        <  4:328098/143‑1 (MQ=255)                                      cAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCg                                                                                                                  >  1:398033/1‑148 (MQ=255)                                       aGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGaa                                                                                                                <  7:218534/149‑1 (MQ=255)                                       aGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGaa                                                                                                                >  2:169636/1‑149 (MQ=255)                                       aGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGaa                                                                                                                <  4:422484/149‑1 (MQ=255)                                         aaCGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCTGTCGAACa                                                                                                              <  5:196434/149‑1 (MQ=255)                                          aCGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACa                                                                                                              <  7:45223/148‑1 (MQ=255)                                            gcggcggcgATCAGGGCGATCCAGAAAAGAAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGAc                                                                                                           <  2:216840/149‑1 (MQ=255)                                              ggcggcgATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCTTGGGCCTCCCt                                                                                                                                           <  6:61164/115‑1 (MQ=255)                                              ggcggcgATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCTTGGGCCTCCCt                                                                                                                                           >  5:61164/1‑115 (MQ=255)                                                 ggcgATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAg                                                                                                      <  7:209498/149‑1 (MQ=255)                                                    gATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCa                                                                                                                      >  1:181806/1‑130 (MQ=255)                                                    gATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCa                                                                                                                      <  2:181806/130‑1 (MQ=255)                                                          ggCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCa                                                                                                                      >  8:105946/1‑124 (MQ=255)                                                          ggCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCa                                                                                                                      <  7:105946/124‑1 (MQ=255)                                                           gCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGc                                                                                            <  2:423178/149‑1 (MQ=255)                                                               tCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCGTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGaa                                                                                        >  3:228069/1‑149 (MQ=255)                                                                  aGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGACGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAAc                                                                                       >  7:217744/1‑147 (MQ=255)                                                                  aGAAAAGTAGGCGAACCATGGTGTCCAGATTGTCGATTGAAGAAGCCAACAGTCTAGGACGAGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGAACAAAGCTCATGTAGCTGAAc                                                                                       >  6:257322/1‑147 (MQ=255)                                                                     aaaGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAgg                                                                                   >  4:386540/1‑148 (MQ=255)                                                                         tAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCAAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACgg                                                                              >  8:18666/1‑149 (MQ=255)                                                                           ggCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTTATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTc                                                                            <  1:341722/149‑1 (MQ=255)                                                                             cGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCtt                                                                          >  6:368788/1‑149 (MQ=255)                                                                                   aTGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCgg                                                                    >  7:183445/1‑149 (MQ=255)                                                                                    tGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCgg                                                                    >  1:325666/1‑149 (MQ=255)                                                                                    tGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGc                                                                   <  2:423514/149‑1 (MQ=255)                                                                                          ccAGCTTGTCGATTGAAGAAGACAACAGTCAAGGACGCGACCAAT‑ATAAAAGGGAGGCCCAAGTCATCCATTGCCTTCACCGCTGAAGGTCCGACCAACAGACGACAGGTCATGTCGCTTAAACATTATCGGTCTGACTCCGCTTCttt                                                             >  8:150892/1‑149 (MQ=255)                                                                                           cAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGaa                                                                                 <  6:136350/128‑1 (MQ=255)                                                                                           cAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGaa                                                                                 >  5:136350/1‑128 (MQ=255)                                                                                           cAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCttt                                                             <  7:150829/148‑1 (MQ=255)                                                                                            aGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCATTCg                                                           >  5:103955/1‑149 (MQ=255)                                                                                             gCTTGTCGATTGAAGACGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCgg                                                          <  4:108354/149‑1 (MQ=255)                                                                                             gCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCgg                                                          <  4:114369/149‑1 (MQ=255)                                                                                             gCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCgg                                                          <  1:14964/149‑1 (MQ=255)                                                                                                    gATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTcc                                                   >  2:176990/1‑149 (MQ=255)                                                                                                         aagaagCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCAc                                              >  2:109979/1‑149 (MQ=255)                                                                                                         aagaagCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGACCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCAc                                              >  6:248848/1‑149 (MQ=255)                                                                                                           gaagCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGa                                            >  7:137944/1‑149 (MQ=255)                                                                                                           gaagCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCACGGCTTCTTTCGGGCCTTCCGGAACGa                                            >  7:395347/1‑149 (MQ=255)                                                                                                            aagCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGAt                                           >  4:38821/1‑149 (MQ=255)                                                                                                              gCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCgg                                                          >  5:143405/1‑132 (MQ=255)                                                                                                              gCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCgg                                                          <  6:143405/132‑1 (MQ=255)                                                                                                              gCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTg                                         >  8:354037/1‑149 (MQ=255)                                                                                                                 aaCAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGt                                      <  1:176990/149‑1 (MQ=255)                                                                                                                    aGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGc                                   >  6:23637/1‑149 (MQ=255)                                                                                                                      tCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCg                                 >  7:166224/1‑149 (MQ=255)                                                                                                                                  ccAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTa                     <  7:354037/149‑1 (MQ=255)                                                                                                                                  ccAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTa                     <  8:183445/149‑1 (MQ=255)                                                                                                                                   cAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAg                    >  3:301179/1‑149 (MQ=255)                                                                                                                                      t‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGgca                 <  7:198512/149‑1 (MQ=255)                                                                                                                                         tgaaaGGGAGGCCTTGGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGgcag                >  8:292227/3‑148 (MQ=255)                                                                                                                                         aaaaaGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGgcagc               <  3:386540/149‑1 (MQ=255)                                                                                                                                         aaaaaGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGgcag                >  2:307282/1‑148 (MQ=255)                                                                                                                                          aaaaGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGgcagca              <  4:410100/149‑1 (MQ=255)                                                                                                                                           aaaGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGcagcag             <  5:125249/149‑1 (MQ=255)                                                                                                                                             aGGGAGGCCTTGGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGCagcagc            >  7:91088/1‑148 (MQ=255)                                                                                                                                             aGGGAGGCCTTGGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGCAGCAGct           >  7:40687/1‑149 (MQ=255)                                                                                                                                             aGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACC‑AAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTGTAGGCAGCAGctc          >  4:346505/1‑149 (MQ=255)                                                                                                                                              gggAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGt                          <  8:370359/133‑1 (MQ=255)                                                                                                                                              gggAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGt                          >  7:370359/1‑133 (MQ=255)                                                                                                                                                 aGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGCAGCAGCTCTcg       <  2:325666/149‑1 (MQ=255)                                                                                                                                                    cccAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGCAGCAGCTCTCGCgg    <  1:169636/149‑1 (MQ=255)                                                                                                                                                     ccAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGCAGCAGCTCTCGCgg    >  5:225356/1‑148 (MQ=255)                                                                                                                                                      cAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCTTGATTGGGTCGCCGTTCTCGTCTTTAGGCAGCAGCTCTCGCGGCa  >  6:422963/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                     GTGCGACTGCTGGCTGACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAAT‑AAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAACAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGCAGCAGCTCTCGCGGCA  >  NC_002947/727167‑727463

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