Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,391,330	C→T	8.2%	L83F (CTC→TTC)	PP_5707 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,391,330	0	C	T	8.2%	212.2 / 6.6	73	L83F (CTC→TTC)	PP_5707	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/34);  new base T (3/3);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

CCAAAACGCCCGCGAACAGGCCTCGGAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGCGCCGCAGCATCGCCATCAA  >  NC_002947/5391185‑5391476                                                                                                                                                  |                                                                                                                                                    ccAAAACGCCCGCGAACAGGCCTCGGAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTc                                                                                                                                                   <  2:238975/148‑1 (MQ=255)  cAAAACGCCCGCGAACAGGCCTCGGAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCgg                                                                                                                                                 <  3:42574/149‑1 (MQ=255)   aaaaCGCCCGCGAACAGGCCTCGGAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGc                                                                                                                                                <  1:28332/149‑1 (MQ=255)    aaaCGCCCGCGAACAGGCCTCGGAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCa                                                                                                                                               <  1:108122/149‑1 (MQ=255)          ccGCGAACAGGCCTCGGAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAGTTCGGcatcatcc                                                                                                                                         <  4:175928/149‑3 (MQ=255)          ccGCGAACAGGCCTCGGAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAGTTCGGcatcatcc                                                                                                                                         <  4:28820/149‑3 (MQ=255)          ccGCGAACAGGCCTCGGAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGc                                                                                                                                         >  8:32606/1‑149 (MQ=255)          ccGCGAACAGGCCTCGGAACTGATGCGCTATGCGCGGGACCTGACCATCGACACACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGCGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGc                                                                                                                                         >  1:310846/1‑149 (MQ=255)            gcgAACAGGCCTCGGAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGt                                                                                                                                       <  6:381050/149‑1 (MQ=255)                       tCGGAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCa                                                                                                                            >  1:187386/1‑149 (MQ=255)                          gAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAACCTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCCCGCACcg                                                                                                                         <  2:310846/149‑1 (MQ=255)                               gATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATg                                                                                                                    >  5:133778/1‑149 (MQ=255)                                aTGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATg                                                                                                                    <  2:391542/148‑1 (MQ=255)                                  gggCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCTCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCt                                                                                                                  <  8:58291/148‑1 (MQ=255)                                     cTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTaaa                                                                                                              <  7:52747/149‑1 (MQ=255)                                        tGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATACGTTGAATCACGCACCGCGATGCTTAAAggg                                                                                                           <  6:203126/149‑1 (MQ=255)                                         gcgcGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAgggg                                                                                                          >  8:286547/1‑149 (MQ=255)                                            cgGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCa                                                                                                       >  4:240883/1‑149 (MQ=255)                                                 ccTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGt                                                                                                   >  6:76271/1‑148 (MQ=255)                                                     aCCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAg                                                                                                <  6:96310/147‑1 (MQ=255)                                                        aTCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTaaa                                                                                                              <  1:91014/130‑1 (MQ=255)                                                        aTCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTaaa                                                                                                              >  2:91014/1‑130 (MQ=255)                                                              ggACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGccc                                                                                     <  8:194052/149‑1 (MQ=255)                                                                 cTGATGGAAGACAAGCCCAAACTGATCTGGGCGACGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCatgat                                                                                                                                           >  2:299928/1‑92 (MQ=255)                                                                    aTGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCtttt                                                                                <  7:286547/148‑1 (MQ=255)                                                                      ggAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGAt                                                                                                                     >  1:20183/1‑109 (MQ=255)                                                                      ggAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGAt                                                                                                                     <  2:20183/109‑1 (MQ=255)                                                                      ggAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACg                                                                             >  8:217562/1‑149 (MQ=255)                                                                       gAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGTATGATGCCGAACACGGAATGATGCGTTGAATCACCCACCGCGATGCTAAAAGGGGTAATCGACAGCCGATTGGCCCCTTTTAccc                                                                            >  5:344188/1‑147 (MQ=255)                                                                       gAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACg                                                                             <  6:133778/148‑1 (MQ=255)                                                                       gAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACg                                                                             <  3:240883/148‑1 (MQ=255)                                                                         aGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCca                                                                          >  1:157117/1‑149 (MQ=255)                                                                               gCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGAt                                                                   >  2:115872/1‑149 (MQ=255)                                                                                cccAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATc                                                                  <  5:76271/149‑1 (MQ=255)                                                                                    aaCTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACg                                                              >  3:181843/1‑149 (MQ=255)                                                                                         aTCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGa                                                          <  7:217562/148‑1 (MQ=255)                                                                                             gggCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGAt                                                     >  1:353537/1‑149 (MQ=255)                                                                                             gggCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGAt                                                     >  1:347937/1‑149 (MQ=255)                                                                                              ggcggcGCACTACCTCTGTTCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGTTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAgg                                                        <  1:99148/145‑1 (MQ=255)                                                                                              ggcggcGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGAt                                                     <  7:32606/148‑1 (MQ=255)                                                                                              ggcggcGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATg                                                    <  7:61881/149‑1 (MQ=255)                                                                                                cggcgCACTACCTCTGTGCTCTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCg                                                  >  5:142388/1‑149 (MQ=255)                                                                                                   cgcACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCaa                                                                >  4:382109/1‑132 (MQ=255)                                                                                                   cgcACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCaa                                                                <  3:382109/132‑1 (MQ=255)                                                                                                      aCTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCg                                                  >  3:370542/1‑143 (MQ=255)                                                                                                      aCTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCg                                                  <  4:370542/143‑1 (MQ=255)                                                                                                         aCCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTtct                                          >  3:203766/1‑148 (MQ=255)                                                                                                          ccTCTGTGCTCTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCa                                        <  6:142388/149‑1 (MQ=255)                                                                                                           ctctGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTtct                                          >  3:250092/1‑146 (MQ=255)                                                                                                           ctctGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAg                                       >  3:44276/1‑149 (MQ=255)                                                                                                                   tttGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTc                               >  6:157557/1‑149 (MQ=255)                                                                                                                   tttGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTc                               >  5:128593/1‑149 (MQ=255)                                                                                                                        gCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTc                                                                                                        >  3:98522/1‑72 (MQ=255)                                                                                                                        gCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTc                                                                                                        <  4:98522/72‑1 (MQ=255)                                                                                                                        gCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCAtt                          >  5:348026/1‑149 (MQ=255)                                                                                                                         cAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCAttt                         >  6:276844/1‑149 (MQ=255)                                                                                                                         cAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACCCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCAtt                          >  7:211672/1‑148 (MQ=255)                                                                                                                          aaaGCGCTACTTGATGATGCCTAACTCGTCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTa                        <  6:344188/149‑1 (MQ=255)                                                                                                                          aaaGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTa                        <  6:128593/149‑1 (MQ=255)                                                                                                                            aGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGc                                                                                       <  8:153859/85‑1 (MQ=255)                                                                                                                            aGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGc                                                                                       >  7:153859/1‑85 (MQ=255)                                                                                                                             gcgcTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATgc                     <  1:53399/149‑1 (MQ=255)                                                                                                                              cgcTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATgcg                    >  8:36017/1‑149 (MQ=255)                                                                                                                                   acgcatcatgCCGAGTTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGcgccg                 >  5:357175/8‑147 (MQ=255)                                                                                                                                    tGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGCGCCgcag              >  5:161280/1‑149 (MQ=255)                                                                                                                                      gatgatGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCACCACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGCGCCgcagc             <  1:115872/149‑1 (MQ=255)                                                                                                                                       atcatgCCGAGTTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGCGCcaca               <  4:288477/142‑1 (MQ=255)                                                                                                                                       atcatgCCGAGTTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGCGCcaca               >  3:288477/4‑145 (MQ=255)                                                                                                                                        tgatgCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGCGCCGCAGCATc          <  2:353537/149‑1 (MQ=255)                                                                                                                                        tgatgCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGCGCCGCAGCATc          <  2:347937/149‑1 (MQ=255)                                                                                                                                        tcatgCCGAGTTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGCGCCGCAGCATc          >  3:380018/3‑149 (MQ=255)                                                                                                                                                aaCTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGCGCCGCAGCATCGCCATCaa  <  4:203766/149‑1 (MQ=255)                                                                                                                                                aaCTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGCGCCGCAGCATCGCCATCaa  >  8:48410/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                    CCAAAACGCCCGCGAACAGGCCTCGGAACTGATGGGCTATGCGCGGGACCTGACCATCGACGGACTGATGGAAGACAAGCCCAAACTGATCTGGGCGGCGCACTACCTCTGTGCTTTGGGCAAAGCGCTACTGGATGATGCCGAACTCGGCATGATGCGTTGAATCACGCACCGCGATGCTTAAAGGGGTCATCGTCAGCCGCTTGGCCCCTTTTACGCCA‑CACGATCAACGCCGAGGGATGCGGTCTTTCTCAGCAACCGTCGCATTTATGCGCCGCAGCATCGCCATCAA  >  NC_002947/5391185‑5391476

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