Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,506,377	C→A	5.8%	A5S (GCG→TCG)	PP_5574 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,506,377	0	C	A	5.8%	225.4 / 2.2	69	A5S (GCG→TCG)	PP_5574	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/32);  new base A (2/2);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

ACGGGTATAGGAATCAGACTTGAAACAATACCGCTCATCATTGATCCAGCACACCTGCCAGGCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCAAAAGATGGCGCAATGCCATCTCATTCATGACGGCTAATGCTAAGATGAGCGATATTCATAGTTGACAGAGTTGCGCATTCC  >  NC_002947/3506231‑3506518                                                                                                                                                   |                                                                                                                                              aCGGGTATAGGAATCAGACTTGAAACAATACCGCTCATCATTGATCCAGCACACCTGCCAGGCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTAt                                                                                                                                            >  2:776917/1‑150 (MQ=255) aCGGGTATAGGAATCAGACTTGAAACAATACCGCTCATCATTGATCCAGCACACCTGCCAGGCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTAt                                                                                                                                            >  1:340323/1‑150 (MQ=255) aCGGGTATAGGAATCAGACTTGAAACAATACCGCTCATCATTGATCCAGCACACCTGCCAGGCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTAt                                                                                                                                            >  2:392196/1‑150 (MQ=255)      tataGGAATCAGACTTGAAACAATACCGCTCATCATTGATCCAGCACACCTGCCAGGCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGa                                                                                                                                       <  2:673072/150‑1 (MQ=255)       ataGGAATCAGACTTGAAACAATACCGCTCATCATTGATCCAGCACACCTGCCAGGCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACGTGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGaa                                                                                                                                      <  2:155198/150‑1 (MQ=255)        taGGAATCAGACTTGAAACAATACCGCTCATCATTGATCCAGCACACCTGCCAGGCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAAc                                                                                                                                     >  1:1115220/1‑150 (MQ=255)          ggAATCAGACTTGAAACAATACCGCTCATCATTGATCCAGCACACCTGCCAGGCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACAt                                                                                                                                   <  1:528337/150‑1 (MQ=255)           gAATCAGACTTGAAACAATACCGCTCATCATTGATCCAGCACACCTGCCAGGCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACAtt                                                                                                                                  >  1:289630/1‑150 (MQ=255)             aTCAGACTTGAAACAATACCGCTCATCATTGATCCAGCACACCTGCCAGGCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGaa                                                                                                                                      >  1:1359071/1‑144 (MQ=255)             aTCAGACTTGAAACAATACCGCTCATCATTGATCCAGCACACCTGCCAGGCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGaa                                                                                                                                      <  2:1359071/144‑1 (MQ=255)             aTCAGACTTGAAACAATACCGCTCATCATTGATCCAGCACACCTGCCAGGCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGACCATTGa                                                                                                                                <  2:478513/150‑1 (MQ=255)                 gACTTGAAACAATACCGCTCATCATTGATCCAGCACACCTGCCAGGCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACcac                                                                                                                            >  2:1131071/1‑150 (MQ=255)                    ttGAAACAATACCGCTCATCATTGATCCAGCACACCTGCCAGGCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAgt                                                                                                                         >  2:888815/1‑150 (MQ=255)                       aaaCAATACCGCTCATCATTGATCCAGCACACCTGCCAGGCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTc                                                                                                                      >  1:1274384/1‑150 (MQ=255)                          cAATACCGCTCATCATTGATCCAGCACACCTGCCAGGCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCtg                                                                                                                   <  2:472750/150‑1 (MQ=255)                              aCCGCTCATCATTGATCCAGCACACCTGCCAGGCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCCGTGaa                                                                                                               <  2:1274384/150‑1 (MQ=255)                                   tcatcaTTGATCCAGCACACCTGCCAGGCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTc                                                                                                          <  2:1115220/150‑1 (MQ=255)                                    catcatTGATCCAGCACACCTGCCAGGCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCt                                                                                                         <  1:220746/150‑1 (MQ=255)                                     atcatTGATCCAGCACACCTGCCAGGCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCtt                                                                                                        >  1:1475240/1‑150 (MQ=255)                                        atTGATCCAGCACACCTGCCAGGCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGt                                                                                                     >  2:19169/1‑150 (MQ=255)                                          tGATCCAGCACACCTGCCAGGCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTggcgg                                                                                                   <  1:875762/150‑4 (MQ=255)                                            aTCCAGCACACCTGCCAGGCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCtgt                                                                                                                  <  1:1030987/133‑1 (MQ=255)                                            aTCCAGCACACCTGCCAGGCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCtgt                                                                                                                  >  2:1030987/1‑133 (MQ=255)                                             tCCAGCACACCTGCCAGGCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTc                                                                                                <  2:1080803/150‑1 (MQ=255)                                               cAGCACACCTGCCAGGCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTg                                                                                              <  2:289630/150‑1 (MQ=255)                                                 gCACACCTGCCAGGCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGgt                                                                                            <  1:19169/150‑1 (MQ=255)                                                     acCTGCCAGGCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTAt                                                                                        <  2:1131713/150‑1 (MQ=255)                                                     acCTGCCAGGCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTAt                                                                                        >  2:1234129/1‑150 (MQ=255)                                                        tGCCAGGCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACg                                                                                                             >  2:520739/1‑126 (MQ=255)                                                        tGCCAGGCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACg                                                                                                             <  1:520739/126‑1 (MQ=255)                                                            aGGCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCaaa                                                                                 >  2:1111451/1‑150 (MQ=255)                                                              gCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCAAAAg                                                                               >  1:1072179/1‑150 (MQ=255)                                                                  aTCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCAAAAGATgg                                                                           <  2:340323/150‑1 (MQ=255)                                                                           tCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCAAAAGATGGCGCAATGcc                                                                  >  2:29430/1‑150 (MQ=255)                                                                             aaCAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGt                                                                                                     >  1:305201/1‑113 (MQ=255)                                                                             aaCAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGt                                                                                                     <  2:305201/113‑1 (MQ=255)                                                                             aaCAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTAt                                                                                        <  1:113099/126‑1 (MQ=255)                                                                             aaCAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTAt                                                                                        >  2:113099/1‑126 (MQ=255)                                                                             aaCAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCAAAAGATGGCGCAATGCCAt                                                                >  1:1211993/1‑150 (MQ=255)                                                                             aaCAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCAAAAGATGGCGCAATGCCAt                                                                >  2:938447/1‑150 (MQ=255)                                                                               cAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCAAAAGATGGCGCAATGCCAtct                                                              <  1:776917/150‑1 (MQ=255)                                                                                  cTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGAATCCCAAGCGTAGACGCTATGTAGAACATTGACCACAGTGTCCTCTGAACGTTCTTGGTAATTCTGGTTTATGGGCAAAAGATGGCGCAATGCCATCgcat                                                           >  1:707275/1‑150 (MQ=255)                                                                                   ttCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCAAAAGATGGCGCAATGCCATCtcatt                                                          <  1:392196/150‑1 (MQ=255)                                                                                      tGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCAAAAGATGGCGCAATGCCATCtcattcat                                                       >  1:302338/1‑150 (MQ=255)                                                                                                  tACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGGACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCAAAAGATGGCGCAATGCCATCTCATTCATGACGGCTAATGc                                           >  1:1293252/1‑150 (MQ=255)                                                                                                    ccccTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCAAAAGATGGCGCAATGCCa                                                                 >  2:649504/1‑126 (MQ=255)                                                                                                    ccccTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCAAAAGATGGCGCAATGCCa                                                                 <  1:649504/126‑1 (MQ=255)                                                                                                     cccTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCAAAAGATgg                                                                           <  1:1064749/115‑1 (MQ=255)                                                                                                     cccTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCAAAAGATgg                                                                           >  2:1064749/1‑115 (MQ=255)                                                                                                            cccACAAGCCCAGCGTCACAAGTATTACAACCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTTGCAATTCTGGTGTATGGGCAAAAGATGGCGCAATGCCATCTCATTCATGACGGCTAATGCTAAGATGAGc                                 <  2:707275/150‑1 (MQ=255)                                                                                                             ccACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCAAAAGATGGCGCAATGCCATCTCATTCATGACGGCTAATGCTAAGATGAGCg                                >  2:863107/1‑150 (MQ=255)                                                                                                              cacaAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCAAAAGATGGCGCAATGCCAt                                                                <  2:1019402/117‑1 (MQ=255)                                                                                                              cacaAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCGTGGTAATTCTGGTGTATGGGCAAAAGATGGCGCAATGCCAt                                                                >  1:1019402/1‑117 (MQ=255)                                                                                                               acaAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCAAAAGATGGCGCAATGCCATCTCATTCATGACGGCTAATGCTAAGATGAGCGat                              >  2:645695/1‑150 (MQ=255)                                                                                                                    cccAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCAAAAGATGGCGCAATGCCATCTCATTCATGACGGCTAATGCTAAGATGAGCGATATTCa                         <  1:888815/150‑1 (MQ=255)                                                                                                                           tCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCAAAAGATGGCGCAATGCCATCtcatt                                                          >  1:316728/1‑110 (MQ=255)                                                                                                                           tCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCAAAAGATGGCGCAATGCCATCtcatt                                                          <  2:316728/110‑1 (MQ=255)                                                                                                                            cacaAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCAAAAGATGGCGCAATGCCATCTCATTCATGACGGCTAAt                                             <  1:493173/122‑1 (MQ=255)                                                                                                                            cacaAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCAAAAGATGGCGCAATGCCATCTCATTCATGACGGCTAAt                                             >  2:493173/1‑122 (MQ=255)                                                                                                                              caAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCAAAAGATGGCGCAATGCCATCTCATTCATGACGGCTAATGCTAAGATGAGCGATATTCATAGTTGACag               >  1:1494716/1‑150 (MQ=255)                                                                                                                               aaGTATCCCAAGCGTAGGCGATATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCAAAAGATGGCGCAATGCCATCtcatt                                                          >  2:531744/1‑106 (MQ=255)                                                                                                                               aaGTATCCCAAGCGTAGGCGATATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCAAAAGATGGCGCAATGCCATCtcatt                                                          <  1:531744/106‑1 (MQ=255)                                                                                                                                 gTATCCCAAGCGTAGGCGATATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCAAAAGATGGCGCAATGCCATCTCATTCATGACGGCTAATGCTAAGATGAGCGATATTCATa                       <  2:595366/139‑1 (MQ=255)                                                                                                                                 gTATCCCAAGCGTAGGCGATATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCAAAAGATGGCGCAATGCCATCTCATTCATGACGGCTAATGCTAAGATGAGCGATATTCATa                       >  1:595366/1‑139 (MQ=255)                                                                                                                                  tATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCAAAAGATGGCGCAATGCCATCTCATTCATGACGGCTAATGCTAAGATGAGCGATATTCATAGTTGACAGAGtt           <  2:1211993/150‑1 (MQ=255)                                                                                                                                    tCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCAAAAGATGGCGCAATGCCATCTCATTCATGACGGCTAATGCTAAGATGAGCGATATTCATAGTTGACAGAGTTgc         <  1:1131071/150‑1 (MQ=255)                                                                                                                                    tCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCAAAAGATGGCGCAATGCCATCTCATTCATGACGGCTAATGCTAAGATGAGCGATATTCATAGTTGACAGAGTTgc         <  1:1234129/150‑1 (MQ=255)                                                                                                                                      ccAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCAAAAGATGGCGCAATGCCATCTCATTCATGACGGCTAATGCTAAGATGAGCGATATTCATAGTTGACAGAGTTgcgc       >  2:1504515/1‑150 (MQ=255)                                                                                                                                           cGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCAAAAGATGGCGCAATGCCATCTCATTCATGACGGCTAATGCTAAGATGAGCGATATTCATAGTTGACAGAGTTGCGCATTcc  <  1:1504515/150‑1 (MQ=255)                                                                                                                                                   cTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGTTGTATGGGCAAAAGATGGCGCAATGCCATCTCATTCATGACGGCTAATGCTAAGATGAGCGATATTCATAg                      >  1:976630/1‑122 (MQ=255)                                                                                                                                                   cTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGTTGTATGGGCAAAAGATGGCGCAATGCCATCTCATTCATGACGGCTAATGCTAAGATGAGCGATATTCATAg                      <  2:976630/122‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                              ACGGGTATAGGAATCAGACTTGAAACAATACCGCTCATCATTGATCCAGCACACCTGCCAGGCAGATCCTGAAGTCAACAGCTTCTGGCAGGGTGTGTACCCCTTCGCCCACAAGCCCAGCGTCACAAGTATCCCAAGCGTAGGCGCTATGTAGAACATTGACCACAGTGTCCTGTGAACGTTCTTGGTAATTCTGGTGTATGGGCAAAAGATGGCGCAATGCCATCTCATTCATGACGGCTAATGCTAAGATGAGCGATATTCATAGTTGACAGAGTTGCGCATTCC  >  NC_002947/3506231‑3506518

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