Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	968,878	G→T	8.3%	intergenic (‑308/+23)	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,878	0	G	T	8.3%	191.2 / 4.4	60	intergenic (‑308/+23)	PP_0831/PP_t68	hypothetical protein/tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base G (29/26);  new base T (2/3);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

AACAATACTTTCGGAATATATCTGACGAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGCCGCTGGTCGGT  >  NC_002947/968730‑969020                                                                                                                                                      |                                                                                                                                                aaCAATAATTTCGGAATATATCTGACGAATGACAGCTTAATGTCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGt                                                                                                                                                <  1:382797/150‑1 (MQ=255)  aCAATACTTCCGGAATATATCTAAAGAATTCAATCATCATTGCTATCGCAAGAGAAGGCCAATAGAATGTAGGTATGGGGGGAGGCGCGGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTa                                                                                                                                               <  2:1873835/150‑1 (MQ=255)         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                                                                                                                                     >  2:1424989/1‑148 (MQ=255)              gAATATATCTGACGAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAGAAAAACCGACCTTACGGTCGGttt                                                                                                                                     <  1:1424989/148‑1 (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                                                                                                              >  1:1370558/1‑130 (MQ=255)                                                      gggCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAg                                                                                                              <  2:1370558/130‑1 (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                                                                        <  2:963510/150‑1 (MQ=255)                                                                        aGGTATGGGGGGAGGCGCGGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGAtttt                                                                        >  1:562751/1‑150 (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                                  >  1:1691873/3‑125 (MQ=255)                                                                                                                                        ggaaaaaaCCGACCTTACGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTg                                  <  2:1691873/123‑1 (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                           <  1:114955/122‑1 (MQ=255)                                                                                                                                                  gaccgTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCa                           >  2:114955/1‑122 (MQ=255)                                                                                                                                                      gTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTg                                  >  2:87803/1‑111 (MQ=255)                                                                                                                                                      gTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTg                                  <  1:87803/111‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                AACAATACTTTCGGAATATATCTGACGAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCA‑GGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCG‑GATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGCCGCTGGTCGGT  >  NC_002947/968730‑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: ‑