Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,934,036	A→T	14.1%	L474Q (CTG→CAG)	PP_3467 ←	sigma‑54 dependent transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,934,036	0	A	T	14.1%	165.1 / 18.3	64	L474Q (CTG→CAG)	PP_3467	sigma‑54 dependent transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (25/30);  new base T (4/5);  total (29/35)
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.21e-01

GCAGGCCCTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATGGGCCTGCATCAGCCGC  >  NC_002947/3933890‑3934185                                                                                                                                                   |                                                                                                                                                      gcagGCCCTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAgcg                                                                                                                                                    >  2:956891/1‑150 (MQ=255) gcagGCCCTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAgcg                                                                                                                                                    <  2:1409562/150‑1 (MQ=255)    ggCCCTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCTGATCTGCGCCa                                                                                                                                                 <  1:650776/150‑1 (MQ=255)         tGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGATGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCg                                                                                                                                            <  1:1402882/150‑1 (MQ=255)          ggATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGa                                                                                                                                           >  2:972239/1‑150 (MQ=255)           gatgatCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGAt                                                                                                                                          <  1:618827/150‑1 (MQ=255)                tCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTc                                                                                                                                             <  1:1132728/142‑1 (MQ=255)                tCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTc                                                                                                                                             >  2:1132728/1‑142 (MQ=255)                tCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCt                                                                                                                                     >  1:496462/1‑150 (MQ=255)                      ccAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGt                                                                                                                               >  2:568305/1‑150 (MQ=255)                       cAGGTCGTCGCGTTCACGCAGGGCCGGCCGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTc                                                                                                                              >  1:1256465/1‑150 (MQ=255)                           gcgtcgCGTTCACGCAGGGCCGGCAGGCGCAGCGAGACACCGTGGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTAGATGTCCTGCTCGTCAccc                                                                                                                          <  2:231432/149‑1 (MQ=255)                              tcgCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGcc                                                                                                                       >  2:1194458/1‑150 (MQ=255)                                gcgTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGcccc                                                                                                                     >  1:1366855/1‑150 (MQ=255)                                     cACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAgg                                                                                                                <  1:637677/150‑1 (MQ=255)                                        gCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCTGATCTGCGCCa                                                                                                                                                 >  2:63281/1‑114 (MQ=255)                                        gCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCTGATCTGCGCCa                                                                                                                                                 <  1:63281/114‑1 (MQ=255)                                               cGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTAcgcg                                                                                                      <  2:1366855/150‑1 (MQ=255)                                               cGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTAcgcg                                                                                                      <  1:940432/150‑1 (MQ=255)                                                      cgcAGCGACACACCGTTACCTCGGCAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACGAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTAGATGTCCTGCTCGGCACCCGCCCCCAAGGGTGCTACGCGGCGCTcc                                                                                               <  2:566763/150‑1 (MQ=255)                                                       gcAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCt                                                                                              >  2:1317555/1‑150 (MQ=255)                                                          gCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCa                                                                                           <  1:1208070/150‑1 (MQ=255)                                                          gCGACACACCGTTGACCCGGGAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCa                                                                                           <  1:718938/150‑1 (MQ=255)                                                            gACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGt                                                                                         >  1:512918/1‑150 (MQ=255)                                                                       gACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAACGCCACGTCCATGTCCTGCTCGTCCCCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAg                                                                              <  1:357986/150‑1 (MQ=255)                                                                         ccgatctgAACAGGTCTTCACGGAAGTGCTCCTCATGGACCAGTCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGc                                                                            <  2:1045720/143‑1 (MQ=255)                                                                         cccGGTAGAACAGGTCTTCACGGAAGTGCTGCCCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTAGCCGTAACAGGc                                                                            <  2:1304007/150‑1 (MQ=255)                                                                          ccGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACGAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCg                                                                           <  2:1256465/150‑1 (MQ=255)                                                                          ccGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGc                                                                                                            >  2:95406/1‑117 (MQ=255)                                                                          ccGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGc                                                                                                            <  1:95406/117‑1 (MQ=255)                                                                          ccGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGc                                                                                            >  1:460890/1‑133 (MQ=255)                                                                          ccGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGc                                                                                            <  2:460890/133‑1 (MQ=255)                                                                           cGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGt                                                                          >  1:107540/1‑150 (MQ=255)                                                                                gAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCTGATCTGCGCCa                                                                                                                                                 >  2:254942/1‑74 (MQ=39)                                                                                gAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCTGATCTGCGCCa                                                                                                                                                 <  1:254942/74‑1 (MQ=39)                                                                                 aaCAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCaa                                                                                                                  <  1:853334/104‑1 (MQ=255)                                                                                 aaCAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCaa                                                                                                                  >  2:853334/1‑104 (MQ=255)                                                                                 aaCAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCgagatcgg                                                                    >  1:1045720/1‑143 (MQ=255)                                                                                  aCAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGAt                                                                                                                                          <  2:1223707/79‑1 (MQ=255)                                                                                  aCAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGAt                                                                                                                                          >  1:1223707/1‑79 (MQ=255)                                                                                  aCAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTAcgcg                                                                                                      <  2:279488/115‑1 (MQ=255)                                                                                  aCAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTAcgcg                                                                                                      >  1:279488/1‑115 (MQ=255)                                                                                  aCAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCa                                                                   <  1:1194458/150‑1 (MQ=255)                                                                                      gTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGc                                                               <  1:1317555/150‑1 (MQ=255)                                                                                        cTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCTATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGa                                                             >  1:1389841/1‑150 (MQ=255)                                                                                                     ctgctCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGcc                                                <  2:1169364/150‑1 (MQ=255)                                                                                                         tCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGAtt                                            >  2:582552/1‑150 (MQ=255)                                                                                                               aCCAGGCGCTTGAGGTCGCGGTGGATGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTtcgtcg                                      >  2:80499/1‑150 (MQ=255)                                                                                                                 cAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAg                                    <  2:496462/150‑1 (MQ=255)                                                                                                                 cAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAg                                    <  1:1210151/150‑1 (MQ=255)                                                                                                                          gAGGTCGCGGTGGGTTGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGGACAGCGt                           <  1:838937/150‑1 (MQ=255)                                                                                                                                 cgGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTg                                                                                             >  2:81097/1‑77 (MQ=255)                                                                                                                                 cgGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTg                                                                                             <  1:81097/77‑1 (MQ=255)                                                                                                                                        tggCGCTGATCTGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGc                                                           >  2:616529/1‑104 (MQ=255)                                                                                                                                        tggCGCTGATCTGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGc                                                           <  1:616529/104‑1 (MQ=255)                                                                                                                                        tggCGCTGATCTGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTtcgtc                                       >  2:118178/1‑124 (MQ=255)                                                                                                                                        tggCGCTGATCTGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTtcgtc                                       <  1:118178/124‑1 (MQ=255)                                                                                                                                          gcgcAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATGGGCCTGc           >  2:707769/1‑150 (MQ=255)                                                                                                                                            gcAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATGGGCCTGCAt         >  1:1279926/1‑150 (MQ=255)                                                                                                                                              aGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATGGGCCTGCATCa       <  2:107540/150‑1 (MQ=255)                                                                                                                                              aGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATGGGCCTGCATCa       <  2:1279926/150‑1 (MQ=255)                                                                                                                                                 tCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCa                                                                                           <  1:44505/63‑1 (MQ=255)                                                                                                                                                 tCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCa                                                                                           >  2:44505/1‑63 (MQ=255)                                                                                                                                                 acAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATGGGCCTGCATCAgcc    >  2:978027/2‑150 (MQ=255)                                                                                                                                                   aGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATGGGCCTGCATCAgccgc  >  2:692233/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                      GCAGGCCCTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATGGGCCTGCATCAGCCGC  >  NC_002947/3933890‑3934185

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