Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,921,633	C→G	5.9%	A666G (GCC→GGC)	PP_4330 →	toxin‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,921,633	0	C	G	5.9%	219.0 / 5.0	69	A666G (GCC→GGC)	PP_4330	toxin‑like protein
Reads supporting (aligned to +/- strand):  ref base C (28/36);  new base G (2/2);  total (30/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.25e-01

AAAGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGC  >  NC_002947/4921491‑4921778                                                                                                                                                |                                                                                                                                                  aaaGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAggg                                                                                                                                            <  1:1629886/150‑1 (MQ=255) aaaGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAgga                                                                                                                                            <  2:1332377/150‑2 (MQ=255)   aGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCt                                                                                                                                          >  2:200674/1‑150 (MQ=255)       aTTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAg                                                                                                                                      <  1:1838356/150‑1 (MQ=255)                 cATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAg                                                                                                                                      <  1:1005821/140‑1 (MQ=255)                 cATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAg                                                                                                                                      >  2:1005821/1‑140 (MQ=255)                 cATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGc                                                                                                                            <  1:830114/150‑1 (MQ=255)                 cATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGc                                                                                                                            >  2:928128/1‑150 (MQ=255)                   ttGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCtt                                                                                                                          <  2:1433484/150‑1 (MQ=255)                      cTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTGTTGGTAGCCA‑GCGCCGCAGGGCTGCAGGTATTGCCGCTTGaa                                                                                                                       <  2:21692/150‑1 (MQ=255)                         gCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAagag                                                                                                                    >  1:320949/1‑150 (MQ=255)                         gCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAagag                                                                                                                    <  2:1671866/150‑1 (MQ=255)                               cTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGaa                                                                                                                       <  2:1587531/141‑1 (MQ=255)                               cTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGaa                                                                                                                       >  1:1587531/1‑141 (MQ=255)                                    cgccAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGcttctt                                                                                                         >  1:936671/1‑150 (MQ=255)                                       cAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAg                                                                                                                                      <  2:648887/118‑1 (MQ=255)                                       cAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAg                                                                                                                                      >  1:648887/1‑118 (MQ=255)                                        aaCGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCACGGCTTCTTCCgg                                                                                                     >  2:1860644/1‑150 (MQ=255)                                        aaCGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCgg                                                                                                     >  2:1093054/1‑150 (MQ=255)                                               cgcgcCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCt                                                                                              <  1:1489144/150‑1 (MQ=255)                                               cgcgcCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCt                                                                                              <  2:483528/150‑1 (MQ=255)                                                     tGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCaga                                                                                        <  2:1151490/150‑1 (MQ=255)                                                           aCAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCt                                                                                  >  2:58471/1‑150 (MQ=255)                                                            cAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCa                                                                                          >  2:1648777/1‑141 (MQ=255)                                                            cAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCa                                                                                          <  1:1648777/141‑1 (MQ=255)                                                               tGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAg                                                                                                                >  2:1239630/1‑116 (MQ=255)                                                               tGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAg                                                                                                                <  1:1239630/116‑1 (MQ=255)                                                                      aGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGAAAGTCTTCTCCGGGGTCTTCTTCCAGATCCGCTATTACTCGCAg                                                                       >  2:610822/1‑150 (MQ=255)                                                                          tGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCggg                                                                                                    <  2:486188/117‑1 (MQ=255)                                                                          tGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCggg                                                                                                    >  1:486188/1‑117 (MQ=255)                                                                           gcgcCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGa                                                                  <  1:848163/150‑1 (MQ=255)                                                                                 gTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCaga                                                                                        >  1:1379166/1‑122 (MQ=255)                                                                                 gTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCaga                                                                                        <  2:1379166/122‑1 (MQ=255)                                                                                    tCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCt                                                         <  1:200674/150‑1 (MQ=255)                                                                                       gTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGa                                                      <  1:58471/150‑1 (MQ=255)                                                                                       gTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGG‑CCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGAt                                                     <  2:1584301/150‑1 (MQ=255)                                                                                        tCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCaa                                                                                                                 >  2:454749/1‑90 (MQ=255)                                                                                        tCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCaa                                                                                                                 <  1:454749/90‑1 (MQ=255)                                                                                        tCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGAt                                                     <  2:320949/150‑1 (MQ=255)                                                                                          gAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTAcgc                                                                           >  2:402405/1‑126 (MQ=255)                                                                                          gAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTAcgc                                                                           <  1:402405/126‑1 (MQ=255)                                                                                           aaGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGc                                                                                                                               <  2:1831716/73‑1 (MQ=255)                                                                                           aaGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGc                                                                                                                               >  1:1831716/1‑73 (MQ=255)                                                                                              gCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAAc                                               <  2:936671/150‑1 (MQ=255)                                                                                                     cTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCa                                        >  2:352503/1‑150 (MQ=255)                                                                                                           cTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTAt                                                                                                                                  <  2:159221/54‑1 (MQ=255)                                                                                                           cTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTAt                                                                                                                                  >  1:159221/1‑54 (MQ=255)                                                                                                            ttACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGccac                                 >  2:1778073/1‑150 (MQ=255)                                                                                                                gAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGt                                            >  2:1852312/1‑135 (MQ=255)                                                                                                                gAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGt                                            <  1:1852312/135‑1 (MQ=255)                                                                                                                     gAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTAcgc                        <  1:1778073/150‑1 (MQ=255)                                                                                                                     gAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTAcgc                        <  1:928128/150‑1 (MQ=255)                                                                                                                       aagcggcaatacctGCAGCCCTGGGGCCC‑GGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATgg                                                   <  2:1330399‑M1/107‑1 (MQ=255)                                                                                                                       aagcggcaatacctGCAGCCC‑TGGGGCCCGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATgg                                                   >  1:1330399‑M1/15‑121 (MQ=255)                                                                                                                           aaGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGCCGGCAGCGCCACCAATTACGCGTCttc                  >  2:932943/1‑148 (MQ=255)                                                                                                                             ggCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGt                                            <  2:582320/122‑1 (MQ=255)                                                                                                                             ggCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGt                                            >  1:582320/1‑122 (MQ=255)                                                                                                                                tATTGGCAGCCT‑GGGCCACAGGACTGCAGGGATTCACGCTTGAAGAGGAAGGCTTCTTCCGGGTCTTCAGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCGGGATGGCAACGGTGGCATCGCCACTAATTACGCGTCTGCGACCg             >  2:437367/1‑150 (MQ=255)                                                                                                                                tATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGt                                            >  1:203789/1‑119 (MQ=255)                                                                                                                                tATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGt                                            <  2:203789/119‑1 (MQ=255)                                                                                                                                 aTTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGa                                                                  >  1:632774/1‑96 (MQ=255)                                                                                                                                 aTTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGa                                                                  <  2:632774/96‑1 (MQ=255)                                                                                                                                    tGCAGCCCTGGGGCCC‑GGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGc                                                                                      >  2:1450864‑M1/2‑73 (MQ=255)                                                                                                                                    tGCAGCCCTGGGGCCC‑GGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGc                                                                                      <  1:1450864‑M1/72‑1 (MQ=255)                                                                                                                                     GCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCGAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGAtgtg        <  2:1103772‑M1/150‑1 (MQ=255)                                                                                                                                        GCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGtgctg     >  1:1657749‑M1/1‑150 (MQ=255)                                                                                                                                        GCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGtgctg     <  1:249828‑M1/150‑1 (MQ=255)                                                                                                                                          CC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTgcg   <  2:1409679‑M1/150‑1 (MQ=255)                                                                                                                                           C‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTgcgc  <  1:8340‑M1/150‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                  AAAGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCC‑GGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGC  >  NC_002947/4921491‑4921778

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