Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,004,529	A→T	11.8%	Q184L (CAG→CTG)	PP_0866 →	transmembrane sensor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,004,529	0	A	T	11.8%	181.8 / 13.7	67	Q184L (CAG→CTG)	PP_0866	transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base A (34/25);  new base T (3/5);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

CCGGATGGCTCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCATATAGAT  >  NC_002947/1004381‑1004677                                                                                                                                                     |                                                                                                                                                     ccGGATGGCTCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAg                                                                                                                                                     <  1:847863/150‑1 (MQ=255)         cTCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTC‑ACACTGCCCTGGGCAgtgtc                                                                                                                                            <  2:1280913/150‑2 (MQ=255)          tCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGtgtg                                                                                                                                            >  1:1021833/1‑150 (MQ=255)          tCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGtgtg                                                                                                                                            >  1:1212781/1‑150 (MQ=255)             cATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTTGTCGACACTGCCCAGGGCAGTGTGATc                                                                                                                                         <  2:104421/150‑1 (MQ=255)              aTCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCTGGGCAgtgtcgacc                                                                                                                                        <  1:1010990/150‑6 (MQ=255)                   gAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGCCGTTACGCGCCGCACAGGACGGAAACCGGCCGTTTGTGGTAGACACTGCCCAGGGCAGTGTGAGCGTcaca                                                                                                                                   >  1:1427214/1‑150 (MQ=255)                    aGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAg                                                                                                                                  >  2:237897/1‑150 (MQ=255)                     gCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAgg                                                                                                                                 >  2:604446/1‑150 (MQ=255)                              cAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGGGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCGCAGGTAAGAGCGGTCTGCGCAGATACCCGCtt                                                                                                                        >  2:760273/1‑150 (MQ=255)                               aaCACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCTGCGCCCATGACAGTACCCGGCAGATGGTTGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCttt                                                                                                                       >  1:774428/1‑150 (MQ=255)                               aaCACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCttt                                                                                                                       >  1:930969/1‑150 (MQ=255)                                  acGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGAc                                                                                                                    <  2:408096/150‑1 (MQ=255)                                       gTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTAAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTATGATCGTCACAGGTACCCGCTTTGACGTCCg                                                                                                               >  2:1304006/1‑150 (MQ=255)                                               tCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGAcc                                                                                                       <  2:1154365/150‑1 (MQ=255)                                               tCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGGCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGAcc                                                                                                       >  1:1293950/1‑150 (MQ=255)                                                   gTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCTGGGCAgtgtcga                                                                                                                                          <  2:38248/111‑4 (MQ=255)                                                   gTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCTGGGCAgtgtcga                                                                                                                                          >  1:38248/1‑108 (MQ=255)                                                        gcgcCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCa                                                                                              >  1:1355138/1‑150 (MQ=255)                                                               ccAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTAcc                                                                                                                             >  2:554836/1‑112 (MQ=255)                                                               ccAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTAcc                                                                                                                             <  1:554836/112‑1 (MQ=255)                                                                cAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGtgtg                                                                                      >  2:1569636/1‑150 (MQ=255)                                                                    ggCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCgg                                                                                  <  1:682293/150‑1 (MQ=255)                                                                     gCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCCGCCAGCACCCGTGTGGCGGt                                                                                 <  2:482302/150‑1 (MQ=255)                                                                     gCACGTGCCCCTGAACGCTGCTGAAGTGATGTTCAGCGTCGCCCATGACGGTAACCGGCCGTTTGTTGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTCACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGt                                                                                 <  1:1205688/150‑1 (MQ=255)                                                                          tGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAAc                                                                            >  2:965418/1‑150 (MQ=255)                                                                                     gCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCg                                                                    >  1:413606/1‑147 (MQ=255)                                                                                     gCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCg                                                                    <  2:413606/147‑1 (MQ=255)                                                                                        ggTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCgg                                                              <  1:376352/150‑1 (MQ=255)                                                                                        ggTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCgg                                                              <  2:774428/150‑1 (MQ=255)                                                                                          tGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGt                                                            >  1:1273518/1‑150 (MQ=255)                                                                                             aGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGcc                                                         >  2:1422720/1‑149 (MQ=255)                                                                                             aGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCa                                                         <  2:1021833/150‑1 (MQ=255)                                                                                              gTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAg                                                        <  2:1212781/150‑1 (MQ=255)                                                                                               tGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAgg                                                       >  2:45050/1‑150 (MQ=255)                                                                                                 aTGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGGCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGCACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGc                                                     >  1:1571913/1‑150 (MQ=255)                                                                                                    ttCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAg                                                  >  1:1283065/1‑150 (MQ=255)                                                                                                    ttCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACATCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCgg                                                              >  1:815098/1‑138 (MQ=255)                                                                                                    ttCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACATCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCgg                                                              <  2:815098/138‑1 (MQ=255)                                                                                                     tCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAgg                                                 <  2:1427214/150‑1 (MQ=255)                                                                                                      cAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAggg                                                <  2:1273518/150‑1 (MQ=255)                                                                                                       aGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGCGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGc                                               >  2:147732/1‑150 (MQ=255)                                                                                                         cGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGa                                             <  2:1355138/150‑1 (MQ=255)                                                                                                          gTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGAc                                            <  2:930969/150‑1 (MQ=255)                                                                                                               ccATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCATCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGcc                                       >  1:398598/1‑150 (MQ=255)                                                                                                                 aTGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACcca                                                                                                     >  1:1111819/1‑86 (MQ=255)                                                                                                                 aTGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACcca                                                                                                     <  2:1111819/86‑1 (MQ=255)                                                                                                                            cggtcttttGTGGTTGTCACTGCCCAGTGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACg                          <  2:348996/144‑1 (MQ=255)                                                                                                                            cGGCCGTTTGTGTTCGACACGGCCCTGGGCAGTGTGATCGTCACAGGTACTCTCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACATGGGTGGGTGCGGTTGCAGGGCAAGGGCGTCTCGCCGGCGCAACTGACg                          <  2:1571913/150‑1 (MQ=255)                                                                                                                            cGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACg                          >  1:539721/1‑150 (MQ=255)                                                                                                                            cGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACg                          >  2:1301929/1‑150 (MQ=255)                                                                                                                              gCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGc                        >  2:478811/1‑150 (MQ=255)                                                                                                                                     gtgGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAg                 >  1:1207519/1‑150 (MQ=255)                                                                                                                                      tgGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGAACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAgg                <  2:1207519/150‑1 (MQ=255)                                                                                                                                       ggTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAggg               >  2:73700/1‑150 (MQ=255)                                                                                                                                         tCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGt             >  2:356432/1‑150 (MQ=255)                                                                                                                                         atcacacTGCCCTGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCAc                                                                                             >  2:47868/4‑70 (MQ=39)                                                                                                                                         atcacacTGCCCTGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCAc                                                                                             <  1:47868/67‑1 (MQ=255)                                                                                                                                          cGACACTGCCCAGGGCAGTCTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGCGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCTGGCAGGGGTc            <  1:760273/150‑1 (MQ=255)                                                                                                                                             cacTGCCCTGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCa         >  2:1400066/1‑150 (MQ=255)                                                                                                                                              acTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCa                  >  1:187994/1‑140 (MQ=255)                                                                                                                                              acTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCa                  <  2:187994/140‑1 (MQ=255)                                                                                                                                                tGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCatat      >  2:1145927/1‑150 (MQ=255)                                                                                                                                                tGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCatat      >  1:4386/1‑150 (MQ=255)                                                                                                                                                   gcAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCATATAGa   >  2:1126795/2‑150 (MQ=255)                                                                                                                                                    cAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGCCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCATATAGAt  <  1:45050/150‑1 (MQ=255)                                                                                                                                                    cAGGGCAATGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCATATAGAt  <  1:1128846/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CCGGATGGCTCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCATATAGAT  >  NC_002947/1004381‑1004677

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