Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,741,884	G→C	6.7%	D90E (GAC→GAG) ‡	PP_5039 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,741,884	0	G	C	6.7%	194.9 / 4.3	60	D90E (GAC→GAG) ‡	PP_5039	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (25/31);  new base C (2/2);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

GATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAAC  >  NC_002947/5741742‑5742028                                                                                                                                                |                                                                                                                                                 gctaccGCTCCACGCACAGGGTCTGCCAGGCCAGCACTGCGGCATGCTCCAGTACCAGCAGCACGCACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAACCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGccc                                                                                                                                            <  5:209702/144‑1 (MQ=255)  aTAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCt                                                                                                                                           <  6:189379/149‑1 (MQ=255)    aaCGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGCGCGCCTTc                                                                                                                                          >  4:165726/1‑148 (MQ=255)                acaGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGCGCGCCTTCGAGGGCGAccggc                                                                                                                             <  3:165726/149‑1 (MQ=255)                    gTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCgtgt                                                                                                                         >  3:263488/1‑149 (MQ=255)                     ttCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTc                                                                                                                        >  8:258830/1‑149 (MQ=255)                        tGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTcgcg                                                                                                                     >  5:81974/1‑149 (MQ=255)                         gCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCgg                                                                                                                    >  7:129081/1‑149 (MQ=255)                         gCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCgg                                                                                                                    >  2:127127/1‑149 (MQ=255)                          ccAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCgg                                                                                                                    >  8:140260/1‑148 (MQ=255)                          ccAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCgg                                                                                                                    <  7:142835/148‑1 (MQ=255)                          ccAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGt                                                                                                                   <  6:81974/149‑1 (MQ=255)                           cAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCgg                                                                                                                    >  6:330628/1‑147 (MQ=255)                            aGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAg                                                                                                                 <  1:127127/149‑1 (MQ=255)                                     cTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGt                                                                                                        <  6:231662/149‑1 (MQ=255)                                      ttCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTg                                                                                                       >  7:163812/1‑149 (MQ=255)                                           cATGCTCCAGTACCAGCAGCACGAACACGCGCTGGACCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCAt                                                                                                            <  1:253474/139‑1 (MQ=255)                                           cATGCTCCAGTACCAGCAGCACGAACACGCGCTGGACCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCAt                                                                                                            >  2:253474/1‑139 (MQ=255)                                             tGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTg                                                                                                <  6:228630/149‑1 (MQ=255)                                                  cAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGa                                                                                           <  2:320213/149‑1 (MQ=255)                                                    gTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAg                                                                                          >  3:203164/1‑148 (MQ=255)                                                       ccAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCaaa                                                                                      <  7:258830/149‑1 (MQ=255)                                                        cagcagCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCaaaa                                                                                     >  5:65916/1‑149 (MQ=255)                                                         agcagcACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAAt                                                                                    <  1:112964/149‑1 (MQ=255)                                                         agcagcACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGACCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGt                                                                                                                   >  5:324440/1‑118 (MQ=255)                                                         agcagcACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGACCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGt                                                                                                                   <  6:324440/118‑1 (MQ=255)                                                          gcagcaCGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGCACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTcg                                                                                                                       <  1:67075/114‑1 (MQ=255)                                                          gcagcaCGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTcg                                                                                                                       >  2:67075/1‑113 (MQ=255)                                                               aCGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATcagcag                                                                              <  5:330628/149‑1 (MQ=255)                                                                cGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAgg                                                                             <  3:286974/149‑1 (MQ=255)                                                                cGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAgg                                                                             >  1:7930/1‑149 (MQ=255)                                                                 gAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGc                                                                            >  5:142204/1‑149 (MQ=255)                                                                 gAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGc                                                                            >  2:315713/1‑149 (MQ=255)                                                                  aaCACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGcc                                                                           >  4:232125/1‑149 (MQ=255)                                                                          cTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCcagccagc                                                                   >  1:56444/1‑149 (MQ=255)                                                                                 ggCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCgg                                                            <  6:142204/149‑1 (MQ=255)                                                                                       agcTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTtg                                                      >  3:153401/1‑149 (MQ=255)                                                                                       agcTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTtg                                                      <  1:315713/149‑1 (MQ=255)                                                                                          tGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTtgctg                                                   >  3:256195/1‑149 (MQ=255)                                                                                          tGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTtgctg                                                   >  5:227909/1‑149 (MQ=255)                                                                                            gTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTtgctgc                                                  <  7:147497/148‑1 (MQ=255)                                                                                            gTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTtgctgc                                                  <  7:140260/148‑1 (MQ=255)                                                                                            gTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGcc                                                 <  4:91173/149‑1 (MQ=255)                                                                                                   cAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTc                                          <  6:65916/149‑1 (MQ=255)                                                                                                     gCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCATGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGa                                        <  1:233720/149‑1 (MQ=255)                                                                                                       ccTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGc                                      <  4:263488/149‑1 (MQ=255)                                                                                                           cTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCAc                                  >  5:81154/1‑149 (MQ=255)                                                                                                              cAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCg                               >  8:180943/1‑149 (MQ=255)                                                                                                                       cATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGc                        <  8:163812/147‑1 (MQ=255)                                                                                                                            cG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCg                 <  2:7930/149‑1 (MQ=255)                                                                                                                             g‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCTGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGa                <  4:203164/149‑1 (MQ=255)                                                                                                                             g‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGa                <  6:81154/149‑1 (MQ=255)                                                                                                                             g‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGa                <  8:129081/149‑1 (MQ=255)                                                                                                                                 caccTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCa             <  6:227909/149‑1 (MQ=255)                                                                                                                                 caccTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCa             <  4:256195/149‑1 (MQ=255)                                                                                                                                    cTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGt          >  8:79231/1‑149 (MQ=255)                                                                                                                                      cgacACGCCGCGCGCCTTCGAGGGCGACCGGCGTGCCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGa                                        >  4:310367/4‑117 (MQ=37)                                                                                                                                      cgacACGCCGCGCGCCTTCGAGGGCGACCGGCGTGCCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGa                                        <  3:310367/114‑1 (MQ=37)                                                                                                                                          aCGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCAAGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCaa    >  3:245319/1‑149 (MQ=255)                                                                                                                                            gCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAAc  <  7:180943/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                 GATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCG‑ACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAAC  >  NC_002947/5741742‑5742028

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