Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,625,416	G→T	6.1%	G24G (GGC→GGA)	PP_3195 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,625,416	0	G	T	6.1%	210.8 / 3.6	66	G24G (GGC→GGA)	PP_3195	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (28/34);  new base T (2/2);  total (30/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.88e-01

CCTGCAGCGGTGCGCACATTTCCTCGATTTCACCCATGCGGGCACCCATTTCGATGGCCAGATCCACCGAGAGGTTCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCG‑AAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGAT  >  NC_002947/3625274‑3625554                                                                                                                                               |                                                                                                                                             ccTGCAGCGGTGCGCACATTTCCTCGATTTCACCCATGCGGGCACCCATTTCGATGGCCAGATCCACCGAGAGGTTCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCa                                                                                                                                          <  4:408021/147‑1 (MQ=255)   tGCAGCGGTGCGCACATTTCCTCGATTTCACCCATGCGGGCACCCATTTCGATGGCCAGATCCACCGAGAGGTTCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACtt                                                                                                                                       >  7:115008/1‑148 (MQ=255)       gCGGTGCGCACATTTCCTCGATTTCACCCATGCGGGCACCCATTTCGATGGCCAGATCCACCGAGAGGTTCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGc                                                                                                                                   >  3:25171/1‑148 (MQ=255)                acaTTTCCTCGATTTCACCCATGCGGGCACCCATTTCGATGGCCAGATCCACCGAGAGGTTCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGAcaggcagg                                                                                                                         >  7:142440/1‑149 (MQ=255)                  aTTTCCTCGATTTCACCCATGCGGGCACCCATTTCGATGGCCAGATCCACCGAGAGGTTCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGAcaggca                                                                                                                           >  1:403107/1‑145 (MQ=255)                  aTTTCCTCGATTTCACCCATGCGGGCACCCATTTCGATGGCCAGATCCACCGAGAGGTTCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGAcaggca                                                                                                                           <  2:403107/145‑1 (MQ=255)                      ccTCGATTTCACCCATGCGGGCACCCATTTCGATGGCCAGATCCACCGAGAGGTTCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGgc                                                                                                                   >  2:44512/1‑149 (MQ=255)                       cTCGATTTCACCCATGCGGGCACCCATTTTGATGGCCAGATCCACCGAGAGGTTCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGgc                                                                                                                   <  3:245890/148‑1 (MQ=255)                       cTCGATTTCACCCATGCGGGCACCCATTTCGATGGCCAGATCCACCGAGAGGTTCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGAcaggc                                                                                                                            >  7:103692/1‑139 (MQ=255)                       cTCGATTTCACCCATGCGGGCACCCATTTCGATGGCCAGATCCACCGAGAGGTTCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGAcaggc                                                                                                                            <  8:103692/139‑1 (MQ=255)                       cTCGATTTCACCCATGCGGGCACCCATTTCGATGGCCAGATCCACCGAGAGGTTCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGgcg                                                                                                                  >  5:329573/1‑149 (MQ=255)                          gttttCACCCATGCGGGCACCCATTTCGATGGCCAGATCCACCGAGAGGTTCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGcgcg                                                                                                                <  7:138514/146‑1 (MQ=255)                          gATTTCACCCATGCGGGCACCCATTTCGATGGCCAGATCCACCGAGAGGTTCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCgg                                                                                                               <  4:25171/149‑1 (MQ=255)                            tttCACCCATGCGGGCACCCATTTCGATGGCCAGATCCACCTAGAGGTTCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCgg                                                                                                               <  8:142440/147‑1 (MQ=255)                                aCCCATGCGGGCACCCATTTCGATGGCCAGATCCACCGAGAGGTTCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGaa                                                                                                         <  1:44512/149‑1 (MQ=255)                                 cccATGCGGGCACCCATTTCGATGGCCAGATCCACCGAGAGGTTCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAg                                                                                                        <  8:88799/149‑1 (MQ=255)                                 cccATGCGGGCACCCATTTCGATGGCCAGATCCACCGAGAGGTTCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAg                                                                                                        <  6:329573/149‑1 (MQ=255)                                     tGCGGGCACCCATTTCGATGGCCAGATCCACCGAGAGGTTCCAGACGTAAGTGGTGGGGAAGTAACGGAACATTCTGGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGc                                                                                                      >  2:119042/1‑147 (MQ=255)                                       cGGGCACCCATTTCGATGGCCAGATCCACCGAGAGGTTCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAg                                                                                                   <  8:115008/148‑1 (MQ=255)                                           cACCCATTTCGATGGCCAGATCCACCGAGAGGTTCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACg                                                                                              <  6:368433/149‑1 (MQ=255)                                            aCCCATTTCGATGGCCAGATCCACCGAGAGGTTCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCTACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACgc                                                                                             <  7:134743/149‑1 (MQ=255)                                            aCCCATTTCCAAGTCCAGATCCACAGAGACGTACCAGACATAAATGGAGGGGAAGTAACCGAACATTCTTGTAGTCATTCTGCCGTCTTGGAGTCCGGCGCCCACTTCAGCGACAGGAAGCCAGGGCGCGCATGGAAGGCAAGGTACgc                                                                                             >  6:247335/1‑149 (MQ=255)                                              ccATTTCGATGGCCAGATCCACCGAGAGGTTCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGAcaggc                                                                                                                            <  2:35962/116‑1 (MQ=255)                                              ccATTTCGATGGCCAGATCCACCGAGAGGTTCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGAcaggc                                                                                                                            >  1:35962/1‑116 (MQ=255)                                              ccATTTCGATGGCCAGATCCACCGAGAGGTTCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACgc                                                                                             >  7:210771/1‑147 (MQ=255)                                              ccATTTCGATGGCCAGATCCACCGAGAGGTTCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACgc                                                                                             <  8:210771/147‑1 (MQ=255)                                              ccATTTCGATGGCCAGATCCACCGAGAGGTTCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCCACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACgcag                                                                                           >  6:63585/1‑149 (MQ=255)                                                aTTTCGATGGCCAGATCCACCGAGAGGTTCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCgg                                                                                                               <  1:114327/127‑1 (MQ=255)                                                aTTTCGATGGCCAGATCCACCGAGAGGTTCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCgg                                                                                                               >  2:114327/1‑127 (MQ=255)                                                          ccAGATCCACCGAGAGGTTCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCTCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGt                                                                               >  6:55300/1‑149 (MQ=255)                                                          ccAGATCCACCGAGAGGTTCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCTCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGt                                                                               >  6:61442/1‑149 (MQ=255)                                                            aGATCCACCGAGAGGTTCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTgg                                                                             <  3:151465/149‑1 (MQ=255)                                                                      agagGTTCCAGACGTAATTGGTTGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGt                                                                   <  8:49918/149‑1 (MQ=255)                                                                           ttCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCTCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCa                                                              <  5:55300/149‑1 (MQ=255)                                                                           ttCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCTCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCa                                                              <  5:61442/149‑1 (MQ=255)                                                                           ttCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGAt                                                                >  7:344738/1‑147 (MQ=255)                                                                                gACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGa                                                          >  7:339431/1‑148 (MQ=255)                                                                                catgTAATTGGGGGGGAAGTAACGGAACTTTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCAGTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAgg                                                            <  8:339431/143‑1 (MQ=255)                                                                                   gTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACgca                                                                                            <  4:344638/111‑1 (MQ=255)                                                                                   gTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACgca                                                                                            >  3:344638/1‑111 (MQ=255)                                                                                   gTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCCACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCg                                                      <  5:63585/149‑1 (MQ=255)                                                                                      aTTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTg                                                                              >  1:70008/1‑122 (MQ=255)                                                                                      aTTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTg                                                                              <  2:70008/122‑1 (MQ=255)                                                                                      aTTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAg                                                    <  3:138007/148‑1 (MQ=255)                                                                                                   tAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCa                                                              <  6:282615/125‑1 (MQ=255)                                                                                                   tAACGGAACATTATTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCa                                                              >  5:282615/1‑125 (MQ=255)                                                                                                       ggAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATcg                                              <  3:152223/137‑1 (MQ=255)                                                                                                       ggAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCAACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATcg                                              >  4:152223/1‑137 (MQ=255)                                                                                                        gAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCG‑AAAcg                                      <  4:195812/143‑1 (MQ=255)                                                                                                        gAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCG‑AAAcg                                      >  3:195812/1‑143 (MQ=255)                                                                                                          aCATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGAGCAGGGATGCGAGGTATCGCG‑AAACGCGATGccc                              >  2:6497/1‑149 (MQ=255)                                                                                                             ttCTTTTTGTCCTTCTGCGGTCGTGGATTCCGGCGCCCAATTTAGCGACAGGCAGGGATGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCGAATGGTGATCAGGGATGCGATGTATCGCG‑AAACGCGATGCCCCAg                           <  5:247335/149‑1 (MQ=255)                                                                                                                ttgttgTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGa                                                     <  1:189183/121‑1 (MQ=255)                                                                                                                ttgttgTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGa                                                     >  2:189183/1‑121 (MQ=255)                                                                                                                ttgttgTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCG‑AAACGCGATGCCCCAG‑CAc                       >  2:14221/1‑149 (MQ=255)                                                                                                                         ttCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCG‑AAACGCGATGCCCCAG‑CACCGCTGCGgt              >  8:306786/1‑149 (MQ=255)                                                                                                                          tCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGCATGAGAGCTATCGCG‑AAACGCGATGCCCCAGACACCGCAGCGgt              >  2:257580/1‑149 (MQ=255)                                                                                                                           cTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCG‑AAACGCGATGCCCCAG‑CACCGCTGCGgt              <  1:119042/147‑1 (MQ=255)                                                                                                                           cTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCG‑AAACGCGATGCCCCAG‑CACCGCTGCGgt              <  8:344738/147‑1 (MQ=255)                                                                                                                                gTCGTGGTGTCCGGCGCCCACTTGAGCGACAGGGAGGCAGGGCGCGGTTTGTAGGCTAGGTAGGCAGGCAGGGCAGGGTGGGGCCAATGGTGATCAGGGATGCGAGTTATCGCGAAAACGCGATGCCCCAGACACCGCAGCGgtgctgg         <  1:257580/149‑5 (MQ=255)                                                                                                                                gTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGCGATGCGAGGTATCGCG‑AAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTggg       >  3:301893/1‑149 (MQ=255)                                                                                                                                 tCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCg                                                      >  5:367468/1‑103 (MQ=255)                                                                                                                                 tCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCg                                                      <  6:367468/103‑1 (MQ=255)                                                                                                                                     ggAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCG‑AAACGCGATGcccc                             <  6:67355/123‑1 (MQ=255)                                                                                                                                     ggAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCG‑AAACGCGATGcccc                             >  5:67355/1‑123 (MQ=255)                                                                                                                                      gAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCG‑AAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGAt  <  1:14221/148‑1 (MQ=255)                                                                                                                                               |                                                                                                                                             CCTGCAGCGGTGCGCACATTTCCTCGATTTCACCCATGCGGGCACCCATTTCGATGGCCAGATCCACCGAGAGGTTCCAGACGTAATTGGTGGGGAAGTAACGGAACATTCTTGTTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCG‑AAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGAT  >  NC_002947/3625274‑3625554

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