Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,131,860	C→A	6.4%	L72M (CTG→ATG) ‡	PP_4517 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,131,860	0	C	A	6.4%	196.1 / 3.7	62	L72M (CTG→ATG) ‡	PP_4517	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/28);  new base A (2/2);  total (32/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 = 2.78e-01

GACTATGTCGCCGCACACGGCTACCTCACCGCTCTTTCGATCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACC‑TGCGCGGCTGGTGCATCGGCTTCATGGAAGGCGTGTTCCTGCGG  >  NC_002947/5131719‑5132006                                                                                                                                                |                                                                                                                                                    gACTATGTCGCCGCACACGGCTACCTCACCGCTCTTTCGATCATCGCCGAGGATGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCtggtgg                                                                                                                                               <  7:68959/148‑1 (MQ=255) gACTATGTCGCCGCACACGGCTACCTCACCGCTCTTTCGATCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCtggtg                                                                                                                                                <  8:157474/147‑1 (MQ=255) gACTATGTCGCCGCACACGGCTACCTCACCGCTCTTTCGATCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCtggtg                                                                                                                                                <  1:344636/147‑1 (MQ=255) gACTATGTCGCCGCACACGGCTACCTCACCGCTCTTTCGATCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGgc                                                                                                                                              <  1:310957/149‑1 (MQ=255) gACTATGTCGCCGCACACGGCTACCTCACCGCTCTTTCGATCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGgc                                                                                                                                              >  4:165316/1‑149 (MQ=255)   cTATGTCGCCGCACACGGCTACCTCACCGCTCTTTCGATCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGTCCAGCGCACCGAGATCGAGGCCACCAGGGTGgcct                                                                                                                                            <  1:237342/149‑3 (MQ=255)        tCGCCGCACACGGCTACCTCACCGCTCTTTCGATCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAg                                                                                                                                       <  5:169589/149‑1 (MQ=255)          gccgcACACGGCTACCTCACCGCTCTTTCGATCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAgg                                                                                                                                      >  5:102402/1‑148 (MQ=255)               acacGGCTACCTCACCGCTCTTTCGATCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGcacaca                                                                                                                                >  7:153121/1‑149 (MQ=255)                       aCCTCACCGCTCTTTCGATCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGc                                                                                                                        >  4:268187/1‑149 (MQ=255)                          tCACCGCTCTTTCGATCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCACCAg                                                                                                                     >  1:282725/1‑149 (MQ=255)                          tCACCGCTCTTTCGATCAACGCCGAGGACGAGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCACCAg                                                                                                                     >  1:286961/1‑149 (MQ=255)                           cACCGCTATATCGATCAACGCCGAGGACGAGCCAGAGCCCGAATGGAACGACGCCCTCTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGCTGCCACCCTCGAGGCGCTCAAGGCAAACATTGCCCGCCAGc                                                                                                                    >  6:277446/1‑149 (MQ=255)                            aCCGCTCTTTCGATCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCACCAg                                                                                                                     <  2:286961/147‑1 (MQ=255)                            aCCGCTCTTTCGATCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCACCAg                                                                                                                     <  2:282725/147‑1 (MQ=255)                              cGCTCTTTCGATCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGcaca                                                                                                                                  <  7:339606/132‑1 (MQ=255)                              cGCTCTTTCGATCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGcaca                                                                                                                                  >  8:339606/1‑132 (MQ=255)                                 tctTTCGATCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCAGGGTGgcctc                                                                                                                                           >  3:308960/1‑117 (MQ=255)                                 tctTTCGATCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCAGGGTGgcctc                                                                                                                                           <  4:308960/120‑4 (MQ=255)                                   tttCGATCAACGCCGAGGACGTGCCCGAGCGCGAATTGATCGACGCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGcacac                                                                                                                                 <  4:221575/128‑1 (MQ=255)                                   tttCGATCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACACTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGc                                                                                                            >  6:371899/1‑149 (MQ=255)                                   tttCGATCAACGCCGAGGACGTGCCCGAGCGCGAATCGATCGACGCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGcacac                                                                                                                                 >  3:221575/1‑128 (MQ=255)                                         tCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACgaa                                                                                                      <  7:293101/149‑1 (MQ=255)                                           aaCGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCC‑GCACTACGGCCAGCGAAGCCCAGCGCACCGAGATCAAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGcgacg                                                                                                        >  1:107799/1‑146 (MQ=255)                                                    gACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGAcct                                                                                           >  6:371459/1‑149 (MQ=255)                                                        tGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGAcctgccc                                                                                       <  3:165316/149‑1 (MQ=255)                                                         gCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGGCACTAC‑GCCAGCGAAGCCCGGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGAcctgccc                                                                                       >  1:69087/1‑149 (MQ=255)                                                                    aaTGGATCGACGCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGcgacctcga                                                                            <  1:20874/148‑1 (MQ=255)                                                                    aaTGGATCGACGCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGcgacctcga                                                                            >  2:20874/1‑148 (MQ=255)                                                                         aTCGACGCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCGGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGAc                                                                      <  2:69087/149‑1 (MQ=255)                                                                         aTCGACGCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGAc                                                                      <  5:371459/149‑1 (MQ=255)                                                                             aCGCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGAt                                                                  >  3:314980/1‑149 (MQ=255)                                                                              cGCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTAGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGAt                                                                  >  2:18812/1‑148 (MQ=255)                                                                               gCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGAc                                                                      >  7:200472/1‑143 (MQ=255)                                                                               gCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGAc                                                                      <  8:200472/143‑1 (MQ=255)                                                                               gCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGAt                                                                  >  7:136963/1‑147 (MQ=255)                                                                                 ccTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACAGATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGTCCTCGACCTGACCGATGAGc                                                              <  4:314980/149‑1 (MQ=255)                                                                                 ccTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCg                                                                                               >  4:89391/1‑116 (MQ=255)                                                                                 ccTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCg                                                                                               <  3:89391/116‑1 (MQ=255)                                                                                     ttCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGa                                                          <  6:102402/149‑1 (MQ=255)                                                                                                    gCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCAAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACC‑Tgc                                           <  2:107799/149‑1 (MQ=255)                                                                                                    gCACTAC‑GCCAGCGAAGCCCAGCGCACAGAGATCGCGCCCACCCTGGTGGCGCACAACGCACACATTGCCCGCCAGCTGCCCAGCGACGAACAAATCGACCTGCCCTGATACCTCGACCTGACCGATGACCCGGACGCTTCCGAcca                                              >  4:318512/1‑146 (MQ=255)                                                                                                            gCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAg                                                                                                                     <  8:306237/69‑1 (MQ=255)                                                                                                            gCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAg                                                                                                                     >  7:306237/1‑69 (MQ=255)                                                                                                             ccAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACC‑TGCGCGGCTgg                                   <  8:153121/149‑1 (MQ=255)                                                                                                                   aaGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACC‑TGCGCGGCTGGTGCATc                             >  4:355750/1‑149 (MQ=255)                                                                                                                       ccAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACC‑TGCGCGGCTGGTGCATCGGCt                         >  7:293352/1‑149 (MQ=255)                                                                                                                          gcgcACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACC‑TGCGCGGCTGGTGCATCGGCTTCa                      <  8:109530/149‑1 (MQ=255)                                                                                                                              aCCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACC‑TGCGCGGCTGGTGCATCGGCTTCATgg                   >  3:225867/1‑148 (MQ=255)                                                                                                                               ccGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACC‑TGCGCGGCTGGTGCATCGGCTTCATGGaa                 <  8:136963/149‑1 (MQ=255)                                                                                                                                cGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACC‑TGCGCGGCTGGTGCATCGGCTTCATGGAAg                >  8:271494/1‑149 (MQ=255)                                                                                                                                 gagaTCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACC‑TGCGCGGCTGGTGCATCGGCTTCATGGAAgg               >  8:186274/1‑149 (MQ=255)                                                                                                                                  agaTCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACC‑TGCGCGGCTGGTGCATCGGCTTCATGGAAGGc              <  3:268187/149‑1 (MQ=255)                                                                                                                                    ttgagcgCCACCAGGGTGGC‑CTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACC‑TGCGCGGCTGGTGCATCGGCTTCATGGAAGACgtg           >  1:355077/7‑149 (MQ=255)                                                                                                                                    aTCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACC‑TGCGCGGCTGGTGCATCGGCTTCa                      <  3:127333/139‑1 (MQ=255)                                                                                                                                    aTCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACC‑TGCGCGGCTGGTGCATCGGCTTCa                      >  4:127333/1‑139 (MQ=255)                                                                                                                                    aTCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACC‑TGCGCGGCTGGTGCATCGGCTTCATGGAAGGCgt            <  6:322912/149‑1 (MQ=255)                                                                                                                                    aTCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACC‑TGCGCGGCTGGTGCATCGGCTTCATGGAAGGCgt            >  2:221118/1‑149 (MQ=255)                                                                                                                                        atgCCACCCTGGTTGCGCTCAAGGCACACAATTCCCGCCAGCTGGCCATCGACGAAGAATTCGCACTGCCCTGCGACCTCGACCTGAACGATGAGCCGGACGATTCCGACCTTCCCCGGCTGGTGCATCGGCTTCATGGAAGGCgtgt          <  1:221118/146‑1 (MQ=255)                                                                                                                                        aGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACC‑TGCGCGGCTGGTGCATCGGCTTCATGGAAGGCGTGTTc        >  3:93375/1‑149 (MQ=255)                                                                                                                                             aCCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGAc                                                                      <  3:195962/83‑1 (MQ=255)                                                                                                                                             aCCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGAc                                                                      >  4:195962/1‑83 (MQ=255)                                                                                                                                               ccTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACC‑TGCGCGGCTGGTGCATCGGCTTCATGGAAGGCGTGTTCCTGCgg  >  1:82772/1‑148 (MQ=255)                                                                                                                                               ccTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACC‑TGCGCGGCTGGTGCATCGGCTTCATGGAAGGCGTGTTCCTGCgg  >  6:195031/1‑148 (MQ=255)                                                                                                                                               ccTGGTGGCGCACAAGGCACACATTGCCCGCCAGCTCGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACC‑TGCGCGGCTGGTGCATCGGCTTCATGGAAGGCGTGTTCCTGCgg  >  5:200876/1‑148 (MQ=255)                                                                                                                                                |                                                                                                                                                    GACTATGTCGCCGCACACGGCTACCTCACCGCTCTTTCGATCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCC‑GCACTAC‑GCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACC‑TGCGCGGCTGGTGCATCGGCTTCATGGAAGGCGTGTTCCTGCGG  >  NC_002947/5131719‑5132006

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