Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,341,861	G→C	6.2%	V41V (GTG→GTC)	panC →	pantothenate synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,341,861	0	G	C	6.2%	206.3 / 4.8	65	V41V (GTG→GTC)	panC	pantothenate synthetase
Reads supporting (aligned to +/- strand):  ref base G (25/36);  new base C (2/2);  total (27/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 = 7.83e-01

CCAGGCAGCGAACACGGGTTCAGTGCATGAATACAGTCAAGACCGTCCGCGAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCAGTATCTTCGTCAACCCGCTGCAGTTCGGCGCCAACGAAGACCTCGACAAATACCCACGCACCCTGGCAGCCGACCAGGAGCGCCTGGTGCAGGCCGGCTGCAACCTGCTGTT  >  NC_002947/5341713‑5342010                                                                                                                                                     |                                                                                                                                                      ccAGGCAGCGAACACGGGTTCAGTGCATGAATACAGTCAAGACCGTCCGCGAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGa                                                                                                                                                      <  1:680642/150‑1 (MQ=255) ccAGGCAGCGAACACGGGTTCAGTGCATGAATACAGTCAAGACCGTCCGCGAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGa                                                                                                                                                      >  2:704335/1‑150 (MQ=255)    ggCAGCGAACACGGGTTCAGTGCATGAATACAGTCAAGACCGTCCGCGAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCa                                                                                                                                                   <  2:949111/150‑1 (MQ=255)       aGCGAACACGGGTTCAGTGCATGAATACAGTCAAGACCGTCCGCGAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCaatt                                                                                                                                                <  1:347937/150‑3 (MQ=255)       aGCGAACACGGGTTCAGTGCATGAATACAGTCAAGACCGTCCGCGAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAgg                                                                                                                                                <  2:165850/150‑1 (MQ=255)          gAACACGGGTTCAGTGCATGAATACAGTCAAGACCGTCCGCGAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTg                                                                                                                                             >  2:1150506/1‑150 (MQ=255)             cacGGGTTCAGTGCATGAATACAGTCAAGACCGTCCGCGAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGccgcat                                                                                                                                          <  1:827946/150‑6 (MQ=255)                gggTTCAGTGCATGAATACAGTCAAGACCGTCCGCGAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAgc                                                                                                                                       >  2:603229/1‑150 (MQ=255)                    tCAGTGCATGAATACAGTCAAGACCGTCCGCGAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAgcgcgc                                                                                                                                   <  2:547995/150‑1 (MQ=255)                     cAGTGCATGAATACAGTCAAGACCGTCCGCGAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGcc                                                                                                                                  <  1:1101250/150‑1 (MQ=255)                      aGTGCATGAATACAGTCAAGACCGTCCGCGAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGcc                                                                                                                                           >  1:1011143/1‑140 (MQ=255)                      aGTGCATGAATACAGTCAAGACCGTCCGCGAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGcc                                                                                                                                           <  2:1011143/140‑1 (MQ=255)                      aGTGCATGAATACAGTCAAGACCGTCCGCGAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGACATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCg                                                                                                                                 >  1:869635/1‑150 (MQ=255)                         gCATGAATACAGTCAAGACCGTCCGCGAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCTAGCGCGCCGACt                                                                                                                              >  2:82170/1‑150 (MQ=255)                           aTGAATACAGTCAAGACCGTCCGCGAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTc                                                                                                                            <  1:103575/150‑1 (MQ=255)                             gAATACAGTCAAGACCGTCCGCGAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGt                                                                                                                          >  2:685484/1‑150 (MQ=255)                             gAATACAGTCAAGACCGTCCGCGAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGt                                                                                                                          <  2:402512/150‑1 (MQ=255)                              aaTACAGTCAAGACCGTCCGCGAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTg                                                                                                                         <  1:50251/150‑1 (MQ=255)                              aaTACAGTCAAGACCGTCCGCGAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTg                                                                                                                         <  1:838506/150‑1 (MQ=255)                                   aGTCAAGACCGTCCGCGAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAgg                                                                                                                                                <  2:818676/122‑1 (MQ=255)                                   aGTCAAGACCGTCCGCGAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAgg                                                                                                                                                >  1:818676/1‑122 (MQ=255)                                   aGTCAAGACCGTCCGCGAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGc                                                                                                                    >  2:468447/1‑150 (MQ=255)                                     tCAAGACCGTCCGCGAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGc                                                                                                                    <  2:249228/148‑1 (MQ=255)                                     tCAAGACCGTCCGCGAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGc                                                                                                                    >  1:249228/1‑148 (MQ=255)                                     tCAAGACCGTCCGCGAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCa                                                                                                                  <  2:1027239/150‑1 (MQ=255)                                       aaGACCGTCCGCGAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTg                                                                                                                         >  1:896565/1‑141 (MQ=255)                                       aaGACCGTCCGCGAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTg                                                                                                                         <  2:896565/141‑1 (MQ=255)                                          gatctCCGCGAACTGCGCGCGGCAGTGGCCAGCGGCAGTGGTGTTGCTTAAAATGTCGGTTGCATACCTACTATGGGCAACCTGCACAGCTGCCATGCGGCCTTGGTGCCCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGTCAGTATc                                                                                                             <  2:265203/146‑1 (MQ=255)                                          aCCGTCCGCGAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCAGTATc                                                                                                             <  2:717445/150‑1 (MQ=255)                                             gTCCGCGAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCAGTATCTTc                                                                                                          <  1:603229/150‑1 (MQ=255)                                             gTCCGCGAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCAGTATCTTc                                                                                                          <  2:1089903/150‑1 (MQ=255)                                                cgcgAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTAGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCAGTATCtagatc                                                                                                       >  1:265203/1‑145 (MQ=255)                                                 gcgAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCAGTATCTTCGTCa                                                                                                      >  2:294557/1‑150 (MQ=255)                                                         cgcgcgGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCAGTATCTTCGTCAACCCGCTg                                                                                              >  2:966635/1‑150 (MQ=255)                                                            gcgGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTCACCAAGGc                                                                                                                                               <  2:567175/98‑1 (MQ=255)                                                            gcgGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTCACCAAGGc                                                                                                                                               >  1:567175/1‑98 (MQ=255)                                                                  gTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCAGTATCTTCGTCAACCCACTGCAGTTCGgc                                                                                     <  2:359257/150‑1 (MQ=255)                                                                     gcccgcgcccgcGGTGCGGGCAACCGCATCGGTTTCGTACCGACTATGGGCCACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCAGTATCTTCGTCAACCAGCTGCAGTTCGGCGcc                                                                                  >  2:1122474/1‑150 (MQ=255)                                                                     gcccgcgcccgcGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCAATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCAGTATCTTCGTCAACCCGCTGCAGTTCGGCGcc                                                                                  <  2:475305/150‑1 (MQ=255)                                                                       ccgcgcccgcgGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGt                                                                                                                          >  1:882890/1‑108 (MQ=255)                                                                       ccgcgcccgcgGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGt                                                                                                                          <  2:882890/108‑1 (MQ=255)                                                                       ccgcgcccgcgGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCAGTATCTTc                                                                                                          <  2:654872/124‑1 (MQ=255)                                                                       ccgcgcccgcgGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCAGTATCTTc                                                                                                          >  1:654872/1‑124 (MQ=255)                                                                         gcgcccgcgGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCAGTATCTTCGTCAACCCGCTGCAGTTCGGCGCCAACg                                                                              <  1:1135581/150‑1 (MQ=255)                                                                                 ggTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCAGTATCTTCGTCAACCCGCTGCAGTTCGGCGCCAACGAAGACCTc                                                                      <  2:855569/150‑1 (MQ=255)                                                                                  gTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCTATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCAGTATATTCGTCAACCCGCTGCAGTTCGGCGCCAACGAAGACCTCg                                                                     <  1:1150506/150‑1 (MQ=255)                                                                                     aGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCAGTATCTTCGTCAACCCGCTGCAGTTCGGCGCCAACGAAGACCTCGACa                                                                  >  2:787726/1‑150 (MQ=255)                                                                                        gCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCAGTATCTTCGTCAACCCGCTGCAGTTCGGCGCCAACGAAGACCTCGACAAAt                                                               >  2:584319/1‑150 (MQ=255)                                                                                           aaCGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCAGTATCTTCGTCAACCCGCTGCAGTTCGGCGCCAACGAAGACCTCGACAAATAcc                                                            <  1:787726/150‑1 (MQ=255)                                                                                              gCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCg                                                                                                                           >  2:451668/1‑84 (MQ=255)                                                                                              gCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCg                                                                                                                           <  1:451668/84‑1 (MQ=255)                                                                                                    gTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCAGTATCTTCGTCAACCCGCTGCAGTTCGGCGCCAACGAAGACCTCGACAAATACCCACGCAccc                                                   >  1:256854/1‑150 (MQ=255)                                                                                                          tACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCAGTATCTTCGTCAACCCGCTGCAGTTCGgcgc                                                                                   <  1:1108275/112‑1 (MQ=255)                                                                                                          tACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCAGTATCTTCGTCAACCCGCTGCAGTTCGgcgc                                                                                   >  2:1108275/1‑112 (MQ=255)                                                                                                                         aCCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCAGTATCTTCGTCAACCCGCTGCAGTTCGGCGCCAACGAAGACCTCGACAAATACCCACGCACCCTGGCAGCCGACCAGGAGCGcc                              >  1:241484/1‑150 (MQ=255)                                                                                                                         aCCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCAGTATCTTCGTCAACCCGCTGCAGTTCGGCGCCAACGAAGACCCCGACAAATACCCACGCACCCTGGCAGCCGACCAGGAGCGcc                              <  2:686391/150‑1 (MQ=255)                                                                                                                           cTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCAGTATCTTCGTCAACCCGCTGCAGTTCGGCGCCAACGAAGACCTCGACAAATACCCACGCACCCTGGCAGCCGACCAGGAGCGCCTg                            <  1:468447/150‑1 (MQ=255)                                                                                                                           cTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCAGTATCTTCGTCAACCCGCTGCAGTTCGGCGCCAACGAAGACCTCGACAAATACCCACGCACCCTGGCAGCCGACCAGGAGCGCCTg                            >  1:872809/1‑150 (MQ=255)                                                                                                                            tGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCAGTATCTTCGTCAACCCGCTGCAGTTCGGCGCCAACGAAGACCTCGACAAATACCCACGCACCCTGGCAGCCGACCAGGAGCGCCTgg                           <  1:1122474/150‑1 (MQ=255)                                                                                                                                       cATGCGGCCTTGGTGACCAAGGCTGCCTAGCGCGCCGACTTCGTGGTCGCCAGTATCTTCGTCAACCCGCTGCAGTTCGGCGCCAACGAAGACCTCGACAAATACCCACGCACCCTGGCAGCCGACCAGGAGCGCCTGGTGCAGGCCGGc                <  1:82170/150‑1 (MQ=255)                                                                                                                                           cagCCTTGGTCACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCAGTATCTTCGTCAACCCGCTGCAGTTCGGCGCCAACGAAGACCTCGACAAATACCCACGCACCCTGGCAGCCGACCAGGAGCGCCTGGTGCAg                      >  2:1039637/3‑140 (MQ=255)                                                                                                                                           cagCCTTGGTCACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCAGTATCTTCGTCAACCCGCTGCAGTTCGGCGCCAACGAAGACCTCGACAAATACCCACGCACCCTGGCAGCCGACCAGGAGCGCCTGGTGCAg                      <  1:1039637/138‑1 (MQ=255)                                                                                                                                                   gTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCAGTATCTTCGTCAACCCGCTGCAGTTCGgc                                                                                     >  2:665822/1‑69 (MQ=255)                                                                                                                                                   gTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCAGTATCTTCGTCAACCCGCTGCAGTTCGgc                                                                                     <  1:665822/69‑1 (MQ=255)                                                                                                                                                   gTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCAGTATCTTCGTCAACCCGCTGCAGTTCGGCGCCAACGAAGACCTCGACAAATACCCACGCACCCTGGCAGCCGACCAGGAGCGCCTGGTGCAGGCCGGCTGCAACctgctg    <  1:294557/150‑1 (MQ=255)                                                                                                                                                     gACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCAGTATCTTCGTCAACCCGCTGCAGTTCGGCGCCAACGAAGACCTCGACAAATACCCACGCACCCTGGCAGCCGACCAGGAGCGCCTGGTGCAGGCCGGCTGCAACCTGCTGtt  <  2:241484/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      CCAGGCAGCGAACACGGGTTCAGTGCATGAATACAGTCAAGACCGTCCGCGAACTGCGCGCGGCAGTGGCCCGCGCCCGCGGTGAGGGCAAACGCATCGGTTTCGTACCGACTATGGGCAACCTGCACAGCGGCCATGCGGCCTTGGTGACCAAGGCTGCCCAGCGCGCCGACTTCGTGGTCGCCAGTATCTTCGTCAACCCGCTGCAGTTCGGCGCCAACGAAGACCTCGACAAATACCCACGCACCCTGGCAGCCGACCAGGAGCGCCTGGTGCAGGCCGGCTGCAACCTGCTGTT  >  NC_002947/5341713‑5342010

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑