Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	968,881	A→C	7.9%	intergenic (‑311/+20)	PP_0831 ← / ← PP_t68	hypothetical protein/tRNA‑Leu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	968,881	0	A	C	7.9%	167.1 / 6.7	63	intergenic (‑311/+20)	PP_0831/PP_t68	hypothetical protein/tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base A (31/27);  new base C (2/3);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

TTTCGGAATATATCTGACGAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGCCGCTGGTCGGT  >  NC_002947/968738‑969020                                                                                                                                                 |                                                                                                                                             tttCGGAATATATCTGACGAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGCAA‑GACGGGCAGAGTCAAAGAGGGGCGTGAAGCGCAGGCACAAAAAAAACGACCGTAAGGCCgg                                                                                                                                        >  2:966855/1‑150 (MQ=255) ttccGGAATATATCGGACGAATGACAGCTTAAGTGCTAGGTGAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCgg                                                                                                                                        <  1:966855/147‑1 (MQ=255)   tCGGAATATATCTGACGAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCAAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATAGCAGAAACAAAAAAACCGACCGTAAGGTCGGtt                                                                                                                                      <  2:548665/150‑1 (MQ=255)      gAATATATCTGACGAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAGAAAAACCGACCTTACGGTCGGttt                                                                                                                                     <  1:1424989/148‑1 (MQ=255)      gAATATATCTGACGAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAGAAAAACCGACCTTACGGTCGGttt                                                                                                                                     >  2:1424989/1‑148 (MQ=255)      gAATATATCTGACGAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCTTACGGTCGGttttt                                                                                                                                   <  2:2336054/150‑1 (MQ=255)            atCTGACGAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCAAAAAAAAACCGACCGTAAGGTCGGTTTTTTCCg‑g                                                                                                                             >  2:904371/1‑150 (MQ=255)                             ttAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGga                                                                                                           >  2:1294051/1‑150 (MQ=255)                                 ttGCTACCGCAAGGCGCTGCCAAGAGAATGCAGGTATGGGGGGAGGCGCGGTAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCCAATCTACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACt                                                                                                       <  2:2069831/149‑1 (MQ=255)                                       gcgcAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGT‑GATCGATAATGTGATTTGGTGCCCCGAGGGAGACTCGAACt                                                                                                 >  2:1355223/1‑150 (MQ=255)                                        cgcAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTc                                                                                                <  2:110528/150‑1 (MQ=255)                                          cAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTccc                                                                                              >  1:963510/1‑150 (MQ=255)                                              gggCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAg                                                                                                              <  2:1370558/130‑1 (MQ=255)                                              gggCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAg                                                                                                              >  1:1370558/1‑130 (MQ=255)                                               ggcggcCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTcc                                                                                         <  2:991671/150‑1 (MQ=255)                                                ggggcCAAGAGAATGCAGGTATGGGGGGAGGCGCGGGAGTAA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCt                                                                                        <  2:114802/148‑1 (MQ=255)                                                      aaTAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGaa                                                                                  >  2:224936/1‑150 (MQ=255)                                                        tAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGaaaa                                                                                >  2:867188/1‑150 (MQ=255)                                                        tAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGaaaa                                                                                >  1:1159297/1‑150 (MQ=255)                                                         aGAATGCAGGTATGGGGGTAGGCGCGGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAAc                                                                               >  2:2273247/1‑150 (MQ=255)                                                                agttATGGGGGGAGGCGCGGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GAttt                                                                                                                         <  1:55162/98‑1 (MQ=255)                                                                aGGTATGGGGGGAGGCGCGGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GAttt                                                                                                                         >  2:55162/1‑101 (MQ=255)                                                                aGGTATGGGGGGAGGCGCGGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGAtttt                                                                        <  1:1294051/150‑1 (MQ=255)                                                                aGGTATGGGGGGAGGCGCGGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGAtttt                                                                        >  1:562751/1‑150 (MQ=255)                                                                aGGTATGGGGGGAGGCGCGGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGAtttt                                                                        <  2:963510/150‑1 (MQ=255)                                                                aGGTAGGGGGGGAGGCGCGGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGAtttt                                                                        <  1:867188/150‑1 (MQ=255)                                                                  gTATGGGGGGAGGCGCGGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGa                                                                      >  2:1069996/1‑150 (MQ=255)                                                                   tATGGGGGGAGGCGCGGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGaa                                                                     <  1:1069996/150‑1 (MQ=255)                                                                   ggtGGGGGTAGGGGCGGGAGTCA‑GGGGCGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGaa                                                                     <  1:1682602/148‑1 (MQ=255)                                                                           gAGGGGCGGAAGTCA‑GGAGCGCGGCGTCAACGACGTGCGTGAACCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGCGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCcgccgc                                                             <  2:66377/150‑1 (MQ=255)                                                                           gAGGCGCGGGAGTCA‑GGGGGGGGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCcgccgc                                                             <  1:1127562/150‑1 (MQ=255)                                                                                gcgGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTa                                                        >  2:82943/1‑150 (MQ=255)                                                                                            gggggCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCGG‑GATTTGGTGCCCCGAGGGAGACTCGAAATCCCACTCCTTTTGAAAACGGATTTTGAACCCGCCGCGTCTACCAATTCCGcc                                             >  1:2276686/1‑150 (MQ=255)                                                                                              gggCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCAt                                           >  2:1930088/1‑150 (MQ=255)                                                                                                  ggcgTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAAc                                                                               >  1:521384/1‑110 (MQ=255)                                                                                                  ggcgTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAAc                                                                               <  2:521384/110‑1 (MQ=255)                                                                                                        aaCGACGTGCGTGAATCGCAGGCACAAAAAAACCGACGGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTgg                                 >  1:1887791/1‑150 (MQ=255)                                                                                                        aaCGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGTTTTTTTCCG‑GATTTGGTGCCCCAACGGAGACTCGAACTCCCACTCCTTTCGAAAAGGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCGGGGCGTgg                                 >  2:1006659/1‑150 (MQ=255)                                                                                                         acgacgTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGgc                                <  1:82943/150‑1 (MQ=255)                                                                                                         acgacgTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGgc                                >  2:1949013/1‑150 (MQ=255)                                                                                                           gacgTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGgcgc                              >  2:1340918/1‑150 (MQ=255)                                                                                                               tgcgtgAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAg                          <  1:1930088/150‑1 (MQ=255)                                                                                                                 cgtgAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGAGTCGACCAATTCCGCCATCAGGGCGTGGCGCGCAGta                        <  1:1355223/150‑1 (MQ=255)                                                                                                                     aaTCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGaa                                                                                                   >  1:1010337/1‑71 (MQ=255)                                                                                                                     aaTCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGaa                                                                                                   <  2:1010337/71‑1 (MQ=255)                                                                                                                      aTCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAag                   <  1:1949013/150‑1 (MQ=255)                                                                                                                         gCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAgg                >  1:1770078/1‑150 (MQ=255)                                                                                                                            ggCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGc             <  1:1340918/150‑1 (MQ=255)                                                                                                                                ggaaaaaaCCGACCTTACGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTg                                  <  2:1691873/123‑1 (MQ=255)                                                                                                                                ggaaaaaaCCGACCTTACGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTg                                  >  1:1691873/3‑125 (MQ=255)                                                                                                                                caAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAAc                                                                               >  2:265295/1‑80 (MQ=255)                                                                                                                                caAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAAc                                                                               <  1:265295/80‑1 (MQ=255)                                                                                                                                 aaaaaaaCCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGCCGCTg        >  2:627410/1‑150 (MQ=255)                                                                                                                                 aaaaaaaCCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGCCGCTg        >  2:971634/1‑150 (MQ=255)                                                                                                                                   aaaaaCCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGaaaa                                                                                >  1:2051673/1‑76 (MQ=255)                                                                                                                                   aaaaaCCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGaaaa                                                                                <  2:2051673/76‑1 (MQ=255)                                                                                                                                   aaaaaCCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGCCGCTggt      >  1:1280169/1‑150 (MQ=255)                                                                                                                                      aaCCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGCCGCTggtcgg   >  1:417816/1‑150 (MQ=255)                                                                                                                                       accgaccgTAAGGTCGGTTTTGCCCG‑GATTTGGTGTCCTGAGGGAGACTCGAACGCCCACTCCTTTCGAACACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGCCGCTggtcggt  <  1:194412/150‑1 (MQ=255)                                                                                                                                          gaccgTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCa                           >  2:114955/1‑122 (MQ=255)                                                                                                                                          gaccgTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCa                           <  1:114955/122‑1 (MQ=255)                                                                                                                                              gTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTg                                  >  2:87803/1‑111 (MQ=255)                                                                                                                                              gTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTg                                  <  1:87803/111‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                             TTTCGGAATATATCTGACGAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGCCGCTGGTCGGT  >  NC_002947/968738‑969020

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