Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,485,238	G→T	6.0%	A243A (GCC→GCA)	PP_3091 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,485,238	0	G	T	6.0%	210.6 / 2.1	66	A243A (GCC→GCA)	PP_3091	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (30/32);  new base T (2/2);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

TGCAGGCCGTTGGTGCGGATTTCCGAAAGGCCAAAGGTCATGCCCAGCGGCTGCTGGCGTTCGTGCAGGTCGAGGCCGTCGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGCGCTCATCGCTGTTGGCCGCCAGCAGTTCGGGCAGGCTGACCACCAGGACCAGGCCATTGATCGCCGCTTTGCCGCGCTGCTTA  >  NC_002947/3485095‑3485382                                                                                                                                                |                                                                                                                                                 tGCAGGCCGTTGGTGCGGATTTCCGAAAGGCCAAAGGTCATGCCCAGCGGCTGCTGGCGTTCGTGCAGGTCGAGGCCGTCGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTTGCATCg                                                                                                                                            >  1:737933/1‑150 (MQ=255) tGCAGGCCGTTGGTGCGGATTTCCGAAAGGCCAAAGGTCATGCCCAGCGGCTGCTGGCGTTCGTGCAGGTCGAGGCCGTAGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTTGCATCg                                                                                                                                            >  1:737888/1‑150 (MQ=255)   cAGGCCGTTGGTGCGGATTTCCGAAAGGCCAAAGGTCATGCCCAGCGGCTGCTGGCGTTCGTGCAGGTCGAGGCCGTCGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTTGCATCGAg                                                                                                                                          <  2:737933/150‑1 (MQ=255)   cAGGCCGTTGGTGCGGATTTCCGAAAGGCCAAAGGTCATGCCCAGCGGCTGCTGGCGTTCGTGCAGGTCGAGGCCGTCGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTTGCATCGAg                                                                                                                                          <  2:737888/150‑1 (MQ=255)   cAGGCCGTTGGTGCGGATTTCCGAAAGGCCAAAGGTCATGCCCAGCGGCTGCTGGCGTTCGTGCAGGTCGAGGCCGTCGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAg                                                                                                                                          <  1:703161/150‑1 (MQ=255)      gCCGTTGGTGCGGATTTCCGAAAGGCCAAAGGTCATGCCCAGCGGCTGCTGGCGTTCGTGCAGGTCGAGGCCGTCGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCa                                                                                                                                       >  1:613845/1‑150 (MQ=255)                  gATTTCCGAAAGGCCAAAGGTCATGCCCAGCGGCTGCTGGCGTTCGTGCAGGTCGAGGCCGTCGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTAttc                                                                                                                           <  1:1145558/150‑1 (MQ=255)                     ttCCGAAAGGCCAAAGGTCATGCCCAGCGGCTGCTGGCGTTCGTGCAGGTCGAGGCCGTCGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTAttcttc                                                                                                                        <  1:1003340/150‑1 (MQ=255)                         gAAAGGCCAAAGGTCATGCCCAGCGGCTGCTGGCGTTCGTGCAGGTCGAGGCCGTCGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCAcc                                                                                                                    >  2:1210462/1‑150 (MQ=255)                         gAAAGGCCAAAGGTCATGCCCAGCGGCTGCTGGCGTTCGTGCAGGTCGAGGCCGTCGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCAcc                                                                                                                    >  2:111592/1‑150 (MQ=255)                         gAAAGGCCAAAGGTCATGCCCAGCGGCTGCTGGCGTTCGTGCAGGTCGAGGCCGTCGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCAcc                                                                                                                    >  1:1262970/1‑150 (MQ=255)                         gAAAGGCCAAAGGTCATGCCCAGCGGCTGCTGGCGTTCGTGCAGGTCGAGGCCGTCGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCAcc                                                                                                                    >  2:889177/1‑150 (MQ=255)                         gAAAGGCCAAAGGTCATGCCCAGCGGCTGCTGGCGTTCGTGCAGGTCGAGGCCGTCGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATGCTTCCAcc                                                                                                                    >  2:310231/1‑150 (MQ=255)                          aaaGGCCAAAGGTCATGCCCAGCGGCTGCTGGCGTTCGTGCAGGTCGAGGCCGTCGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCAccc                                                                                                                   >  1:691146/1‑150 (MQ=255)                            aGGCCAAAGGTCATGCCCAGCGGCTGCTGGCGTTCGTGCAGGTCGAGGCCGTCGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAATCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCgg                                                                                                                 <  2:984888/150‑1 (MQ=255)                            aGGCCAAAGGTCATGCCCAGCGGCTGCTGGCGTTCGTGCAGGTCGAGGCCGTCGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCCGCGTATTCTTCCACCCgg                                                                                                                 <  1:600347/150‑1 (MQ=255)                                cAAAGGTCATGCCCAGCGGCTGCTGGCGTTCGTGCAGGTCGAGGCCGTCGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGcca                                                                                                             >  2:1209574/1‑150 (MQ=255)                                 aaaGGTCATGCCCAGCGGCTGCTGGCGTTCGTGCAGGTCGAGGCCGTCGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGccac                                                                                                            <  1:1012222/150‑1 (MQ=255)                                 aaaGGTCATGCCCAGCGGCTGCTGGCGTTCGTGCAGGTCGAGGCCGTCGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGccac                                                                                                            <  2:691146/150‑1 (MQ=255)                                   aGGTCATGCCCAGCGGCTGCTGGCGTTCGTGCAGGTCGAGGCCGTCGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGccacca                                                                                                          >  2:41537/1‑150 (MQ=255)                                     gTCATGCCCAGCGGCTGCTGGCGTTCGTGCAGGTCGAGGCCAACGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGAATCAGGTAGACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGc                                                                                                        >  2:339041/1‑150 (MQ=255)                                      tCATGCCCAGCGGCTGCTGGCGTTCGTGCAGGTCGAGGCCGTCGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCt                                                                                                       <  1:321882/150‑1 (MQ=255)                                        aTGCCCAGCGGCTGCTGGCGTTCGTGCAGGTCGAGGCCGTCGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTgg                                                                                                     <  2:1053974/150‑1 (MQ=255)                                           cccAGCGGCTGCTGGCGTTCGTGCAGGTCGAGGCCGTCGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCa                                                                                                  <  1:1308057/150‑1 (MQ=255)                                               gCGGCTGCTGGCGTTCGTGCAGGTCGAGGCCGTCGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCa                                                                                              <  1:975286/150‑1 (MQ=255)                                                                tGCAGGTCGAGGCCGTCGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAg                                                                                             <  1:547822/134‑1 (MQ=255)                                                                tGCAGGTCGAGGCCGTCGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAg                                                                                             >  2:547822/1‑134 (MQ=255)                                                                 gCAGGTCGAGGCCGTCGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGCGCTCATCGc                                                                            >  1:1236119/1‑150 (MQ=255)                                                                   aGGTCGAGGCCGTCGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCgg                                                                                                                 <  1:402790/111‑1 (MQ=255)                                                                   aGGTCGAGGCCGTCGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCgg                                                                                                                 >  2:402790/1‑111 (MQ=255)                                                                             cGTCGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCCTCCGGTAAACCCGCGGGTTGGCATCGACGCCTACGGCGTACTCTACCACCAGGCCCATCATCTCCGCAGGCAGGGCATTGCGCTCATCGCTGTTGGCCGCca                                                                <  1:1097118/150‑1 (MQ=255)                                                                                cGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGCGCTCATCGCTGTTGGCCGCca                                                                <  1:4890/147‑1 (MQ=255)                                                                                cGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGCGCTCATCGCTGTTGGCCGCca                                                                >  2:4890/1‑147 (MQ=255)                                                                                 gAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAg                                                                                                                                          <  2:872345/72‑1 (MQ=255)                                                                                 gAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAg                                                                                                                                          >  1:872345/1‑72 (MQ=255)                                                                                     gCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGCGCTCATCGCTGTTGGCCGCCAGCAGTTCg                                                        <  1:889177/150‑1 (MQ=255)                                                                                      ccTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGCGCTCATCGCTGTTGGCCGCCAGCAGTTCgg                                                       >  2:109515/1‑150 (MQ=255)                                                                                        tGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGCGCTCATCGCTGTTGGCCGCCAGCAGTTCGGGc                                                     >  2:174853/1‑150 (MQ=255)                                                                                          gCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGCGCTCATCGCTGTTGGCCGCCAGCAGTTCGGGCAg                                                   >  2:1395810/1‑150 (MQ=255)                                                                                            tGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGCGCTCATCGCTGTTGGCCGCCAGCAGTTCGGGCAgg                                                  <  2:21586/149‑1 (MQ=255)                                                                                            tGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGCGCTCATCGCTGTTGGCCGCCAGCAGTTCGGGCAgg                                                  >  1:21586/1‑149 (MQ=255)                                                                                                gCCCGGCAACCGGACCGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATGGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGTGCTCATCGCTGTTGGCCGCCAGCAGTTCGGGCAGGCTGac                                             <  1:1359968/150‑1 (MQ=255)                                                                                                 cccGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGCGCTCATCGATGTTGGCCTCCAGCAGTTCGTGCAGGCGGagc                                            >  2:1016461/1‑148 (MQ=255)                                                                                                  ccGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGCGCTCATCGCTGTTGGCCGCCAGCAGTTCGGGCAGGCTGacca                                           >  1:1419452/1‑150 (MQ=255)                                                                                                   cGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGCGCTCATCGCTGTTGGCCGCCAGCAGTTCgg                                                       <  2:644174/137‑1 (MQ=255)                                                                                                   cGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGCGCTCATCGCTGTTGGCCGCCAGCAGTTCgg                                                       >  1:644174/1‑137 (MQ=255)                                                                                                    ggCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGCGCTCATCGCTGTTGGCCGCCAGCAGTTCGGGCAGGCTGaccacc                                         <  1:109515/150‑1 (MQ=255)                                                                                                       aaCTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGCGCTCATCGCTGTTGGCCGCCAGCAGTTCGGGCAGGCTGACCaccagg                                      <  1:1209574/150‑1 (MQ=255)                                                                                                         cTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGCGCTCATCGCTGTTGGCCGCCAGCAGTTCGGGCAGGCTGACCaccaggac                                    >  1:254830/1‑150 (MQ=255)                                                                                                            gtcggtcTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAg                                                                                             >  1:294613/1‑90 (MQ=255)                                                                                                            gtcggtcTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAg                                                                                             <  2:294613/90‑1 (MQ=255)                                                                                                            gtcggtcTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGCGCTCATCGCTGTTGGCAGCCAGCAGTTCGGGCAGGCTGACCaccaggaccag                                 <  1:1101976/150‑1 (MQ=255)                                                                                                             tcggtcTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGCGCTCATCGCTGTTGGCCGCCAGCAGTTCGGGCAGGCTGACCaccaggaccagg                                <  2:613845/150‑1 (MQ=255)                                                                                                             tcggtcTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGCGCTCATCGCTGTTGGCCGCCAGCAGTTCGGGCAGGCTGACCaccaggaccagg                                <  2:1419452/150‑1 (MQ=255)                                                                                                               ggtcTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGCGCTCATCGCTGTTGGCCGCCAGCAGTTCGGGCAGGCTGACCACCAGGACCAGGcc                              <  1:111592/150‑1 (MQ=255)                                                                                                               ggtcTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGCGCTCATCGCTGTTGGCCGCCAGCAGTTCGGGCAGGCTGACCACCAGGACCAGGcc                              >  1:3180/1‑150 (MQ=255)                                                                                                                       tCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGCGCTCATCGCTGTTGGCCGCCAGCAGTTCGGGCAGGCTGACCACCAGGACCAGGCCATTGATcg                      >  1:351170/1‑150 (MQ=255)                                                                                                                             tCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGCGCTCATCGCTGTTGGCCGCCAGCAGTTc                                                         >  2:751623/1‑109 (MQ=255)                                                                                                                             tCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGCGCTCATCGCTGTTGGCCGCCAGCAGTTc                                                         <  1:751623/109‑1 (MQ=255)                                                                                                                             gcatttAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGCGCTCATCGCTGTTGGCCGCCAGCAGTTCGGGCAGGCTGACCACCAGGACCAGGCCATTGATCGCCGCtt                >  2:1423220/6‑150 (MQ=255)                                                                                                                                      ccGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGCGCTCATCGCTGTTGGCCGCCAGCAGTTCGGGCAGGCTGACCACCAGGACCAGGCCATTGATCGCCGCTTTGCCGCgct       <  1:1395810/150‑1 (MQ=255)                                                                                                                                        ggcggGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGCGCTCATCGCTGTTGGCCGCCAGCAGTTCGGGCAGGCTGACCACCAGGACCAGGCCATTGATCGCCGCTTTGCCGCgctgc     <  2:3180/150‑1 (MQ=255)                                                                                                                                         gcggGTTGGCATCGAGGCATTCGGCGTCTTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGCGCTCATCGCTGTTGGCCGCCAGCAGTTCGGGCAGGCTGACCACCAGGACCAGGCCATTGATCGCCGCTTTGCCGCgctgct    >  2:471051/1‑150 (MQ=255)                                                                                                                                           gggTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGCGCTCATCGCTGTTGGCCGCca                                                                >  1:817106/1‑88 (MQ=255)                                                                                                                                           gggTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGCGCTCATCGCTGTTGGCCGCca                                                                <  2:817106/88‑1 (MQ=255)                                                                                                                                           gggTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGCGCTCATCGCTGTTGGCCGCCAGCAGTTCGGGCAGGCTGACCACCAGGACCAGGCCATTGATCGCCGCTTTGCCGCGCTGCTTa  <  1:1016461/150‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                 TGCAGGCCGTTGGTGCGGATTTCCGAAAGGCCAAAGGTCATGCCCAGCGGCTGCTGGCGTTCGTGCAGGTCGAGGCCGTCGAAAGCCTGGCTGAAGCCCGGCAACTGGTCGGTCTTGCTCAGCATCAGGTAAACCGGCGGGTTGGCATCGAGGCATTCGGCGTATTCTTCCACCCGGGCCACCAGCTGGGCAGCCAGGGCATTGCGCTCATCGCTGTTGGCCGCCAGCAGTTCGGGCAGGCTGACCACCAGGACCAGGCCATTGATCGCCGCTTTGCCGCGCTGCTTA  >  NC_002947/3485095‑3485382

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