Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,226,666	T→C	6.8%	I78V (ATC→GTC)	PP_5623 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,226,666	0	T	C	6.8%	167.9 / 2.4	59	I78V (ATC→GTC)	PP_5623	membrane protein
Reads supporting (aligned to +/- strand):  ref base T (26/29);  new base C (2/2);  total (28/31)
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.40e-01

AAAGTACGGAACGCCCAAGCTCCTGCTCGTAGCAGATAAGCAACTGCCAAAGCAAAGATCAGCAGACCTATACATCCCAGCCATTTGCCTGCCTTCTCTCCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGACTCTTCAAGCGCAACCTCGCTGTCCTTAGACCAATCACTCGACTGGACAGACAATCCCATATCGCGTAT  >  NC_002947/4226518‑4226810                                                                                                                                                       |                                                                                                                                                  aaaGTACGGAACGCCCAAGCTCCTGCTCGTAGCAGATAAGCAACTGCCAAAGCAAAGATCAGCAGACCTATACATCCCAGCCATTTGCCTGCCTTCTCTCCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGAt                                                                                                                                                   >  1:189613/1‑149 (MQ=255)        ggAACGCCCAAGCTCCTGCTCGTAGCAGATAAGCAACTGCCAAAGCAAAGATCAGCAGACCTATACATCCCAGCCATTTGCCTGCCTTCTCTCCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATcgccgc                                                                                                                                             >  2:330925/1‑148 (MQ=255)             gCCCAAGCTCCTGCTCGTAGCAGATAAGCAACTGCCAAAGCAAAGATCAGCAGACCTATACATCCCAGCCATTTGCCTGCCTTCTCTCCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGt                                                                                                                                       >  4:173188/1‑148 (MQ=255)                    cTCGTGCTCGTAGCAGATAAGCAACTGCCATAGCTATGATCAGCAGACCTATACATCGCAGCCATTTGCCTGCCTTCTGTCCCCACAACGATGGCAGCAGGGAAACCAAAACACGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTATATcccc                                                                                                                               <  6:235185/149‑1 (MQ=255)                       cTGCTCGTAGCAGATAAGCAACTGCCAAAGCAAAGATCAGCAGACCTATACATCCCAGCCATTTGCCTGCCTTCTCTCCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGcc                                                                                                                            >  4:272399/1‑149 (MQ=255)                          cTCGTAGCAGATAAGCAACTGCCAAAGCAAAGATCAGCAGACCTATACATCCCAGCCATTTGCCTGCCTTCTCTCCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCt                                                                                                                           <  3:312309/147‑1 (MQ=255)                                  aGATAAGCAACTGCCAAAGCAAAGATCAGCAGACCTATACATCCCAGCCATTTGCCTGCCTTCTCTCCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGc                                                                                                                 <  1:330925/149‑1 (MQ=255)                                      aaGCAACTGCCAAAGCAAAGATCAGCAGCCCTATACATCCCAGCCATTTGCCTGCCTTCTCTCCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATACCCGCCTTAGCATCTGCTCCt                                                                                                             >  8:329302/1‑149 (MQ=255)                                           aCTGCCAAAGCAAAGATCAGCAGACCTATACATCCCAGCCATTTGCCTGCCTTCTCTCCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCg                                                                                                        <  4:115484/149‑1 (MQ=255)                                                 aaagcaaagATCAGCAGACCTATACATCCCAGCCATTTGCCTGCCTTCTCTCCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGtt                                                                                                   >  4:77423/1‑148 (MQ=255)                                                    gcaaagATCAGCAGACCTATACATCCCAGCCATTTGCCTGCCTTCTCTCCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAAt                                                                                               >  1:188103/1‑149 (MQ=255)                                                        agATAAGCAGACCTATACATACCAGCCATTTGACTGCCTTCTCTACCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCTCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCg                                                                                           <  4:229592/149‑1 (MQ=255)                                                          aTCAGCAGACCTATACATCCCAGCCATTTGCCTGCCTTCTCTCCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGc                                                                                         <  3:173188/149‑1 (MQ=255)                                                                agACCTATACATCCCAGCCATTTGCCTGCCTTCTCTCCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCaaa                                                                                   <  2:188103/149‑1 (MQ=255)                                                                  aCCTATACATCCCAGCCATTTGCCTGCCTTCTCTCCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGAATATGCTCCTTTTCGGCGTTCAATCTCGGCATCaaa                                                                                   >  8:124890/1‑147 (MQ=255)                                                                      ataCATCCCAGCCATTTGCCTGCCTTCTCTCCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCGAATAGGGATCGACGGATAAGTAAATAACCGCCTTAGGATCATCTACTTTTAGGCGTTCAATCTAGGCATCAAaaag                                                                                >  1:282157/1‑149 (MQ=255)                                                                          aTCCCAGCCATTTGCCTGCCTTCTCTCCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGcgccgc                                                                         >  8:159126/1‑149 (MQ=255)                                                                          aTCCCAGCCATTTGCCTGCCTTCTCTCCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGACCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCg                                                                                           <  6:280850/131‑1 (MQ=255)                                                                          aTCCCAGCCATTTGCCTGCCTTCTCTCCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGACCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCg                                                                                           >  5:280850/1‑131 (MQ=255)                                                                          aTCCCAGCCATTTGCCTGCCTTCTCTCCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGACCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCg                                                                                           <  4:267696/131‑1 (MQ=255)                                                                          aTCCCAGCCATTTGCCTGCCTTCTCTCCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGACCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCg                                                                                           >  3:267696/1‑131 (MQ=255)                                                                             caaGCAATTTTCCTGCCTTCTCTCCCCAGAACGATGGCTGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGa                                                                       <  2:282157/146‑1 (MQ=255)                                                                              cAGCCATTTGCCTGCCTTCTCTCCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGAct                                                                     >  7:152493/1‑149 (MQ=255)                                                                              cAGCCATTTGCCTGCCTTCTCTCCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGAct                                                                     <  8:152493/149‑1 (MQ=255)                                                                                gCCATTTGCCTGCCTTCTCTCCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGActct                                                                   <  3:272399/149‑1 (MQ=255)                                                                                   aTTTGCCTTCCTTCTCTCCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAaga                                                                                 <  3:204026/132‑1 (MQ=255)                                                                                   aTTTGCCTGCCTTCTCTCCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCg                                                                                                        >  6:311972/1‑109 (MQ=255)                                                                                   aTTTGCCTGCCTTCTCTCCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCg                                                                                                        <  5:311972/109‑1 (MQ=255)                                                                                   aTTTGCCTGCCTTCTCTCCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAaga                                                                                 >  4:204026/1‑132 (MQ=255)                                                                                       gcctgcctTCTCTCCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGACTCTTCAAgcg                                                            >  7:159878/1‑149 (MQ=255)                                                                                         ctgcctTCTCTCCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCAGCGACTCTTCAAgcgc                                                           >  3:50247/1‑148 (MQ=255)                                                                                          tgcctTCTCTCCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGACTCTTCAAGCGCa                                                          >  2:93841/1‑148 (MQ=255)                                                                                           gcctTCTCTCCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGACTCTTCAAGCGCAAc                                                        <  7:329302/149‑1 (MQ=255)                                                                                                   tcCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGACTCTTCAAGCGCAACCTCGCt                                                  <  8:159878/147‑1 (MQ=255)                                                                                                   tcCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGACTCTTCAAGCGCAACCTCGCTGt                                                >  1:347083/1‑149 (MQ=255)                                                                                                    ccccACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGACTCTTCAAgc                                                             >  2:307883/1‑135 (MQ=255)                                                                                                    ccccACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGACTCTTCAAgc                                                             <  1:307883/135‑1 (MQ=255)                                                                                                     cccACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGAATCTGCTCCTTTTCGGCGTTCAATCACGGCATCAAAGAGGCGCCGCGACTCTTCAAGCGCAACCTCGCTGTcc                                              >  3:245374/1‑149 (MQ=255)                                                                                                     cccACAAAGATGGCAGCAGGGAAACCAATACAGGAATTAGCG‑‑CAGTGATCGCCTCA‑AGGTAAATCCCCGCCTTAGCATCTGATCCTTATCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGACTCTTCAAGCGCAACCTCGCTGTc                                               <  7:124890/148‑1 (MQ=255)                                                                                                      ccACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGACTCTTCAAGCGCAACCTCGCTGTCCt                                             <  3:77423/149‑1 (MQ=255)                                                                                                        acaACGATGTCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGACTCTTCAAGCGCAACCTCGCTGTCCTTa                                           <  1:93841/149‑1 (MQ=255)                                                                                                        acaACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGACTCTTCAAGCGCAACCTCGCTGTCCt                                             >  7:280475/1‑147 (MQ=255)                                                                                                        acaACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGACTCTTCAAGCGCAACCTCGCTGTCCt                                             >  7:30836/1‑147 (MQ=255)                                                                                                        acaACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGACTCTTCAAGCGCAACCTCGCTGTCCt                                             >  5:2346/1‑147 (MQ=255)                                                                                                        acaACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGACTCTTCAAGCGCAACCTCGCTGTCCTTa                                           <  6:2346/149‑1 (MQ=255)                                                                                                        acaACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGACTCTTCAAGCGCAACCTCGCTGTCCTTa                                           <  8:30836/149‑1 (MQ=255)                                                                                                                    gctgggAAACCATAACAGGAATTAGCT‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCATTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGACTCTTCAAGCGCAACCTCTCTGTCCTTAGACCaatcaatc                               <  4:245374/146‑1 (MQ=255)                                                                                                                     cagGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGACTCTTCAAGCGCAACCTCGCTGTCCTTa                                           <  7:244848/136‑1 (MQ=255)                                                                                                                     cagGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGACTCTTCAAGCGCAACCTCGCTGTCCTTa                                           >  8:244848/1‑136 (MQ=255)                                                                                                                       gggAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGACTCTTCAAGCGCAACCTCGCTGTCCTTAGACCAATCACTCGAc                            >  4:206952/1‑149 (MQ=255)                                                                                                                           aaCCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGACTCTTCAAGCGCAACCTCGCTGTCCTTAGACCAATCACTCGACTGga                        <  7:159126/149‑1 (MQ=255)                                                                                                                              cAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGACTCTTCAAGCGCAACCTCGCTGTCCTTAGACCAATCACTCGACTGga                        >  7:188753/1‑146 (MQ=255)                                                                                                                               aaaaCAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGACTCTTCAAGCGCAACCTCGCTGTCCTTa                                           <  2:180487/126‑1 (MQ=255)                                                                                                                               aaaaCAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGACTCTTCAAGCGCAACCTCGCTGTCCTTa                                           >  1:180487/1‑126 (MQ=255)                                                                                                                                    aGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGACTCTTCAAGCGCAACCTCGCTGTCCTTAGACCAATCACTCGACTGGACAGACAATc               <  2:189613/149‑1 (MQ=255)                                                                                                                                       aaTTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGACTCTTCAAGCGCAACCTCGCTGTCCTTAGACCAAt                                    <  7:1541/125‑1 (MQ=255)                                                                                                                                            gcg‑‑cAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGACTCTTCAAGCGCAACCTCGCTGTCCTTAGACCAATCACTCGACTGGACAGACAATCCCATATc        <  8:280475/148‑1 (MQ=255)                                                                                                                                             cg‑‑cAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGACTCTTCAAGCGCAACCTCGCTGTCCTTAGACCAATCACTCGACTGGACAGACAATCCCATATcgc      <  4:50247/149‑1 (MQ=255)                                                                                                                                             cg‑‑cAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGACTCTTCAAGCGCAACCTCGCTGTCCTTAGACCAATCACTCGACTGGACAGACAATCCCATATcgc      <  2:347083/149‑1 (MQ=255)                                                                                                                                                   gTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGACTCTTCAAGCGCAACCTCGCTGTCCTTAGACCAATCACTCGACTGGACAGACAATCCCATATCGCGTAt  >  4:55975/1‑149 (MQ=255)                                                                                                                                                       |                                                                                                                                                  AAAGTACGGAACGCCCAAGCTCCTGCTCGTAGCAGATAAGCAACTGCCAAAGCAAAGATCAGCAGACCTATACATCCCAGCCATTTGCCTGCCTTCTCTCCCCACAACGATGGCAGCAGGGAAACCAAAACAGGAATTAGCG‑‑CAGTGATCGCCGCA‑AGGTAAATCCCCGCCTTAGCATCTGCTCCTTTTCGGCGTTCAATCTCGGCATCAAAGAGGCGCCGCGACTCTTCAAGCGCAACCTCGCTGTCCTTAGACCAATCACTCGACTGGACAGACAATCCCATATCGCGTAT  >  NC_002947/4226518‑4226810

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