Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,839,503	C→A	6.3%	P49H (CCC→CAC)	PP_4254 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,839,503	0	C	A	6.3%	199.8 / 3.2	63	P49H (CCC→CAC)	PP_4254	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (31/28);  new base A (2/2);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.24e-01

AGTGGGCTTCGCTCCGAAAGGACCGATGCAGACGCGGCGCATGGCCTGTCGCCGCCCCAGA‑AAAAGCTTGACCGATCGATCAGAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTGGTTATCAATACTGGCGCCGATTTTTCGACCAACAAGAACATCAAAACCGTGGAACCTTAGAATGAGCACAGCAATCA  >  NC_002947/4839357‑4839650                                                                                                                                                    |                                                                                                                                                     aGTGGGCTTCGCTCCGAAAGGACCGATGCAGACGCGGCGCATGGCCTGTCGCCGCCCCAGA‑AAAAGCTTGACCGATCGATCAGAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCt                                                                                                                                                    >  8:234178/1‑149 (MQ=255)     ggCTTCGCTCCGAAAGGACCGATGCAGACGCGGCGCATGGCCTGTCGCCGCCCCAGA‑AAAAGCTTGACCGATCGATCAGAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCgg                                                                                                                                                <  3:286647/149‑1 (MQ=255)          cGCTCCGAAAGGACCGATGCAGACGCGGCGCAAGGCCTGTCGCCGCCCCAGAGAAAAGCTTGACCGATCGATCAGAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACa                                                                                                                                            >  1:80347/1‑149 (MQ=255)             tCCGAAAGGACCGATGCAGACGCGGCGCATGGCCTGTCGCCGCCCCAGA‑AAAAGCTTGACCGATCGATCAGAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAg                                                                                                                                        >  1:17229/1‑149 (MQ=255)               cGAAAGGACCGATGCAGACGCGGCGCATGGCCTGTCGCCGCCCCAGA‑AAAAGCTTGACCGATCGATCAGAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGt                                                                                                                                       <  3:211087/148‑1 (MQ=255)                  aaGGACCGATGCAGACGCGGCGCATGGCCTGTCGCCGCCCCAGA‑AAAAGCTTGACCGATCGATCAGAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAc                                                                                                                                   >  1:170043/1‑149 (MQ=255)                         gATGCAGACGCGGCGCATGGCCTGTCGCCGCCCCAGA‑AAAAGCTTGACCGATCGATCAGAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTg                                                                                                                           <  8:19661/149‑1 (MQ=255)                          aTGCAGACGCGGCGCATGGCCTGTCGCCGCCCCAGA‑AAAAGCTTGACCGATCGATCAGAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCACTGCGGGACATCAGTGAAC‑GGGTCTgg                                                                                                                          >  3:59681/1‑149 (MQ=255)                           tGCAGACGCGGCGCATGGCCTGTCGCCGCCCCAGA‑AAAAGCTTGACCGATCGATCAGAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGTTCACTGATGTCCCGCAGGGATC‑GGGTCTggg                                                                                                                         <  8:43902/149‑2 (MQ=37)                             cAGACGCGGCGCATGGCCTGTCGCCGCCCCAGA‑AAAAGCTTGACCGATCGATCAGAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGaca                                                                                                                       >  3:222342/1‑149 (MQ=255)                                aCGCGGCGCATGGCCTGTCGCCGCCCCAGA‑AAAAGCTTGACCGATCGATCAGAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACgg                                                                                                                    <  7:234178/149‑1 (MQ=255)                                  gcggcgCATGGCCTGTCGCCGCCCCAGA‑AAAAGCTTGACCGATCGATCAGAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGt                                                                                                                  >  2:149393/1‑149 (MQ=255)                                       gcATGGCCTGTCGCCGCCCCAGA‑AAAAGCTTGACCGATCGATCAGAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCg                                                                                                             <  4:222342/149‑1 (MQ=255)                                         aTGGCCTGTCGCCGCCCCAGA‑AAAAGTTTGACCGATCGATCAGAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGTTCACTGATGTCCCGCAGGGATC‑GGGTCTgg                                                                                                                          >  2:284842/1‑134 (MQ=25)                                         aTGGCCTGTCGCCGCCCCAGA‑AAAAGTTTGACCGATCGATCAGAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGTTCACTGATGTCCCGCAGGGATC‑GGGTCTgg                                                                                                                          <  1:284842/134‑1 (MQ=25)                                              cTGTCGCCGCCCCAGA‑AAAAGCTTGACCGATCGATCAGAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTTCCa                                                                                                        >  7:134707/1‑147 (MQ=255)                                              cTGTCGCCGCCCCAGA‑AAAAGCTTGACCGATCGATCAGAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCa                                                                                                        >  8:58919/1‑147 (MQ=255)                                                 tCGCCGCCCCAGA‑AAAAGCTTGACCGATCGATCAGAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAgt                                                                                                    <  6:63901/148‑1 (MQ=255)                                                        cccAGA‑AAAAGCTTGACCGATCGATCAGAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAgttgcgttgc                                                                                            >  7:19819/1‑149 (MQ=255)                                                        cccAGA‑AAAAGCTTGACCGATCGATCAGAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAgttg                                                                                                  >  6:217790/1‑143 (MQ=255)                                                        cccAGA‑AAAAGCTTGACCGATCGATCAGAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAgttg                                                                                                  <  5:217790/143‑1 (MQ=255)                                                            gA‑AAAAGCTTGACCGATCGATCAGAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCATGCGGGACATCAGTGAAC‑GGGTATGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGt                                                                                         >  3:93252/1‑148 (MQ=255)                                                                       gACCGATCGATCAGAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCt                                                                                <  2:17229/147‑1 (MQ=255)                                                                       gACCGATCGATCAGAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCt                                                                                <  7:58919/147‑1 (MQ=255)                                                                       gACCGATCGATCAGAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCt                                                                                <  5:35182/147‑1 (MQ=255)                                                                       gACCGATCGATCAGAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTg                                                                               <  1:149393/148‑1 (MQ=255)                                                                         ccGATCGATCAGAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTGGt                                                                             >  7:58198/1‑148 (MQ=255)                                                                         ccGATCGATCAGAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTGGt                                                                             >  3:64425/1‑148 (MQ=255)                                                                            atcgatcAGAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAgt                                                                                                    >  1:15017/1‑122 (MQ=255)                                                                            atcgatcAGAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAgt                                                                                                    <  2:15017/122‑1 (MQ=255)                                                                             tcgatcATAAAAAGATGAAAAACGGTACACATTACAAAAATTTATTTGTGTACAATCGCCCAGACCCTATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTGGTTATCa                                                                        <  3:267107/149‑1 (MQ=255)                                                                                atcAGAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTGGTTATCAAt                                                                      >  1:36845/1‑148 (MQ=255)                                                                                    gAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACggg                                                                                                                   >  7:128466/1‑99 (MQ=255)                                                                                    gAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACggg                                                                                                                   <  8:128466/99‑1 (MQ=255)                                                                                              aaaaaCGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTGGTTATCAATACt                                                                   >  3:161376/1‑137 (MQ=255)                                                                                              aaaaaCGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTGGTTATCAATACt                                                                   <  4:161376/137‑1 (MQ=255)                                                                                               aaaaCGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGTCACGGCTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTGGTTATCTATTGTGGCGCCGAttttt                                                      <  2:80347/149‑1 (MQ=255)                                                                                                   cGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTGGTTATCAATACTGGCGCCGATTTTTCGAc                                                  <  4:140428/149‑1 (MQ=255)                                                                                                      tACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTGGTTATCAATACTGGCGCCGATTTTTCGACcaa                                               <  8:58198/149‑1 (MQ=255)                                                                                                          caTTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTGGTTATCAATACTGGCGCCGATTTTTCGACcaac                                              >  8:21981/1‑146 (MQ=255)                                                                                                            ttACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTGGTTATCAATACTGGCGCCGATTTTTCGACCAACAAGa                                          >  5:75986/1‑148 (MQ=255)                                                                                                             tACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTGGTTATCAATACTGGCGCCGATTTTTCGACCAACAAGaa                                         >  6:148921/1‑148 (MQ=255)                                                                                                               cAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTGGTTATCAATACTGGCGCCGATTTTTCGACCAACAAGAACAt                                      <  4:64425/149‑1 (MQ=255)                                                                                                                    atttatttGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTGGTTATCAATACTGGCGCCGATTTTTCGa                                                   >  4:106105/1‑131 (MQ=255)                                                                                                                    atttatttGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTGGTTATCAATACTGGCGCCGATTTTTCGa                                                   <  3:106105/131‑1 (MQ=255)                                                                                                                        atttGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTGGTTATCAATACTGGCGCCGATTTTTCGACCAACAAGAACATCAAAACCGt                             >  1:152817/1‑149 (MQ=255)                                                                                                                         tttGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTGGTTATCAATACTGGCGCCGATTTTTCGACCAACAAGAACATCAAAACCGt                             <  2:170043/148‑1 (MQ=255)                                                                                                                          ttGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTGGTTATCAATACTGGCGCCGATTTTTCGACCAACAAGAACATCAAAACCGTg                            <  7:21981/148‑1 (MQ=255)                                                                                                                           tgtgTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTGGTTATCAATACTGGCGCCGATTTTTCGACCAACAAGAACATCAAAACCGTgg                           >  5:154586/1‑148 (MQ=255)                                                                                                                            gtgtACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTGGTTATCAATACTGGCGCCGATTTTTCGACCAACAAGAACATCAAAACCGTGGaa                         >  7:270156/1‑149 (MQ=255)                                                                                                                               tACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTGGTTATCAATACTGGCGGCGATTTTTCGACCAACAAGAACATCAAAACCGTGGAACCt                      >  4:281512/1‑149 (MQ=255)                                                                                                                               tACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTGGTTATCAATACTGGCGCCGATTTTTCGACCAACAAGAACATCAAAACCGTGGAACCt                      >  5:218824/1‑149 (MQ=255)                                                                                                                               tACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTGGTTATCAATACTGGCGCCGATTTTTCGACCAACAAGAACATCAAAACCGTGGAACCt                      >  4:265844/1‑149 (MQ=255)                                                                                                                                    tCGCCCAGACTCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTg                                                                               >  4:54207/1‑87 (MQ=255)                                                                                                                                    tCGCCCAGACTCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTg                                                                               <  3:54207/87‑1 (MQ=255)                                                                                                                                    tCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTGGTTATCAATACTGGCGCCGATTTTTCGACCAACAAGAACATCAAAACCGTGGAACCTTAGaa                 >  2:278320/1‑149 (MQ=255)                                                                                                                                      gCCCAGACCCGATCCCTGCGGGACATCAGTGAACGGGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTGGTTATCAATACTGGCGCCGATTTTTCGACCAACAAGAACATCAAAACCGTGGAACCTTAGAAt                >  1:210798/1‑149 (MQ=255)                                                                                                                                      gCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTGGTTATCAATACTGGCGCCGATTTTTCGACCAACAAGAACATCAAAACCGTGGAACCTTAGAATg               <  6:75986/149‑1 (MQ=255)                                                                                                                                       cccAGAGCCGATCCCTGCGGGACATAAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCTGTAACCTTGCTGGTTATCAATACTGGCGCCGATTTTTCGACaaaca                                             <  5:329685/118‑1 (MQ=255)                                                                                                                                         cAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTGGTTATCAATACTGGCGCCGATTTTTCGACCAACAAGAACATCAAAACCGTGGAACCTTAGAATGAGc            <  8:270156/149‑1 (MQ=255)                                                                                                                                            aCCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTGGTTATCAATACTGGCGCCGATTTTTCGACCAACAAGAACATCAAAACCGTGGAACCTTAGAATGAGcaca         <  6:218824/149‑1 (MQ=255)                                                                                                                                              ccGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTGGTTATCAATACTGGCGCCGATTTTTCGACCAACAAGAACATCAAAACCGTGGAACCTTAGAATGAGCACAg        >  5:101631/1‑148 (MQ=255)                                                                                                                                                 aTCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTGGTTATCAATACTGGCGCCGATTTTTCGACCAACAAGAACATCAAAACCGTGGAACCTTAGAATGAGCACAGCAAt    <  2:36845/149‑1 (MQ=255)                                                                                                                                                    ccTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTGGTTATCAATACTGGCGCCGATTTTTCGACCAACAAGAACATCAAAACCGTGGAACCTTAGAATGAGCACAGCAATCa  <  8:19819/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     AGTGGGCTTCGCTCCGAAAGGACCGATGCAGACGCGGCGCATGGCCTGTCGCCGCCCCAGA‑AAAAGCTTGACCGATCGATCAGAAAAAGATGAAAAACGGTACACATTACAAATATTTATTTGTGTACAATCGCCCAGACCCGATCCCTGCGGGACATCAGTGAAC‑GGGTCTGGACACGGGTCGGCGTACCAAAGTTGCGTTGCGGTAACCTTGCTGGTTATCAATACTGGCGCCGATTTTTCGACCAACAAGAACATCAAAACCGTGGAACCTTAGAATGAGCACAGCAATCA  >  NC_002947/4839357‑4839650

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