Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,625,483	A→C	9.2%	L110R (CTG→CGG) ‡	PP_4943 ←	glycosyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,625,483	0	A	C	9.2%	178.5 / 10.0	67	L110R (CTG→CGG) ‡	PP_4943	glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (34/26);  new base C (2/4);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 3.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

GGCAAGCAACGGCGCATATCGTCACGTACCGCGTCCGAAACACCCAACAGGCTCAGACGCTTGCTGAACAGGTTCGCGAACAGGCGGCGGCCTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGGCCGCAATGCGCCGCAATGCGCGGATGGCGCCGAGTTTCAGGCCACGCACGGCCTTGGAACTGAACTCAAGGAAC  >  NC_002947/5625345‑5625629                                                                                                                                           |                                                                                                                                                   ggCAAGCAACGGCGCATATCGTCACGTACCGCGTCCGAAACACCCAACAGGCTCAGACGCTTGCTGAACAGGTTCGCGAACAGGCGGCGGCCTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCCTGCCGGtgatt                                                                                                                                         <  2:656028/150‑5 (MQ=255)  gcaagcaaCGGCGCATATCGTCACGTACCGCGTCCGAAACACCCAACAGGCTCAGACGCTTGCTGAACAGGTTCGCGAACAGGCGGCGGCCTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCCTGCCGGtgattg                                                                                                                                        <  2:334478/150‑6 (MQ=255)  gcaagcaaCGGCGCATATCGTCACGTACCGCGTCCGAAACACCCAACAGGCTCAGACGCTTGCTGAACAGGTTCGCGAACAGGCGGCGGCCTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCCTGCCGGtga                                                                                                                                           <  1:742857/147‑3 (MQ=255)  gcaagcaaCGGCGCATATCGTCACGTACCGCGTCCGAAACACCCAACAGGCTCAGACGCTTGCTGAACAGGTTCGCGAACAGGCGGCGGCCTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCCTGCCGGtga                                                                                                                                           >  2:742857/1‑145 (MQ=255)    aagcaaCGGCGCATATCGTCACGTACCGCGTCCGAAACACCCAACAGGCTCAGACGCTTGCTGAACAGGTTCGCGAACAGGCGGCGGCCTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCCTGCCGGtgat                                                                                                                                          <  1:1088505/146‑4 (MQ=255)    aagcaaCGGCGCATATCGTCACGTACCGCGTCCGAAACACCCAACAGGCTCAGACGCTTGCTGAACAGGTTCGCGAACAGGCGGCGGCCTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCCTGCCGGtgat                                                                                                                                          >  2:1088505/1‑143 (MQ=255)    aagcaaCGGCGCATATCGTCACGTACCGCGTCCGAAACACCCAACAGGCTCAGACGCTTGCTGAACAGGTTCGCGAACAGGCGGCGGCCTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCcagc                                                                                                                                      >  1:530975/1‑150 (MQ=255)     agcaaCGGCGCATATCGTCACGTACCGCGTCCGAAACACCCAACAGGCTCAGACGCTTGCTGAACAGGTTCGCGAACAGGCGGCGGCCTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCcagca                                                                                                                                     <  1:986843/150‑1 (MQ=255)                        aCGTACCGCGTCCGAAACACCCAACAGGCTCAGACGCTTGCTGAACAGGTTCGCGAACAGGCGGCGGCCTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCtt                                                                                                                  >  2:26465/1‑150 (MQ=255)                             ccGCGTCCGAAACACCCAACAGGCTCAGACGCTTGCTGAACAGGTTCGCGAACAGGCGGCGGCCTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGc                                                                                                             >  1:477203/1‑150 (MQ=255)                                cgTCCGAAACACCCAACAGGCTCAGATGCTTGCTGAACAGGTTCGCGAACAGGCGGCGGCCTTTGAGCTATTAGTCGCCAAATGAATGGTGAACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGt                                                                                                          <  1:450163/150‑1 (MQ=255)                                cgTCCGAAACACCCAACAGGCTCAGACGCTTGCTGAACAGGTTCGCGAACAGGCGGCGGCCTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGt                                                                                                          <  1:1422317/150‑1 (MQ=255)                                    cGAAACACCCAACAGGCTCAGACGCTTGCTGAACAGGTTCGCGAACAGGCGGCGGCCTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGc                                                                                                      <  1:1342218/150‑1 (MQ=255)                                    cGAAACACCCAACAGGCTCAGACGCTTGCTGAACAGGTTCGCGAACAGGCGGCGGCCTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGc                                                                                                      >  1:1309665/1‑150 (MQ=255)                                          acCCAACAGGCTCAGACGCTTGCTGAACAGGTTCGCGAACAGGCGGCGGCCTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCa                                                                                                <  2:747052/150‑1 (MQ=255)                                          acCCAACAGGCTCAGACGCTTGCTGAACAGGTTCGCGAACAGGCGGCGGCCTTTGCGCTGGTAGTCGCCAAAGAAATCCTGCACCCCCATCACCACCTACCCGTGCGTGAACAGAGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCa                                                                                                <  1:26465/150‑1 (MQ=255)                                           cccAACAGGCTCAGACGCTTGCTGAACAGGTTCGCGAACAGGCGGCGGCCTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAg                                                                                               >  2:1026796/1‑150 (MQ=255)                                            ccAACAGGCTCAGACGCTTGCTGAACAGGTTCGCGAACAGGCGGCGGCCTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGCCTTGAAGCGGTGGGCGATGCAGa                                                                                              >  2:711845/1‑150 (MQ=255)                                               aCAGGCTCAGACGCTTGCTGAACAGGTTCGCGAACAGGCGGCGGCCTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAg                                                                                           <  2:530975/150‑1 (MQ=255)                                                   gCTCAGACGCTTGCTGAACAGGTTCGCGAACAGGCGGCGGCCTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGaa                                                                                       >  2:949160/1‑150 (MQ=255)                                                     tCAGACGCTTGCTGAACAGGTTCGCGAACAGGCGGCGGCCTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAAt                                                                                     >  1:407806/1‑150 (MQ=255)                                                        gACGCTTGCTGAACAGGTTCGCGAACAGGCGGCGGCCTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTgc                                                                                  >  2:874814/1‑150 (MQ=255)                                                                   aaCAGGTTCGCGAACAGGCGGCGGCCTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGgccgca                                                                       >  2:862516/1‑150 (MQ=255)                                                                   aaCAGGTTCGCGAACAGGCGGCGGCCTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGgccgca                                                                       >  1:332627/1‑150 (MQ=255)                                                                        gTTCGCGAACAGGCGGCGGCCTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTggg                                                                                                       <  2:754501/113‑1 (MQ=255)                                                                        gTTCGCGAACAGGCGGCGGCCTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGgccgcaatgcg                                                                  >  2:23897/1‑150 (MQ=255)                                                                        gTTCGCGAACAGGCGGCGGCCTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGgccgcaatgcg                                                                  >  2:1382288/1‑150 (MQ=255)                                                                        gTTCGCGAACAGGCGGCGGCATTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTggg                                                                                                       >  1:754501/1‑113 (MQ=255)                                                                            gcgAACAGGCGGCGGCCTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGc                                                                                                                    >  1:1264768/1‑96 (MQ=255)                                                                            gcgAACAGGCGGCGGCCTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGc                                                                                                                    <  2:1264768/96‑1 (MQ=255)                                                                                    gcggcggcCTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGGCCGCAATGCGCCGCAATGcgcg                                                      >  1:176506/1‑150 (MQ=255)                                                                                         ggcCTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTGCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGaaa                                                                                      >  2:973333/1‑113 (MQ=255)                                                                                         ggcCTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTGCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGaaa                                                                                      <  1:973333/113‑1 (MQ=255)                                                                                         ggcCTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGGCCGCAATGCGCCGCAATGCGCGGATgg                                                 >  1:791025/1‑150 (MQ=255)                                                                                            cTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGGCCGCAATGCGCCGCAATGCGCGGATGgcgc                                              >  2:28498/1‑150 (MQ=255)                                                                                             tttGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCt                                                                                                                   >  1:773368/1‑80 (MQ=255)                                                                                             tttGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCt                                                                                                                   <  2:773368/80‑1 (MQ=255)                                                                                               tGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGGCCGCAATGCGCCGCAATGCGCGGATGGCGCCGa                                           <  1:862516/150‑1 (MQ=255)                                                                                                    tGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGGCCGCAATGCGCCGCAATGCGCGGATGGCGCCGAGTTTc                                      >  1:847800/1‑150 (MQ=255)                                                                                                    tGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGGCCGCAATGCGCCGCAATGCGCGGATGGCGCCGAGTTTc                                      <  1:23897/150‑1 (MQ=255)                                                                                                      gTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGGCCGCAATGCGCCGCAATGCGCGGATGGCGCCGAGTTTCAg                                    <  2:1201128/150‑1 (MQ=255)                                                                                                      gTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGGCCGCAATGCGCCGCAATGCGCGGATGGCGCCGAGTTTCAg                                    >  2:319134/1‑150 (MQ=255)                                                                                                       tAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGGCCGCAATGCGCCGCAATGCGCGGATGGCGCCGAGTTTCAgg                                   <  1:843669/150‑1 (MQ=255)                                                                                                        ccagcgcAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGGCCGCAATGCGCCGCAATGCGCGGATGGCGCCGAGTTTCAGGc                                  <  1:1114568/144‑1 (MQ=255)                                                                                                          tCGCTAAAGTCGTGGTGCACACGAATCACCGGCAGTCCGGTCCCCAGCAGGTAAACGTAGCAGAGATTGACGCGGTGGGCGATGCAGAAAGCAAAATTGCGCTCGGACGCAATGCGCTGCAATGCGCGGATGGTGCCGAGATTCAGGCca                                >  2:605927/1‑150 (MQ=255)                                                                                                            gCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGGCCGCAATGCGCCGCAATGCGCGGATGGCGCCGAGTTTCAGGCcacg                              >  1:257804/1‑150 (MQ=255)                                                                                                             ccAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGGCCGCAATGCGCCGCAATGCGCGGATGGCGCCGATTTTCAGGCcacgc                             >  2:683344/1‑150 (MQ=255)                                                                                                             ccAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGGCCGCAATGCGCCGCAATGCGCGGATGGCGCCGAGTTTCAGGCcacgc                             >  2:20647/1‑150 (MQ=255)                                                                                                               aaaGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGGCCGCAATGCGCCGCAATGCGCGGATGGCGCCGAGTTTCAGGCcacgcac                           >  1:908043/1‑150 (MQ=255)                                                                                                                         tGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGGCCGCAATGCGCCGCAATGCGCGGATGGCGCCGAGTTTCAGGCCACGCACGGCCTTGgaa                 <  2:908043/150‑1 (MQ=255)                                                                                                                           cacaCCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGGCCGCAATGCGCCGCAATGCGCGGATGGCGCCGAGTTTCAGGCCACGCACGGCCTTGgaact               <  1:1311784/150‑1 (MQ=255)                                                                                                                            acacCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGGCCGCAATGCGCCGCAATGCGCGGATGGCGCCGAGTTTCAGGCcacgcacg                          <  2:339325/138‑1 (MQ=255)                                                                                                                            acacCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGGCCGCAATGCGCCGCAATGCGCGGATGGCGCCGAGTTTCAGGCcacgcacg                          >  1:339325/1‑138 (MQ=255)                                                                                                                            acacCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGGCCGCAATGCGCCGCAATGCGCGGATGGCGCCGAGTTTCAGGCCACGCACGGCCTTGgaactg              <  2:1309665/150‑1 (MQ=255)                                                                                                                            acacCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGGCCGCAATGCGCCGCAATGCGCGGATGGCGCCGAGTTTCAGGCCACGCACGGCCTTGgaactg              <  1:20647/150‑1 (MQ=255)                                                                                                                            acacCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGGCCGCAATGCGCCGCAATGCGCGGATGGCGCCGAGTTTCAGGCCACGCACGGCCTTGgaactg              >  1:123677/1‑150 (MQ=255)                                                                                                                             cacCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGgccgcaatgcgcc                                                                >  1:301375/1‑99 (MQ=255)                                                                                                                             cacCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGgccgcaatgcgcc                                                                <  2:301375/99‑1 (MQ=255)                                                                                                                              acCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCCGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGGCCGCAATGCGCCGCAATGCGCGGATGGCGCCGAGTTTCAGGCCACGCACGGCCTTGgaactgaa            >  2:457509/1‑150 (MQ=255)                                                                                                                               ccAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGGCCGCAATGCGCCGCAATGCGCGGATGGCGCCGAGTTTCAGGCCACGCACGGCCTTGgaactgaac           >  2:319688/1‑150 (MQ=255)                                                                                                                                cAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGGCCGCAATGCGCCGCAATGCGCGGATGGCGCCGAGTTTCAGGCCACGCACGGCCTTGgaactgaact          >  1:418935/1‑150 (MQ=255)                                                                                                                                   tCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGGCCGCAATGCGCCGCAATGCGCGGATGGCGCCGAGTTTCAGGCCACGCACGGCCTTGGAACTGAACTCaa       <  1:874814/150‑1 (MQ=255)                                                                                                                                   tCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGGCCGCAATGCGCCGCAATGCGCGGATGGCGCCGAGTTTCAGGCCACGCACGGCCTTGGAACTGAACTCaa       <  1:949160/150‑1 (MQ=255)                                                                                                                                    cACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGGCCGCAATGCGCCGCAATGCGCGGATGGCGCCGAGTTTCAGGCCACGCACGGCCTTGGAACTGAACTCAAg      <  2:477203/150‑1 (MQ=255)                                                                                                                                     aCCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGGCCGCAATGCGCCGCAATGCGCGGATGGCGCCGAGTTTCAGGCCACGCACGGCCtt                     >  2:969028/1‑134 (MQ=255)                                                                                                                                     aCCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGGCCGCAATGCGCCGCAATGCGCGGATGGCGCCGAGTTTCAGGCCACGCACGGCCtt                     <  1:969028/134‑1 (MQ=255)                                                                                                                                        ggCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGGCCGCAATGCGCCGCAATGCGCGGATGGCGCCGAGTTTCAGGCCACGCACGGCCTTGGAACTGAACTCAAGGAAc  <  2:176506/150‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                   GGCAAGCAACGGCGCATATCGTCACGTACCGCGTCCGAAACACCCAACAGGCTCAGACGCTTGCTGAACAGGTTCGCGAACAGGCGGCGGCCTTTGCGCTGGTAGTCGCCAAAGGCATGGTGCACACCAATCACCGGCAGGCCGGTCCCCAGCAGGGCAACGTAGATCGGCTTGAAGCGGTGGGCGATGCAGAAAGCGAAATTGCGCTCGGCCGCAATGCGCCGCAATGCGCGGATGGCGCCGAGTTTCAGGCCACGCACGGCCTTGGAACTGAACTCAAGGAAC  >  NC_002947/5625345‑5625629

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