Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,921,589	C→G	11.7%	A666A (GCC→GCG) ‡	PP_4330 →	toxin‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,921,589	0	C	G	11.7%	180.8 / 11.2	66	A666A (GCC→GCG) ‡	PP_4330	toxin‑like protein
Reads supporting (aligned to +/- strand):  ref base C (29/29);  new base G (3/5);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGACAAAGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACC  >  NC_002947/4921441‑4921737                                                                                                                                                     |                                                                                                                                                     gCGACAAAGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAATGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACACAAGGCTATTGGCAGCCCGGGcc                                                                                                                                                      >  1:222534/1‑149 (MQ=255)   gACAAAGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGccc                                                                                                                                                     >  4:167908/1‑148 (MQ=255)      aaaGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAgg                                                                                                                                                 >  5:192076/1‑149 (MQ=255)      aaaGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAgg                                                                                                                                                 >  1:49612/1‑149 (MQ=255)       aaGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAggg                                                                                                                                                <  7:192960/149‑1 (MQ=255)          tCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCTGGGGCCCGGGCTg                                                                                                                                             <  1:9142/149‑1 (MQ=255)            aTTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGc                                                                                                                                            <  6:23639/148‑1 (MQ=255)             ttAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAg                                                                                                                                          <  1:22321/149‑1 (MQ=255)             ttAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAg                                                                                                                                          <  5:85793/149‑1 (MQ=255)                gCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGAGAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTTCGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGt                                                                                                                                        <  4:3297/148‑1 (MQ=255)                gCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCACTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTa                                                                                                                                       <  4:98641/149‑1 (MQ=255)                 cGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCCGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCTGGGGCCCGGGCTGccaata                                                                                                                                       <  3:195491/148‑6 (MQ=255)                 cGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCCGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCTGGGGCCCGGGCTGccaa                                                                                                                                         >  4:195491/1‑143 (MQ=255)                    gCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGGTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGc                                                                                                                                   >  4:43703/1‑149 (MQ=255)                    gCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCTGGGGCTCGGGCTGc                                                                                                                                            >  2:76740/1‑140 (MQ=255)                    gCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCTGGGGCTCGGGCTGc                                                                                                                                            <  1:76740/140‑1 (MQ=255)                       aTTGCTGGCCAACCTTGCCGCCAACGCCGCGCGACTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTAGGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGc                                                                                                                                <  2:102901/149‑1 (MQ=255)                        ttGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCt                                                                                                                               <  1:65860/149‑1 (MQ=255)                          gCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCtt                                                                                                                              <  3:167908/148‑1 (MQ=255)                          gCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCtt                                                                                                                              >  2:153752/1‑148 (MQ=255)                            tGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGcccc                                                                                                                                                    >  4:1978/1‑124 (MQ=255)                            tGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGcccc                                                                                                                                                    <  3:1978/124‑1 (MQ=255)                                   ccTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAg                                                                                                                    <  7:176905/149‑1 (MQ=255)                                         cgccAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCACTTGAAGAGCAAGGcttct                                                                                                              <  6:192076/149‑1 (MQ=255)                                           ccAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGcttctt                                                                                                             >  2:172387/1‑148 (MQ=255)                                           ccAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCttcttc                                                                                                            >  3:270315/1‑149 (MQ=255)                                           ccAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGCTATTGCCGCTTGAAGAGCAAGGcttctt                                                                                                             <  1:172387/148‑1 (MQ=255)                                              aggccccGCGCCTGCCCTACAGTGAACACCGCATGCGCCGGTTTCTGTCGAAGGCTGCGCCTACGGCTTACGACCAGACCCCAAGGCTATTGGCAGCCCGGGCCCCTGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCgg                                                                                                         <  6:184287/147‑1 (MQ=255)                                                gccgcGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAaga                                                                                                                         <  1:8728/131‑1 (MQ=255)                                                gccgcGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAaga                                                                                                                         >  2:8728/1‑131 (MQ=255)                                                     gcgcCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCt                                                                                                  <  7:32973/149‑1 (MQ=255)                                                       gcCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGc                                                                                                <  5:277941/149‑1 (MQ=255)                                                        ccTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGcc                                                                                               <  1:153752/149‑1 (MQ=255)                                                                aCAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGc                                                                                       <  4:270315/149‑1 (MQ=255)                                                                  aGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGTTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTa                                                                                     <  1:213735/149‑1 (MQ=255)                                                                  aGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGc                                                                                       <  5:86433/147‑1 (MQ=255)                                                                  aGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGc                                                                                       >  6:86433/1‑147 (MQ=255)                                                                    tGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTAt                                                                                    >  5:132525/1‑148 (MQ=255)                                                                        tACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGC                                                                               >  3:91491‑M1/1‑149 (MQ=255)                                                                             cATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGCTCCGCTATTACGCGCAG                                                                           >  3:185551‑M1/1‑148 (MQ=255)                                                                             cATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAG                                                                           >  7:253410‑M1/1‑148 (MQ=255)                                                                                gcgcCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCGAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTAc                                                                        >  7:258576/1‑148 (MQ=255)                                                                                     tGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTCGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTCAAGAGCAAGGCTTCTTCCGGCCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCa                                                                  >  5:30639/1‑149 (MQ=255)                                                                                                aaGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATgg                                                       >  3:76288/1‑149 (MQ=255)                                                                                                    cTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGcacatcggtcgcagacgcgtaatt                                                   >  8:87097‑M1/1‑125 (MQ=255)                                                                                                    cTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAAc                                                   >  2:167260/1‑149 (MQ=255)                                                                                                     ttCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTCCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACg                                                  >  2:141710/1‑149 (MQ=255)                                                                                                       tttttACCTCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCATGGTTGCAGGTATTGCCGTTTGAAGAGCAAGGCTTCTTCCGTGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGt                                                <  7:111290/145‑1 (MQ=255)                                                                                                        gCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTg                                               >  3:241799/1‑149 (MQ=255)                                                                                                           tACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCa                                            >  3:298453/1‑149 (MQ=255)                                                                                                            aCCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCATCGCTGTATGGCAACTGttata                                            >  8:111290/1‑144 (MQ=255)                                                                                                               gCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGc                                        <  1:167260/149‑1 (MQ=255)                                                                                                                    cGAACAGTACGCAAGGCTATTGGGAGCCCGGTCGCCAGGGCTGCAGGTATTGGCGCTTTAAGAGCATGGCTTCTTGCGGGCCTTCTGCCAGAGCCGCTATTACGCGGAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATGGccacc                                    <  6:30639/148‑1 (MQ=255)                                                                                                                     gAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCGATGccaccgt                                  <  2:49612/149‑3 (MQ=255)                                                                                                                         aGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCG                                                                              >  5:259281‑M1/1‑101 (MQ=255)                                                                                                                         aGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCG                                                                              <  6:259281‑M1/101‑1 (MQ=255)                                                                                                                                aaGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTg                       >  4:210179/1‑149 (MQ=255)                                                                                                                                  ggCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGcacatcggtcgcagacgcgtaattggtggcgatgccaccgttgccatccagcgc                     >  5:258070‑M1/1‑95 (MQ=255)                                                                                                                                  ggCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGc                      >  7:11899/1‑148 (MQ=255)                                                                                                                                    atacctgCAGCCCTGGGGCCCGGGCTGCCGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCt                        <  1:126748/138‑1 (MQ=255)                                                                                                                                    atacctgCAGCCCTGGGGCCCGGGCTGCCGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCt                        >  2:126748/7‑144 (MQ=255)                                                                                                                                         ggCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGAt               <  4:241799/148‑1 (MQ=255)                                                                                                                                          gCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGAtgt             >  3:248353/1‑149 (MQ=255)                                                                                                                                          gCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGAtgt             >  8:11911/1‑149 (MQ=255)                                                                                                                                           cAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGAtgtg            <  4:248353/149‑1 (MQ=255)                                                                                                                                            aGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGc           >  4:91491‑M2/1‑148 (MQ=255)                                                                                                                                              cccGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGctg         >  4:76288‑M2/1‑146 (MQ=255)                                                                                                                                               ccGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGctgc        <  8:78322‑M2/149‑5 (MQ=255)                                                                                                                                                 gggCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGcacatcggtcgcagacgcgtaattggtggcgatgccaccgttgccatccagcgctgcctcgtactgcgc      <  6:258070‑M1/149‑70 (MQ=255)                                                                                                                                                 gggCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGcacatcggtcgcagacgcgtaattggtggcgatgccaccgttgccatccagcgctgcctcgtactgcgc      <  7:87097‑M1/149‑70 (MQ=255)                                                                                                                                                     cccAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTAcc  <  1:141710/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GCGACAAAGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACC  >  NC_002947/4921441‑4921737

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