Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,817,899	G→T	6.9%	G199C (GGT→TGT)	PP_4242 →	conserved protein of unknown function with Rhs repeat‑associated core domain

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,817,899	0	G	T	6.9%	206.5 / 4.0	70	G199C (GGT→TGT)	PP_4242	conserved protein of unknown function with Rhs repeat‑associated core domain
Reads supporting (aligned to +/- strand):  ref base G (30/35);  new base T (3/2);  total (33/37)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GGTGGCATCAGGAGCTCTACGTGGTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGAAAGTTATTTAATTGAAGAACATGCATTTACTAATATTCGAGTAC  >  NC_002947/4817752‑4818046                                                                                                                                                    |                                                                                                                                                    ggtggCATCAGGAGCTCTACGTGGTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATTGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAg                                                                                                                                                     <  7:151222/148‑1 (MQ=255)   tggCCTCAGGAGCTCTACGTGGTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGTCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTg                                                                                                                                                  <  2:103589/149‑1 (MQ=255)   tggCATCAGGAGCTCTACGTGGTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGt                                                                                                                                                   >  4:275245/1‑148 (MQ=255)   tggCATCAGGAGCTCTACGTGGTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGt                                                                                                                                                   >  2:199079/1‑148 (MQ=255)       aTCAGGAGCTCTACGTGGTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGc                                                                                                                                              >  4:270125/1‑149 (MQ=255)          agtaaatgcatgttcttcaattaaataactttccctggttgcctcatttttcccatcccatatttttttgatcttgcccgcagccacaagcacatcgaaatccttgatgtacacaaccaaatacaCGAGCTTCACAAAGGTGAAGCTc                                                                                                                                            >  1:238524‑M1/126‑148 (MQ=255)          aGGAGCTCTACGTGGTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTc                                                                                                                                            <  5:272431/148‑1 (MQ=255)          aGGAGCTCTACGTGGTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCg                                                                                                                                           <  8:297683/149‑1 (MQ=255)          aGGAGCTCTACGTGGTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCg                                                                                                                                           >  3:268578/1‑149 (MQ=255)            gtaaatgcatgttcttcaattaaataactttccctggttgcctcatttttcccagcccatattattttgatcttgcccgcagccacaagcacatcgaaatccttgatgtacacaccaaatacaCGAGCTTCACAAAGGTGAAGCTCgt                                                                                                                                          >  5:187421‑M1/124‑148 (MQ=255)                   aCGTGGTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTgg                                                                                                                                  >  4:49941/1‑149 (MQ=255)                    cGTGGTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTgg                                                                                                                                  <  3:270125/148‑1 (MQ=255)                    cGTGGTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGgt                                                                                                                                 <  6:59582/149‑1 (MQ=255)                    cGTGGTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGgt                                                                                                                                 <  3:275245/149‑1 (MQ=255)                     gTGGTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGgtg                                                                                                                                <  2:288882/149‑1 (MQ=255)                      tGGTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGgtgt                                                                                                                               >  7:26464/1‑149 (MQ=255)                       ggTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGgtgtg                                                                                                                              >  6:110190/1‑149 (MQ=255)                        gTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGgtgtgt                                                                                                                             <  6:99959/149‑1 (MQ=255)                        gTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGgtgtg                                                                                                                              <  7:290736/148‑1 (MQ=255)                        gTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGgtgtg                                                                                                                              >  1:68520/1‑148 (MQ=255)                             tAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTAc                                                                                                                           <  1:50890/146‑1 (MQ=255)                              agagTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATc                                                                                                                        <  6:79495/148‑1 (MQ=255)                                agTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAg                                                                                                                     <  8:320677/149‑1 (MQ=255)                                  tCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATAGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGa                                                                                                                   >  3:276246/1‑149 (MQ=255)                                    tGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGAtt                                                                                                                 >  5:31408/1‑149 (MQ=255)                                    tGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGAtt                                                                                                                 >  7:338572/1‑149 (MQ=255)                                        gTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGa                                                                                                             >  6:382648/1‑149 (MQ=255)                                         ttCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGAt                                                                                                            <  5:110190/149‑1 (MQ=255)                                         ttCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGAt                                                                                                            <  4:268578/149‑1 (MQ=255)                                         ttCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGAt                                                                                                            <  6:31408/149‑1 (MQ=255)                                           cAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGGTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGAtgt                                                                                                          <  8:185039/149‑1 (MQ=255)                                                      aCATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTGGATGTGCTTGGGGCTg                                                                                               >  6:282208/1‑149 (MQ=255)                                                      aCATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTg                                                                                               <  1:199079/149‑1 (MQ=255)                                                           aaaTTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCggg                                                                                           >  3:34003/1‑148 (MQ=255)                                                            aaTTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGc                                                                                          >  2:8148/1‑148 (MQ=255)                                                            aaTTGTAGAACAAGGGCACGGCTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAATGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGc                                                                                          >  4:301665/1‑148 (MQ=255)                                                                gttgATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGa                                                                                      <  3:130369/148‑1 (MQ=255)                                                                gttgATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGa                                                                                      <  3:49941/148‑1 (MQ=255)                                                                     tCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCaaa                                                                                 >  8:124946/1‑148 (MQ=255)                                                                        tGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCaaaaaaa                                                                             <  1:22353/149‑1 (MQ=255)                                                                            cAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATatg                                                                         <  8:26464/149‑1 (MQ=255)                                                                              ggggTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATatggg                                                                       >  7:294740/1‑149 (MQ=255)                                                                                ggTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATatgggat                                                                     >  1:45089/1‑149 (MQ=255)                                                                                     tattttaCACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATAtgggatggga                                                                 >  2:81854/1‑148 (MQ=255)                                                                                             acGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAgg                                                        >  1:371009/1‑149 (MQ=255)                                                                                              cGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGc                                                       <  1:8148/149‑1 (MQ=255)                                                                                               cacaagcacatcgaaatccttgatgtacacaccaaatacaCGAGCTTCACCTTTGTGAAGCTTGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGc                                                       <  4:289217‑M1/108‑1 (MQ=255)                                                                                                 aaCTATCGCGTGGCCGAGCCCTATGCGAATGCAGCCGCCAAGCATCACAAAGCTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCGTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAAc                                                    <  6:190951/149‑1 (MQ=255)                                                                                                    tATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAg                                                 >  4:209079/1‑149 (MQ=255)                                                                                                      tGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAATAATGAGGCAACCAggg                                               >  4:364539/1‑149 (MQ=255)                                                                                                        gCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGa                                              >  8:182645/1‑148 (MQ=255)                                                                                                           tGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGAAAGt                                          <  4:34003/149‑1 (MQ=255)                                                                                                               cGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGAAAGTTAtt                                      <  3:209079/149‑1 (MQ=255)                                                                                                                gAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGAAAGTTAttt                                     <  7:124946/149‑1 (MQ=255)                                                                                                                 cttgatgtacacaccaaatacaCGAGCTTCACCTTTGTGAAGCTTGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGAAAGTTAttt                                     <  6:214684‑M1/126‑1 (MQ=255)                                                                                                                  ttgatgtacacaccaaatacaCGAGCTTCACCTTTGTGAAGCTTGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGAAAGTTATTTaa                                   >  5:214684‑M1/22‑149 (MQ=255)                                                                                                                     cTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGAAAGTTATTTAAtt                                 >  7:129103/1‑148 (MQ=255)                                                                                                                     cTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGAAAGTTATTTAATTg                                <  8:338572/149‑1 (MQ=255)                                                                                                                      tgtacacaccaaatacaCGAGCTTCACCTTTGTGAAGCTTGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGAAAGTTATTTAATTga                               >  3:289217‑M1/18‑149 (MQ=255)                                                                                                                        tGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGAAAGTTATTTAATTgaa                              <  7:182645/148‑1 (MQ=255)                                                                                                                           gAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCCTGTATTTGGTGTGTAAATAAAGGATTTCGATGTGCTTGTGGCTGCCGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGAAAGTTATTTAATTgaagaa                           >  2:313641/1‑148 (MQ=255)                                                                                                                                 aGCCGCCAAGCTTCACAAAGGTGAAGCTCGAGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGAAAGTTATTTAATTGAAGAACATGCAt                    >  1:30100/1‑149 (MQ=255)                                                                                                                                       CAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGAAAGTTATTTAATTGAAGAACATGCATTTACt               <  4:52421‑M1/148‑1 (MQ=255)                                                                                                                                       CAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGAAAGTTATTTAATTGAAGAACATGCATTTACTa              <  7:319914‑M1/149‑1 (MQ=255)                                                                                                                                         accatCACAACGGTCAAGCTCGTTTATTTGGGGTGTAAATCAACGATATCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGAAAGTTATTTAATTGAAGAACATGCATTTACTAat            <  4:276246/145‑1 (MQ=255)                                                                                                                                         aGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGAAAGTTATTTAATTGAAGAACATGCATTTACTAat            <  2:371009/149‑1 (MQ=255)                                                                                                                                              cacaAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGAAAGTTATTTAATTGAAGAACATGCATTTACTAATATTc        <  8:294740/148‑1 (MQ=255)                                                                                                                                               acaAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGAAAGTTATTTAATTGAAGAACATGCATTTACTAATATTCGa      <  8:129103/149‑1 (MQ=255)                                                                                                                                                 aaaGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGAAAGTTATTTAATTGAAGAACATGCATTTACTAATATTCGAg     >  8:93468/1‑148 (MQ=255)                                                                                                                                                 aaaGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATCGGATGGGAAAAATGAGGCAACCAGGGAAAGTTATTTAATTGAAGAACATGCATTTACTAATATTCGAGt    <  2:30100/149‑1 (MQ=255)                                                                                                                                                   aGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGAAAGTTAGTTAATTGAAGAACATGCATTTCCTAATATTCGAGTac  >  5:269032/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GGTGGCATCAGGAGCTCTACGTGGTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGAAAGTTATTTAATTGAAGAACATGCATTTACTAATATTCGAGTAC  >  NC_002947/4817752‑4818046

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