Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,073,610	T→C	7.6%	intergenic (‑148/+476)	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,610	0	T	C	7.6%	176.7 / 4.8	66	intergenic (‑148/+476)	mdtA/PP_3586	multidrug transporter membrane fusion protein/transposase
Reads supporting (aligned to +/- strand):  ref base T (26/35);  new base C (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 = 5.72e-01

AACGAACCTTGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGAAAAAAACAGGCGCGCGGG  >  NC_002947/4073469‑4073756                                                                                                                                               |                                                                                                                                                   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                                                                                                                              <  7:318563/128‑1 (MQ=255)                                     gAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGgtgt                                                                                                                              >  8:318563/1‑128 (MQ=255)                                       aCGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCt                                                                                                        >  3:289112/1‑148 (MQ=255)                                              gCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACa                                                                                                                                    >  8:321235/1‑113 (MQ=255)                                              gCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACa                                                                                                                                    <  7:321235/113‑1 (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                                                                                          <  8:239217/143‑1 (MQ=255)                                                          ttttAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGc                                                                                          >  7:239217/1‑143 (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                                                                                    >  6:193532/1‑137 (MQ=255)                                                                      gCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTTGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTc                                                                                    <  5:193532/137‑1 (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)                                                                                                                                               |                                                                                                                                                   AACGAACCTTGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTT‑AAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGAAAAAAACAGGCGCGCGGG  >  NC_002947/4073469‑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: ‑