Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,338,493	A→T	9.2%	intergenic (+113/+29)	PP_3812 → / ← PP_5632	lipase/allo‑isoleucine metabolism protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,338,493	0	A	T	9.2%	145.0 / 6.0	62	intergenic (+113/+29)	PP_3812/PP_5632	lipase/allo‑isoleucine metabolism protein
Reads supporting (aligned to +/- strand):  ref base A (23/33);  new base T (3/3);  total (26/36)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.24e-01

CAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCA‑‑CAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATC  >  NC_002947/4338353‑4338638                                                                                                                                             |                                                                                                                                                    cAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCcgcg                                                                                                                                             <  7:149765/149‑1 (MQ=255) cAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCcgc                                                                                                                                              >  7:84955/1‑148 (MQ=255)  aaCGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCcgcg                                                                                                                                             >  1:44785/1‑148 (MQ=255)     gAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGAAAg                                                                                                                                         <  2:44785/149‑1 (MQ=255)     gAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAg                                                                                                                                         >  6:79252/1‑149 (MQ=255)            aaGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGc                                                                                                                                  <  3:88313/149‑1 (MQ=255)              gTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCa                                                                                                                                 >  5:293371/1‑148 (MQ=255)                  aTCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAAACCGCGAAAGGGCCGGCAAATcc                                                                                                                            >  3:207613/1‑149 (MQ=255)                    cTGCAGTGAGGCGATGAGCAATAGACCCGCGTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTg                                                                                                                          >  1:29872/1‑149 (MQ=255)                    cTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTGCCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTg                                                                                                                          >  1:245208/1‑149 (MQ=255)                         gTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAg                                                                                                                     >  5:128538/1‑149 (MQ=255)                           gAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAgg                                                                                                                                        <  6:278364/128‑1 (MQ=255)                           gAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAgg                                                                                                                                        >  5:278364/1‑128 (MQ=255)                            aGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGTAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCTAAAGGGCCGGCAAATCCTGCTCAggg                                                                                                                   <  4:158399/148‑1 (MQ=255)                                   gAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAAc                                                                                                            <  3:148948/148‑1 (MQ=255)                                       aaTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGACTCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTc                                                                                                       <  4:203698/149‑1 (MQ=255)                                       aaTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTc                                                                                                       <  1:6212/149‑1 (MQ=255)                                       aaTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTc                                                                                                       <  6:128538/149‑1 (MQ=255)                                        aTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATGCTGCTCAGGGCTCCAACCCTTcc                                                                                                      >  7:227474/1‑149 (MQ=255)                                        aTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTcc                                                                                                      <  8:140479/149‑1 (MQ=255)                                        aTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGGAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTcc                                                                                                      >  7:7406/1‑149 (MQ=255)                                         tAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCca                                                                                                     <  6:293371/149‑1 (MQ=255)                                           gACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCca‑‑c                                                                                                  <  5:68281/148‑1 (MQ=255)                                            aCCGGCCTTGGCGGGTCTTTTGCTTCATTTTGCGCAAACAGCCCTCAGGCTTCAATTATTCTGCCAGCGTTGGGGCAGGTCCGGCCTTTTGGTGGGTAAACCCGCGAATGGGCCGGCAAATCCTGCTCAGGGCTCAAAc                                                                                                            <  8:46504/139‑1 (MQ=255)                                            aCCCGCCTGGGCGGGTGTTTTGCTTCTATTTGCGGAAACCGCCCGCTGGCTGCAATGATTCTGCCAGGGGGGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCGCAGGGCTCAAACCCTTCCA‑‑Caa                                                                                                <  2:116815/149‑1 (MQ=255)                                            aCCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCca                                                                                                     <  5:84653/146‑1 (MQ=255)                                            aCCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCca                                                                                                     >  6:84653/1‑146 (MQ=255)                                            aCCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCA‑‑Caa                                                                                                <  4:20092/149‑1 (MQ=255)                                              ccGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGaaa                                                                                                                                          >  3:118400/1‑107 (MQ=255)                                              ccGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGaaa                                                                                                                                          <  4:118400/107‑1 (MQ=255)                                              ccGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCA‑‑CAATg                                                                                              <  1:1571/149‑1 (MQ=255)                                                gCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTAGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCA‑‑CAATGa                                                                                             >  3:220313/1‑148 (MQ=255)                                                     ggCGGGTCTTTTGCTTCTATTTGCGGAAACATCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAACCCTGCTCAGGGCTCAAACCCTTCCA‑‑CAATGATCAc                                                                                         <  5:79252/147‑1 (MQ=255)                                                     ggCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCA‑‑CAATGATCACCt                                                                                       <  7:22757/149‑1 (MQ=255)                                                              tttGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCTAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCA‑‑CAATGATCACCTCAGCCTTc                                                                               <  6:198663/148‑1 (MQ=255)                                                                  cTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTaaa                                                                                                                                                   <  1:311939/78‑1 (MQ=255)                                                                  cTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTaaa                                                                                                                                                   >  2:311939/1‑78 (MQ=255)                                                                   ttCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCCGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCA‑‑CAATGATCACCTCAGCCTTCGCCAcc                                                                         >  1:242236/1‑149 (MQ=255)                                                                      tATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCA‑‑CAATGATCACCTCAGCCTTCGCCACCCCt                                                                      >  1:20722/1‑149 (MQ=255)                                                                             ggAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCA‑‑CAATGATCACCTCAGCCTTCGCCACCCCTTGGCGAt                                                               <  4:325876/149‑1 (MQ=255)                                                                             ggAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCCGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCATCCA‑‑CAATGATCACCTCAGCCTTCGCCACCCCATGGCGAt                                                               >  1:1427/1‑149 (MQ=255)                                                                             ggAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCATTAGCGGGTAAAACCGCGAAAGGGGCGGCAAAAACTGCTAAGGGCTCAAAACCTTACA‑‑CAATGAAAACCACAGCCATCGCCACCCCATGGCGAt                                                               >  4:153147/1‑149 (MQ=255)                                                                                 aCAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCA‑‑CAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCt                                                           <  2:245208/149‑1 (MQ=255)                                                                                 aCAGCCCGCAGGCTTCAATGAATCTGCCAGCGGGGCGGCAGGGCCGCCCCTTTCGCGGTAAACCCCGCTAATGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCA‑‑CAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCt                                                           <  2:327017/149‑1 (MQ=255)                                                                                   aGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTAGCGGGTAAAACCGCGAAACGGGCGGCAAAACCTGCTGAGGGCTCAAACCCATACA‑‑CAATGATAACCACAGCCTTAGCCACCCCGTGGCGATGGcag                                                          >  2:94937/1‑146 (MQ=255)                                                                                     cccGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCA‑‑CAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCAt                                                       <  4:220313/149‑1 (MQ=255)                                                                                     cccGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCA‑‑CAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCAt                                                       <  4:207613/149‑1 (MQ=255)                                                                                           ggCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCA‑‑CAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCctct                                                 <  2:20722/149‑1 (MQ=255)                                                                                              ttCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGGCGGCCAATCCTGCTCAGGGCTCAAACCCGTGCA‑‑CAATGATCACCTCAGCCTTCGCCACCCCGTGGCGCTGGCTGGATGCCTCTTgg                                              >  1:281548/1‑149 (MQ=255)                                                                                                     aagcTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCA‑‑CAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAgcg                                       <  7:135003/146‑1 (MQ=255)                                                                                                      aTCTGCCAGCGGTGGGGGATGTGCGGGCCTTTTGCGGGTTAAACCGCGAAATGGGCGGGAAATACTTCTCACGGCGCAAAACCTTTCACTCAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAg                                         <  1:94937/148‑1 (MQ=255)                                                                                                           ccAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCA‑‑CAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAg                                  <  8:84955/148‑1 (MQ=255)                                                                                                             aggggTGCGGCTGGTCGGGCCTTTTCGGGGGTAAACCCGCGAATGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCA‑‑CAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCgg                               <  4:26056/146‑1 (MQ=255)                                                                                                                ggTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCA‑‑CAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCCGAGCGGTa                             >  2:32310/1‑148 (MQ=255)                                                                                                                         ggTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCGCAAACCCTTCCC‑‑CAATGATCACCACAGCCTTACCCCCCCCGTGGCGATGGCTGCATGCCGCTTGGTACAGCGCGGAGCGGTAACCGGCGCCg                   >  6:256075/1‑149 (MQ=255)                                                                                                                             cGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCA‑‑CAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCt               >  8:305778/1‑149 (MQ=255)                                                                                                                              ggCCCTTTCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCA‑‑CAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTg              >  6:120075/1‑149 (MQ=255)                                                                                                                               gCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCA‑‑CAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTgc             <  2:242236/149‑1 (MQ=255)                                                                                                                                     tCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCA‑‑CAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTa        <  7:305778/148‑1 (MQ=255)                                                                                                                                      cgcgGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCA‑‑CAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCCGAGCGGTAACAGGCGACGGCCTGCGCGTAc       <  1:32310/148‑1 (MQ=255)                                                                                                                                        cgGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCA‑‑CAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGaa    <  3:153147/149‑1 (MQ=255)                                                                                                                                         gggTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCA‑‑CAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAAt   <  2:29872/149‑1 (MQ=255)                                                                                                                                         gggTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCA‑‑CAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAAt   >  3:273931/1‑149 (MQ=255)                                                                                                                                           gTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCA‑‑CAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATc  >  2:69700/1‑148 (MQ=255)                                                                                                                                             |                                                                                                                                                    CAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCA‑‑CAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATC  >  NC_002947/4338353‑4338638

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