Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,013,339	G→A	9.9%	A38T (GCG→ACG)	PP_5678 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,013,339	0	G	A	9.9%	152.3 / 8.6	61	A38T (GCG→ACG)	PP_5678	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (26/28);  new base A (3/3);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.87e-01

GCAGTGCAACAAGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGACTCCTTACTCTC  >  NC_002947/5013193‑5013484                                                                                                                                                   |                                                                                                                                                  gCAGTGCAACAAGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCgcgg                                                                                                                                                >  2:1437601/1‑150 (MQ=255)        aacaaGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGaa                                                                                                                                         >  2:788598/1‑150 (MQ=255)         acaaGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGCGCTCAGCACGGCGTgctgt                                                                                                                                        <  1:750293/150‑5 (MQ=255)               aacAGAGTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGACATCCACAGAACGACGTGCTGAGCGCGGGGTGGAAATGAGGc                                                                                                                                  >  1:645246/3‑150 (MQ=255)                tCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCt                                                                                                                                 <  1:788598/150‑1 (MQ=255)                          gggAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCaaaa                                                                                                                       >  2:213802/1‑150 (MQ=255)                          gggAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCaaaa                                                                                                                       >  1:513700/1‑150 (MQ=255)                          gggAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGCGCTCAGCacggcg                                                                                                                                              >  2:881382/1‑127 (MQ=255)                          gggAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGCGCTCAGCacggcg                                                                                                                                              <  1:881382/127‑1 (MQ=255)                                 ccAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGaa                                                                                                                >  1:77601/1‑150 (MQ=255)                                      gatgatCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTa                                                                                                           >  2:642446/1‑150 (MQ=255)                                      gatgatCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGCGCTCAGCACGGCGTgctgt                                                                                                                                        <  1:104014/121‑5 (MQ=255)                                      gatgatCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGCGCTCAGCACGGCGTgctgt                                                                                                                                        >  2:518281/1‑117 (MQ=255)                                      gatgatCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGCGCTCAGCACGGCGTgctgt                                                                                                                                        >  2:104014/1‑117 (MQ=255)                                      gatgatCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGCGCTCAGCACGGCGTgctgt                                                                                                                                        <  1:518281/121‑5 (MQ=255)                                           tCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTg                                                                                                                                            >  2:599918/1‑112 (MQ=255)                                           tCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTg                                                                                                                                            <  1:599918/112‑1 (MQ=255)                                           tCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAAc                                                                                                      >  1:458811/1‑150 (MQ=255)                                             aaaGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTc                                                                                                    <  2:77601/150‑1 (MQ=255)                                               agagTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCa                                                                                                  >  1:724760/1‑150 (MQ=255)                                                 agTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATc                                                                                                <  1:1437601/150‑1 (MQ=255)                                                  gTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCAGCCACAGCACGCCGTGCTACCCTTTGCATTGGAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCt                                                                                               >  1:774461/1‑150 (MQ=255)                                                     gACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTGGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTg                                                                                            >  2:305805/1‑150 (MQ=255)                                                      aCCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGc                                                                                           <  2:513700/150‑1 (MQ=255)                                                        ccTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCga                                                                                         <  1:1556934/150‑1 (MQ=255)                                                          tCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCgcag                                                                                                                            >  1:304151/1‑113 (MQ=255)                                                          tCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCgcag                                                                                                                            <  2:304151/113‑1 (MQ=255)                                                              aGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTg                                                                                                                                            >  1:1243369/1‑93 (MQ=255)                                                              aGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTg                                                                                                                                            <  2:1243369/93‑1 (MQ=255)                                                                  gTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGcttc                                                                               >  1:448277/1‑150 (MQ=255)                                                                       ggTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGc                                                                          <  1:642446/150‑1 (MQ=255)                                                                         tCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGaa                                                                                                                <  1:83788/110‑1 (MQ=255)                                                                         tCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGaa                                                                                                                >  2:83788/1‑110 (MQ=255)                                                                                           gTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGc                                                      <  1:213802/150‑1 (MQ=255)                                                                                            ttGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCTGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCg                                                     <  2:387911/150‑1 (MQ=255)                                                                                            ttGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCGGAGCGCTGCGTGGAATTCAGGCGCGCAGGAAAACCCCGAAGACTATGAACTCCATCTTTGCGAGAGCTGGTTCTTTGCCACGCTTTCAACGATAAAGCg                                                     <  2:387716/150‑1 (MQ=255)                                                                                             tGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGt                                                    <  2:448277/150‑1 (MQ=255)                                                                                                ggCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCGTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGaa                                                 >  1:177891/1‑150 (MQ=255)                                                                                                  ccAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGcttctt                                                                             <  2:1091019/120‑1 (MQ=255)                                                                                                  ccAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGcttctt                                                                             >  1:1091019/1‑120 (MQ=255)                                                                                                   cAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGc                                              <  2:458811/150‑1 (MQ=255)                                                                                                   cAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGc                                              >  1:1390348/1‑150 (MQ=255)                                                                                                   cAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGc                                              >  1:1286993/1‑150 (MQ=255)                                                                                                           tCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAg                                                       <  1:633142/133‑1 (MQ=255)                                                                                                           tCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAg                                                       >  2:633142/1‑133 (MQ=255)                                                                                                              ttCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATaaa                                                        >  2:1492466/1‑129 (MQ=255)                                                                                                              ttCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATaaa                                                        <  1:1492466/129‑1 (MQ=255)                                                                                                              ttCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTggg                                          <  2:1589664/143‑1 (MQ=255)                                                                                                              ttCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTggg                                          >  1:1589664/1‑143 (MQ=255)                                                                                                              ttCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAAtt                                   >  2:1033689/1‑150 (MQ=255)                                                                                                                  gatgaCGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATg                               >  1:208316/1‑150 (MQ=255)                                                                                                                   atgaCGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGt                              >  2:1115862/1‑150 (MQ=255)                                                                                                                      aCGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTtct                           <  2:774461/150‑1 (MQ=255)                                                                                                                          ccATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTGGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGa                       <  1:305805/150‑1 (MQ=255)                                                                                                                          ccATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGa                       <  1:1115862/150‑1 (MQ=255)                                                                                                                               cacaGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAg                  >  2:1461069/1‑150 (MQ=255)                                                                                                                                 caGCACGCCGTGCTGAGCGCGTCGTGGAAATCAGGCGCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGGGAACGCTGGGCCCAATTGCTGTTCTCTGAGGAAGCt                <  2:839730/150‑1 (MQ=255)                                                                                                                                 caGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCt                <  1:1033689/150‑1 (MQ=255)                                                                                                                                 caGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCt                <  2:1390348/150‑1 (MQ=255)                                                                                                                                 caGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCt                <  1:1461069/150‑1 (MQ=255)                                                                                                                                   gCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGa              >  1:1200010/1‑150 (MQ=255)                                                                                                                                    cACGCCGCGCTCAGCACGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGAc             >  1:1571032/1‑150 (MQ=255)                                                                                                                                            gcttatcgccgcGTGGAAATCAGGCTCGCACCAAAACCCCGAAGCCTATTAACTCCATCTTTGCGAGATCTGCTTCTTTTCCACACTTTCAACTATAAAGCTTGAACGCTTGGCCCAATTGATGTTCTTTGAGGTATCTGTCTCCTTAct     <  2:645246/144‑1 (MQ=255)                                                                                                                                              tGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGACTCCTTActct   <  2:724760/150‑1 (MQ=255)                                                                                                                                               gAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGACTCCTTACtctc  <  1:948874/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  GCAGTGCAACAAGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGACTCCTTACTCTC  >  NC_002947/5013193‑5013484

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