Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	727,315	C→T	9.3%	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	9.3%	174.7 / 6.4	64	W6* (TGG→TAG) ‡	PP_0621	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (28/30);  new base T (3/3);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

ACTGCTGGCTGACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGCAGCAGCTCTCGCGG  >  NC_002947/727172‑727461                                                                                                                                                |                                                                                                                                                    aCTGCTGGCTGACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGc                                                                                                                                                <  7:298320/149‑1 (MQ=255)            aCTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCTTGGGCCTCCCttttt                                                                                                                                      <  6:258747/148‑3 (MQ=255)            aCTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTttc                                                                                                                                       <  1:200699/147‑1 (MQ=255)              tttGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGTCCTCCCTttcatt                                                                                                                                    >  4:241644/1‑148 (MQ=255)              tttGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTttcatt                                                                                                                                    >  2:171830/1‑148 (MQ=255)                tGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCAc                                                                                                                                 >  4:250615/1‑149 (MQ=255)                    tttGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCTTGGGCCTCCCttttttattg                                                                                                                                 >  4:213509/1‑138 (MQ=255)                    tttGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCTTGGGCCTCCCttttttattg                                                                                                                                 <  3:213509/145‑8 (MQ=255)                    tttGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCTTGGGCCTCCCttt                                                                                                                                        >  8:63180/1‑138 (MQ=255)                    tttGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCTTGGGCCTCCCttt                                                                                                                                        <  7:63180/138‑1 (MQ=255)                    tttGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTttca                                                                                                                                      <  1:129011/140‑1 (MQ=255)                        cgcCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGaa                                                                                                                          <  3:78443/148‑1 (MQ=255)                        cgcCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAg                                                                                                                         <  2:177281/149‑1 (MQ=255)                          ccACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGt                                                                                                                       <  7:261992/149‑1 (MQ=255)                               gCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTc                                                                                                                  <  7:7287/149‑1 (MQ=255)                                      cgcgGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTttca                                                                                                                                      <  3:93162/122‑1 (MQ=255)                                      cgcgGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTttca                                                                                                                                      >  4:93162/1‑122 (MQ=255)                                       gcggcggcgATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGa                                                                                                          >  7:291712/1‑148 (MQ=255)                                       gcggcggcgATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGAc                                                                                                         <  2:116336/149‑1 (MQ=255)                                        cggcggcgATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGAcc                                                                                                        <  2:247557/149‑1 (MQ=255)                                         ggcggcgATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGagc                                                                                                        >  8:214347/1‑146 (MQ=255)                                           cggcgATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCaaa                                                                                                     <  3:241644/149‑1 (MQ=255)                                              cgttctgggCGATCCAGAAAAGTAGCCCAACCTTGGTGGCCTGCTTGTCTATTGAACAACCAAACAGTCTAGGACGCGCCAAATAAAAACGGGCGCCCCAAGGCCTCCCTTTAATTCACCGCTGAAGGTCAGTCTAA‑CACACAAACggc                                                                                                  <  6:45120/143‑2 (MQ=255)                                               gATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGt                                                                                                  >  7:74241/1‑148 (MQ=255)                                                  cAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTAACCGCTGAAGGTCAGTCGAA‑CAGACAAAAGGTCATg                                                                                              >  4:177465/1‑149 (MQ=255)                                                    gggCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCGAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATg                                                                                              >  5:68842/1‑147 (MQ=255)                                                       cGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCt                                                                                         <  8:74241/149‑1 (MQ=255)                                                             aGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTg                                                                                        <  6:191918/144‑1 (MQ=255)                                                             aGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTGAAc                                                                                     >  6:301164/1‑147 (MQ=255)                                                             aGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTGAACCa                                                                                   <  3:88762/149‑1 (MQ=255)                                                             aGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCATTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTGAAc                                                                                     >  5:235455/1‑147 (MQ=255)                                                             aGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTACCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTg                                                                                        >  5:191918/1‑144 (MQ=255)                                                                aaaGTAGGCGAACCATGGTGTCCAGCTTCTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTGAACCAGGa                                                                                <  8:291712/149‑1 (MQ=255)                                                                      ggCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTc                                                                                                 >  2:22977/1‑126 (MQ=255)                                                                      ggCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTc                                                                                                 <  1:22977/126‑1 (MQ=255)                                                                      ggCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTGAACCAGGAACGGt                                                                           <  1:7622/148‑1 (MQ=255)                                                                      ggCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTGAACCAGGAACGGt                                                                           <  3:250615/148‑1 (MQ=255)                                                                       gCGAACCATGGTGTGCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCt                                                                         >  7:274398/1‑149 (MQ=255)                                                                        cGAACCATGGTGTCCAGCTTGTCGATTGAAGAATCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCAACGCTGAAGTTAAGTAGAA‑CAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCtt                                                                        >  1:293771/1‑149 (MQ=255)                                                                                            gTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGcc                                                     <  3:133621/148‑1 (MQ=255)                                                                                                 ttGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCgg                                               <  1:171830/149‑1 (MQ=255)                                                                                                    aagaagCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCAc                                            >  4:162383/1‑149 (MQ=255)                                                                                                       aagCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGAt                                         >  7:283678/1‑149 (MQ=255)                                                                                                        agCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCTCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGAtt                                        <  6:68842/149‑1 (MQ=255)                                                                                                          ccAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTgg                                      >  7:118918/1‑149 (MQ=255)                                                                                                            aaCAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCa                                             <  7:189548/140‑1 (MQ=255)                                                                                                            aaCAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCa                                             >  8:189548/1‑140 (MQ=255)                                                                                                            aaCAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGt                                    >  2:26952/1‑149 (MQ=255)                                                                                                            aaCAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGt                                    >  8:255591/1‑149 (MQ=255)                                                                                                             aCAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTc                                   >  8:82713/1‑149 (MQ=255)                                                                                                              cAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCg                                  <  3:162383/149‑1 (MQ=255)                                                                                                                  cTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCg                               <  8:274398/148‑1 (MQ=255)                                                                                                                   tAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGtt                             >  8:182961/1‑149 (MQ=255)                                                                                                                    aGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAAGGTGACCAAAGGGCATGTAGCAGAACCAGGAACGGTCTTTCTCGGCTTATTTCGGGCCTTCCGGCACGATTGGGTCGCCGt                              >  1:261844/1‑148 (MQ=255)                                                                                                                     ggACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGCA‑CAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGtt                             <  3:177465/147‑1 (MQ=255)                                                                                                                             ccAATAAAAAATGGAGGCCCATGGCCTACCTTTCATTCACCGATGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTTAACCAGGATCGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGTTCGGCGTTCTCGTCTTTa                   <  2:293771/149‑1 (MQ=255)                                                                                                                             ccAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTa                   <  2:19369/149‑1 (MQ=255)                                                                                                                              cAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAg                  <  2:113944/149‑1 (MQ=255)                                                                                                                              cAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAg                  >  1:113944/1‑149 (MQ=255)                                                                                                                               aaTAAAAAAGGCACGCCCAATGCCTCCCTTTCATTCACCGATGAAGGTCAGTCGAA‑CAGACCAAAGGTCATTTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAgg                 <  7:214347/149‑1 (MQ=255)                                                                                                                                  aaaaaaGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGgcag              >  4:104924/1‑149 (MQ=255)                                                                                                                                    aaaaGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGgcag              >  2:169090/1‑147 (MQ=255)                                                                                                                                      aaGGGAGGCCTTGGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGCagcagc          >  8:76199/1‑149 (MQ=255)                                                                                                                                             gCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGCAGCAGCTCTcgcg   >  1:129956/1‑149 (MQ=255)                                                                                                                                               ccAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGCAGCAGCTCTCGCgg  >  6:201563/1‑148 (MQ=255)                                                                                                                                                |                                                                                                                                                    ACTGCTGGCTGACTTTGTATTTGCGCCACAGCCAGAACGCGGCGGCGATCAGGGCGATCCAGAAAAGTAGGCGAACCATGGTGTCCAGCTTGTCGATTGAAGAAGCCAACAGTCTAGGACGCGGCCAATAAAAAAGGGAGGCCCAAGGCCTCCCTTTCATTCACCGCTGAAGGTCAGTCGAA‑CAGACCAAAGGTCATGTAGCTGAACCAGGAACGGTCTTTCTCGGCTTCTTTCGGGCCTTCCGGCACGATTGGGTCGCCGTTCTCGTCTTTAGGCAGCAGCTCTCGCGG  >  NC_002947/727172‑727461

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