Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,282,330	C→T	8.1%	E27K (GAG→AAG)	PP_2883 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,282,330	0	C	T	8.1%	168.2 / 4.7	62	E27K (GAG→AAG)	PP_2883	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (26/31);  new base T (2/3);  total (28/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 = 7.95e-01

CCAGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGAAAC  >  NC_002947/3282187‑3282475                                                                                                                                                |                                                                                                                                                   ccAGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAggtg                                                                                                                                              <  1:1069690/150‑1 (MQ=255)   aGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAggtgg‑t                                                                                                                                           >  1:868102/1‑150 (MQ=255)    gtgtTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCA‑CCACCTTGAggtggctg                                                                                                                                          <  1:829935/150‑4 (MQ=255)       ttCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCAc                                                                                                                                       >  2:682189/1‑150 (MQ=255)       ttCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCAc                                                                                                                                       <  2:396483/150‑1 (MQ=255)           aGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCA‑CCACCTTGAGGTGGCTGCACGc                                                                                                                                     <  1:1058945/148‑1 (MQ=255)           aGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCATTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCA‑CCACCTTGAGGTGGCTGCACGc                                                                                                                                     >  2:1058945/1‑148 (MQ=255)            ggccggccTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCaa                                                                                                                                  >  2:1418529/1‑150 (MQ=255)             gccggccTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAg                                                                                                                                 >  1:956092/1‑150 (MQ=255)              ccggccTGTTTGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGCGGAGCGTCAGTTGCCCGCCCTGGATGCCAACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAgg                                                                                                                                <  1:461596/150‑1 (MQ=255)              ccggccTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAgg                                                                                                                                <  1:682189/150‑1 (MQ=255)              ccggccTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAgg                                                                                                                                <  2:1158894/150‑1 (MQ=255)                   cTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGTACGCCAAGGGCcac                                                                                                                           <  2:1028951/150‑1 (MQ=255)                    tGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCcaca                                                                                                                          <  2:73273/150‑1 (MQ=255)                    tGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCcaca                                                                                                                          <  2:72138/150‑1 (MQ=255)                         gCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCAc                                                                                                                                       <  2:1252080/132‑1 (MQ=255)                         gCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCAc                                                                                                                                       >  1:1252080/1‑132 (MQ=255)                             tCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAgg                                                                                                                 <  1:518113/150‑1 (MQ=255)                             tCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAgg                                                                                                                 <  1:560008/150‑1 (MQ=255)                             tCGGGGGCTTGGCAACGAAAGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAgg                                                                                                                 >  1:877298/1‑150 (MQ=255)                                     ttGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAg                                                                                                                                                 >  1:314352/1‑111 (MQ=255)                                     ttGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAg                                                                                                                                                 <  2:314352/111‑1 (MQ=255)                                          aacgaacgCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACAt                                                                                                    >  2:1010511/1‑150 (MQ=255)                                          aacgaacgCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACAt                                                                                                    >  1:1123214/1‑150 (MQ=255)                                                 gCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGTGc                                                                                             <  2:1123214/150‑1 (MQ=255)                                                      gcgcCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCg                                                                                        <  1:1418529/150‑1 (MQ=255)                                                       cgcCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCgg                                                                                       <  1:703809/150‑1 (MQ=255)                                                         ccAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTg                                                                                     >  1:1274981/1‑150 (MQ=255)                                                            ggCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAAc                                                                                  <  2:956092/150‑1 (MQ=255)                                                             gCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGc                                                                                                                              <  1:426701/105‑1 (MQ=255)                                                             gCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGc                                                                                                                              >  2:426701/1‑105 (MQ=255)                                                                   ggTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCtgtg                                                                           <  1:823337/150‑1 (MQ=255)                                                                        ccAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCa                                                                      >  1:1183690/1‑150 (MQ=255)                                                                          aGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAggtgg                                                                                                                                             <  2:243688/78‑1 (MQ=255)                                                                          aGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAggtgg                                                                                                                                             >  1:243688/1‑78 (MQ=255)                                                                          aGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAg                                                                    >  1:632767/1‑150 (MQ=255)                                                                          aGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAg                                                                    <  2:868102/150‑1 (MQ=255)                                                                                  cGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGaaa                                                            <  1:436621/150‑1 (MQ=255)                                                                                        ttGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGcc                                                      <  1:1100245/150‑1 (MQ=255)                                                                                             aGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCa                                                                <  2:105833/135‑1 (MQ=255)                                                                                             aGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCa                                                                >  1:105833/1‑135 (MQ=255)                                                                                             aGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCgt                                                 >  2:323755/1‑150 (MQ=255)                                                                                               cccTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTg                                                                                                            <  1:594334/89‑1 (MQ=255)                                                                                               cccTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTg                                                                                                            >  2:594334/1‑89 (MQ=255)                                                                                                 cTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTg                                             >  2:1169345/1‑150 (MQ=255)                                                                                                  tGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTgg                                            >  2:1391598/1‑150 (MQ=255)                                                                                                         cGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCtgat                                     >  1:348187/1‑150 (MQ=255)                                                                                                                  tCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGc                            <  2:1183690/150‑1 (MQ=255)                                                                                                                    aGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCt                          <  2:632767/150‑1 (MQ=255)                                                                                                                      gccgccATCGCGGCGTGCAGCCACCTCAAGGAGG‑TGCACGCCAAGCGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCa                        >  2:630915/1‑150 (MQ=255)                                                                                                                             tCGCGGCTTGCAGCCGCCTCAAGGTGG‑TGCACGCCAGGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTGGGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGCTCTCGGCCTCATCCGGGc                 >  2:1122434/1‑150 (MQ=255)                                                                                                                             ccttggcgTGCA‑CCACCTTGAGGTGGCTGCACGCCGCGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGa                                                                   >  2:377190/5‑100 (MQ=21)                                                                                                                             ccttggcgTGCA‑CCACCTTGAGGTGGCTGCACGCCGCGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGa                                                                   <  1:377190/96‑1 (MQ=21)                                                                                                                                 ggcgTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGAt             <  1:64991/150‑1 (MQ=255)                                                                                                                                 ggcgTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGAt             <  2:877298/150‑1 (MQ=255)                                                                                                                                     tGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTc         <  2:348187/150‑1 (MQ=255)                                                                                                                                          ccACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGTGGTTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGaa    >  1:100314/1‑150 (MQ=255)                                                                                                                                            aCCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGAAAc  >  2:534972/1‑150 (MQ=255)                                                                                                                                              cTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGAt             <  1:1083594/137‑1 (MQ=255)                                                                                                                                              cTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGAt             >  2:1083594/1‑137 (MQ=255)                                                                                                                                              cTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGa     >  2:419636/1‑145 (MQ=255)                                                                                                                                              cTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGa     <  1:419636/145‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                   CCAGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGAAAC  >  NC_002947/3282187‑3282475

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