Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,073,609	G→A	7.6%	intergenic (‑147/+477)	mdtA ← / ← PP_3586	multidrug transporter membrane fusion protein/transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,073,609	0	G	A	7.6%	187.1 / 4.0	66	intergenic (‑147/+477)	mdtA/PP_3586	multidrug transporter membrane fusion protein/transposase
Reads supporting (aligned to +/- strand):  ref base G (25/35);  new base A (2/3);  total (28/38)
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.49e-01

ATGGGTTGAACGAACCTTGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGAAAAAAACAGGCGCGCGGG  >  NC_002947/4073461‑4073756                                                                                                                                                      |                                                                                                                                                    aTGGGTTGAACGAACCTTGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACg                                                                                                                                                     >  7:248487/1‑149 (MQ=255) aTGGGTTGAACGAACCTTGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACg                                                                                                                                                     >  3:7337/1‑149 (MQ=255)         aacgaacCTTGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCgg                                                                                                                                             >  3:306833/1‑149 (MQ=255)          acgaccCTTGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGCGCCggg                                                                                                                                            <  1:252782/149‑1 (MQ=255)          acgaacCTTGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCggg                                                                                                                                            <  8:231697/149‑1 (MQ=255)           cgaacCTTGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACACGTTCCGGGc                                                                                                                                           <  4:77147/149‑1 (MQ=255)            gcacCTTGTTCGGAATGACGACTAGACGAGGAGGAACCATAAGCACGGCCATCGTTTAAGCAACTGGCCGTTACCAGAGCTTTACCTTTCCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGcc                                                                                                                                          <  4:331109/147‑1 (MQ=255)            gaacCTTGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGcc                                                                                                                                          <  7:79841/149‑1 (MQ=255)              acCTTGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGt                                                                                                                                        <  2:304234/149‑1 (MQ=255)                   gTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACACGTTCCGGGCCGttttc                                                                                                                                    <  3:360633/148‑4 (MQ=255)                    ttCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAgg                                                                                                                                  <  4:7337/149‑1 (MQ=255)                        gAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGCATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGgtgt                                                                                                                              <  1:285361/149‑1 (MQ=255)                             acgacTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAgg                                                                                                                                  >  1:216895/1‑140 (MQ=255)                             acgacTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAgg                                                                                                                                  <  2:216895/140‑1 (MQ=255)                                    gaggagGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCg                                                                                                                   <  1:49796/148‑1 (MQ=255)                                        aggaTGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGc                                                                                                              <  4:165337/149‑1 (MQ=255)                                            tGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCTGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGGTTGGAGCGGTGGCGGGc                                                                                                          >  5:252781/1‑149 (MQ=255)                                            tGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGc                                                                                                          >  5:258322/1‑149 (MQ=255)                                             gAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGgtgt                                                                                                                              >  8:318563/1‑128 (MQ=255)                                             gAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGgtgt                                                                                                                              <  7:318563/128‑1 (MQ=255)                                               aCGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCt                                                                                                        >  3:289112/1‑148 (MQ=255)                                                      gCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACa                                                                                                                                    <  7:321235/113‑1 (MQ=255)                                                      gCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACa                                                                                                                                    >  8:321235/1‑113 (MQ=255)                                                      gCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTc                                                                                                <  7:254296/149‑1 (MQ=255)                                                              aTCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACACGTTCCGGGCCGtttt                                                                                                                                     <  6:16913/104‑3 (MQ=255)                                                              aTCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACACGTTCCGGGCCGtttt                                                                                                                                     >  5:16913/1‑102 (MQ=255)                                                                  ttttAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTCAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGt                                                                                     >  1:311042/1‑149 (MQ=255)                                                                  ttttAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGc                                                                                          >  7:239217/1‑143 (MQ=255)                                                                  ttttAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGc                                                                                          <  8:239217/143‑1 (MQ=255)                                                                  ttttAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTc                                                                                    >  3:29127/1‑149 (MQ=255)                                                                     tAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTg                                                                                  <  4:29127/148‑1 (MQ=255)                                                                     tAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTg                                                                                  <  4:306833/148‑1 (MQ=255)                                                                       gCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGAGGGCCTGAAGGTTCGGTTGCCTGGTCTGGcc                                                                               >  6:93572/1‑149 (MQ=255)                                                                       gCAAAGGGCCTTTACAATCTGTTTCCCGTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGGCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGcc                                                                               <  7:334799/149‑1 (MQ=255)                                                                        cAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTg                                                                                                       >  1:123447/1‑124 (MQ=255)                                                                        cAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTg                                                                                                       <  2:123447/124‑1 (MQ=255)                                                                              gCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTTGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTc                                                                                    <  5:193532/137‑1 (MQ=255)                                                                              gCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTTGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTc                                                                                    >  6:193532/1‑137 (MQ=255)                                                                              gCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTcgc                                                                        <  2:310581/149‑1 (MQ=255)                                                                                cTTTACCAGAGGTTTACGTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCgg                                                                      <  7:86329/149‑1 (MQ=255)                                                                                   tACCAGAGGTTTACCTTTGCATGACGTTGCCATGAGGCTGAATTTT‑AAATGAAATCGGCCCGGAAAGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGt                                                                    <  6:243451/148‑1 (MQ=255)                                                                                      cAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTgg                                                                                 <  1:45237/132‑1 (MQ=255)                                                                                       agagGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAAcc                                                               >  1:53524/1‑149 (MQ=255)                                                                                       agagGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAAcc                                                               >  7:214144/1‑149 (MQ=255)                                                                                           gTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCt                                                           >  2:228493/1‑149 (MQ=255)                                                                                               aCCTTTGCCTGACGTTGCCAAGCGGCGTAATTTT‑AAATTAAAAGGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCCGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCt                                                        <  6:252781/148‑1 (MQ=255)                                                                                               aCCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCCGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCt                                                        <  6:258322/148‑1 (MQ=255)                                                                                                ccTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTAt                                                      <  8:248487/149‑1 (MQ=255)                                                                                                 cTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTAt                                                      >  3:285256/1‑148 (MQ=255)                                                                                                 cTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATc                                                     <  1:228493/149‑1 (MQ=255)                                                                                                      ccTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATcaaac                                                 >  2:323589/1‑148 (MQ=255)                                                                                                                 cAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATcaaa                                                  <  8:82029/136‑1 (MQ=255)                                                                                                                 cAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATcaaa                                                  >  7:82029/1‑136 (MQ=255)                                                                                                                 cAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATaa                                        >  4:333088/1‑146 (MQ=255)                                                                                                                 cAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATaa                                        <  3:333088/146‑1 (MQ=255)                                                                                                                      ggCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTGGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCGTCGATCAAACTAAAGATAACCCACTGt                                >  8:356840/1‑149 (MQ=255)                                                                                                                         ggAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGc                                                                                                              <  2:357232/68‑1 (MQ=255)                                                                                                                         ggAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGc                                                                                                              >  1:357232/1‑68 (MQ=255)                                                                                                                         ggAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGttct                             >  4:236483/1‑149 (MQ=255)                                                                                                                           aaTTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGttcttt                           <  4:285256/149‑1 (MQ=255)                                                                                                                               tt‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTa                        <  5:93572/148‑1 (MQ=255)                                                                                                                                t‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAg                       <  4:289112/148‑1 (MQ=255)                                                                                                                                      gAAAACGGCCCGGAACTTGTTCAGGTCCGTTACACGGTGTTGCTTGGAGCGGTGGCGGCCCTGAAGGTTCGGTTTCATGTTCAGGACCTTTCGCGGGTGAACCAGCTTCTAGCAAACAACAGATAACCCAATGTTCTTTATTGCGaaa                  >  4:263705/1‑148 (MQ=255)                                                                                                                                      gAAAACCGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACaaaa                                             <  6:288033/121‑1 (MQ=255)                                                                                                                                       gtaaCGGCCCGGAACACGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGaaaa                 >  3:63872/3‑148 (MQ=255)                                                                                                                                       aaaaCGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGaaaa                 >  4:326713/1‑148 (MQ=255)                                                                                                                                             gCCCGGAACGTGTTCCGGGCAGATACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGACCTTTCGCGCGTCAACCCGCATCTATCAAACAAAAGATAAACCACTGTTCTTACTTCCGAAAAAAACAgg          >  4:32037/1‑149 (MQ=37)                                                                                                                                                     cGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGGGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGAAAAAAACAGGCGCGCggg  <  2:311042/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                    ATGGGTTGAACGAACCTTGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGAAAAAAACAGGCGCGCGGG  >  NC_002947/4073461‑4073756

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑