Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,498,831	C→A	8.1%	K391N (AAG→AAT)	PP_3988 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,498,831	0	C	A	8.1%	187.9 / 5.8	62	K391N (AAG→AAT)	PP_3988	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (32/25);  new base A (3/2);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

GCTCAAATCTGGTTTGGCTCGACCCACAAAAACCAACGGCTTACTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCGAAAGGGTGCTTCTTAAGATAGGGC  >  NC_002947/4498684‑4498972                                                                                                                                                      |                                                                                                                                              gCTCAAATCTGGTTTGGCTCGACCCACAAAAACCAACGGCTTACTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCGAAAGAGTGGGACTCTCTGGACt                                                                                                                                              >  4:148718/1‑149 (MQ=255)   tCAAATCTGGTTTGGCTCGACCCACAAAAACCAACGGCTTACTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCAA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACtt                                                                                                                                             >  1:309507/1‑148 (MQ=255)     aaaTCTGGTTTGGCTCGACCCACAAAAACCAACGGCTTACTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAg                                                                                                                                          >  2:102382/1‑149 (MQ=255)      aaTCTGGTTTGGCTCGACCCACAAAAACCAACGGCTTACTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGc                                                                                                                                         <  3:62564/149‑1 (MQ=255)       aTCTGGTTTGGCTCGACCCACAAAAACCAACGGCTTACTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCt                                                                                                                                        <  6:23048/149‑1 (MQ=255)            gTTTGGCTCGACCCACAAAAACCAACGGCTTACTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTtata                                                                                                                                    <  8:221118/148‑1 (MQ=255)            gTTTGGCTCGACCCACAAAAACCAACGGCTTACTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTtata                                                                                                                                    <  1:202822/148‑1 (MQ=255)                 gCTCGACCCACAAAAACCAACGGCTTACTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGTGACTCTCTGGACTTTAGCTTATAGCGCt                                                                                                                               <  2:366050/148‑1 (MQ=255)                     gACCCACAAAAACCAACGGCTTACTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCATCCATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCtgtg                                                                                                                            <  1:102382/149‑1 (MQ=255)                        ccACAAAAACCAACGGCTTACTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATc                                                                                                                       >  3:324918/1‑149 (MQ=255)                          acaAAAACCAACGGCTTACTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAgcg                                                                                                                                 >  7:321637/1‑137 (MQ=255)                          acaAAAACCAACGGCTTACTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAgcg                                                                                                                                 <  8:321637/137‑1 (MQ=255)                           caAAAACCAACGGCTTACTCACGGCATGTTCTAGGGATTTGAACTTCGCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CTTACACATTAAATTCCAAAGAGTGGGACTCTGTGGACTTTAGCTTATAGCGCTGTGCGATCGTc                                                                                                                    <  5:237222/149‑1 (MQ=255)                                   aaCGGCTTACTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGcc                                                                                                            >  7:343755/1‑149 (MQ=255)                                            cTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGaa                                                                                                    <  6:218492/148‑1 (MQ=255)                                              cACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGca                                                                                                 >  1:224416/1‑149 (MQ=255)                                                          tAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTg                                                                                     >  5:136247/1‑149 (MQ=255)                                                          tAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGAGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTg                                                                                     >  6:299031/1‑149 (MQ=255)                                                           aGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTa                                                                                                                                           >  2:220527/1‑94 (MQ=255)                                                           aGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTa                                                                                                                                           <  1:220527/94‑1 (MQ=255)                                                            gggATTTGAACTTCCCGTAGCTCAACGGTCGATAGAGACGAACACCGTCCA‑‑CATACAAATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCTTCTGATCGCCCTCGAGAAGCACATCAACCGCTGcc                                                                                   >  1:161167/1‑149 (MQ=255)                                                            gggATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGcac                                                                                                >  2:275147/1‑136 (MQ=255)                                                            gggATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGcac                                                                                                >  2:280092/1‑136 (MQ=255)                                                            gggATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGcac                                                                                                <  1:280092/136‑1 (MQ=255)                                                            gggATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGcac                                                                                                <  1:275147/136‑1 (MQ=255)                                                                         cccGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGttt                                                                        <  4:324918/147‑1 (MQ=255)                                                                                   aaCGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAAttt                                                             <  2:309507/148‑1 (MQ=255)                                                                                         cGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAAt                                                      >  7:219475/1‑149 (MQ=255)                                                                                           aTAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATc                                                     >  4:297191/1‑148 (MQ=255)                                                                                             agagACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGATGTGCGATCGTATGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGa                                                   >  1:328483/1‑148 (MQ=255)                                                                                                gaCGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGATTGTTATCAATATGAAAAATCGAAGcc                                               >  6:14878/1‑149 (MQ=255)                                                                                                     acacCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGt                                           >  2:47773/1‑148 (MQ=255)                                                                                                      cacCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCAAATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCa                                         >  3:93179/1‑149 (MQ=255)                                                                                                         cGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCTTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTg                                      <  2:224416/149‑1 (MQ=255)                                                                                                         cGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTg                                      >  1:324719/1‑149 (MQ=255)                                                                                                          gTCCA‑‑CATACACATTAAATTCCAAAGAGTCGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTgg                                     <  3:40782/149‑1 (MQ=255)                                                                                                           tCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTgg                                     >  2:90749/1‑148 (MQ=255)                                                                                                           tCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATAgg                                     >  4:83319/1‑148 (MQ=255)                                                                                                            ccA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCa                                   >  7:112426/1‑149 (MQ=255)                                                                                                                  aTACACATTAAATTCCAAAGAGTGGGACTATCTGGACTTTAGCTTATAGCGCTGTGCGATAGTCTGATCGCCATCGAGAAGCACATCAACCGCTGACCCTCGCCCTTTGTTATCAATTTGAAAAATCGATGCCCAGTCATTGGCATCa                                >  8:127220/1‑148 (MQ=255)                                                                                                                   tACACATTAAATTCCAAAGAGTGGGACTCTCAGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTCTTATCAATTTGAAAAATCGATGCCCCGTCATTGCcaacaa                               >  7:65766/1‑148 (MQ=255)                                                                                                                      acaTTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATc                           >  1:23045/1‑149 (MQ=255)                                                                                                                      acaTTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATc                           <  8:112426/149‑1 (MQ=255)                                                                                                                        gctAAAGTCCAGAGAGTCCCACTCTTTGGAATTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCg                          >  8:287651/3‑148 (MQ=38)                                                                                                                        gctAAAGTCCAGAGAGTCCCACTCTTTGGAATTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCg                          >  1:124702/3‑148 (MQ=38)                                                                                                                        aTTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCg                          >  7:115408/1‑148 (MQ=255)                                                                                                                        aTTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCg                          >  1:242082/1‑148 (MQ=255)                                                                                                                            aaTTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCGAAAg                      <  6:136247/148‑1 (MQ=255)                                                                                                                               tCCAGAGAGTCCCACTCTTTGGAATTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCGAAAGGGTg                  >  1:192606/1‑149 (MQ=255)                                                                                                                                 cAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCGAAAGGGTGct                >  3:248689/1‑149 (MQ=255)                                                                                                                                  aaaGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCGAAAGGGTGctt               >  2:118443/1‑149 (MQ=255)                                                                                                                                     gagTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCATATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCGAAAGGGTGcttctt            <  4:93179/149‑1 (MQ=255)                                                                                                                                      agTTGGACTCTCTGGACTTTATCATATAGCGCTGTGCGATCGTCTGATCTCCCTCGAGAATCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCTTCATTGGCATCAATATCGAAAGGGTGCTTCTTa           <  7:127220/149‑1 (MQ=255)                                                                                                                                       gTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCGAAAGGGTGCTTCTTaa          <  4:248689/149‑1 (MQ=255)                                                                                                                                       gTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCGAAAGGGTGCTTCTTaa          <  8:343755/149‑1 (MQ=255)                                                                                                                                        tcccaCTCTTTGGAATTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCGAAAGGGTGCTTCTTAAg         <  7:287651/145‑1 (MQ=255)                                                                                                                                          ccaCTCTTTGGAATTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCGAAAGGGTGCTTCTTAAGAt       <  2:192606/147‑1 (MQ=255)                                                                                                                                            aCTCTCTGGACTTTTGGTTATAGCGCTGTGCGATCGTCTGATCGCCGTCGAGAAGCACATGTACCGCTGCCGCTCGGCGTTTGTTATCATTTTGAATAATCGATGCCCCGTCATTGGGATCAATATCGAAAGGGTGCTTCTTAAGATAg     <  8:65766/149‑1 (MQ=255)                                                                                                                                            aCTCTCTGGACTTTAGCTTATAGCGCTGTGCTATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCGAAAGGGTGCTTCTTAAGATAg     <  8:219475/149‑1 (MQ=255)                                                                                                                                             ctctctGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCACTCGACGTTTTTTATAAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCGAAAGGGTGCTTCTTAAGATAgt    <  2:161167/149‑2 (MQ=255)                                                                                                                                                tctGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCGAAAGGGTGCTTCTTAAGATAGGGc  >  2:276311/1‑148 (MQ=255)                                                                                                                                                tctGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCGAAAGGGTGCTTCTTAAGATAGGGc  >  7:157433/1‑148 (MQ=255)                                                                                                                                                      |                                                                                                                                              GCTCAAATCTGGTTTGGCTCGACCCACAAAAACCAACGGCTTACTCACGGCATGATCTAGGGATTTGAACTTCCCGTAGCTCAACGGGCGATAGAGACGAACACCGTCCA‑‑CATACACATTAAATTCCAAAGAGTGGGACTCTCTGGACTTTAGCTTATAGCGCTGTGCGATCGTCTGATCGCCCTCGAGAAGCACATCAACCGCTGCCCCTCGCCGTTTGTTATCAATTTGAAAAATCGATGCCCCGTCATTGGCATCAATATCGAAAGGGTGCTTCTTAAGATAGGGC  >  NC_002947/4498684‑4498972

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: ‑