Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,791,834	G→C	7.1%	G162R (GGT→CGT)	rseP →	regulatory intramembrane protein RIP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,791,834	0	G	C	7.1%	211.9 / 6.2	71	G162R (GGT→CGT)	rseP	regulatory intramembrane protein RIP
Reads supporting (aligned to +/- strand):  ref base G (35/31);  new base C (3/2);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.11e-01

GCTATCCTGTTTTTCTGGGTGCTGGCGATGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCAAGGGCGCCGACGAGCCGGAC  >  NC_002947/1791687‑1791980                                                                                                                                                     |                                                                                                                                                   gCTATCCTGTTTTTCTGGGTGCTGGCGATGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACgg                                                                                                                                                   <  1:91796/149‑1 (MQ=255)     tCCTGTTTTTCTGGGTGCTGGCGATGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTg                                                                                                                                                <  6:252395/148‑1 (MQ=255)       cTGTTTTTCTGGGTGCTGGCGATGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTc                                                                                                                                             >  5:35506/1‑149 (MQ=255)          tttttCTGGGTGCTGGCGATGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTcg                                                                                                                                            <  6:224245/147‑1 (MQ=255)          tttttCTGGGTGCTGGCGATGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGCCAAGCCGACCAACGGTTGGTcgg                                                                                                                                           >  5:22711/1‑148 (MQ=255)          tttttCTGGGTGCTGGCGATCCTTGGTACCCTGCAAATCCTCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTcgt                                                                                                                                           <  6:22711/148‑2 (MQ=255)            tttCTGGGTGCTGGCGATGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTcggcg                                                                                                                                         <  4:65680/148‑1 (MQ=255)            tttCTGGGTGCTGGCGATGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTcggcg                                                                                                                                         <  5:6725/148‑1 (MQ=255)             ttCTGGGTGCTGGCGATGCTGTGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACCGTTGGTcggctt                                                                                                                                        <  5:98644/148‑3 (MQ=255)                tGGGTGCTGGCGATGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCaa                                                                                                                                    >  3:366644/1‑149 (MQ=255)                     gCTGGCGATGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTg                                                                                                                                >  5:390856/1‑148 (MQ=255)                           gATGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTgg                                                                                                                         >  4:138749/1‑149 (MQ=255)                            ttGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGt                                                                                                                        >  8:276699/2‑149 (MQ=255)                             tGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACTGTGGCAGCCTGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGtt                                                                                                                       >  6:346878/1‑149 (MQ=255)                                tGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAG‑CATCTGCAGGGCTGCCCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGCCCAACGGGTGGTCGGCGGTCAATTTGCAGCTGGTtcgt                                                                                                                    >  6:368629/1‑149 (MQ=255)                                  ggTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTtcgtc                                                                                                                   <  2:78810/148‑1 (MQ=255)                                       ccAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTgg                                                                                                             <  8:142938/149‑1 (MQ=255)                                         aGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGc                                                                                                           >  1:8577/1‑149 (MQ=255)                                         aGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGc                                                                                                           <  3:120678/149‑1 (MQ=255)                                                  gCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCgg                                                                                                   >  8:41040/1‑148 (MQ=255)                                                     cGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTAcc                                                                                               <  5:353484/149‑1 (MQ=255)                                                       gTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCtt                                                                                             >  8:269429/1‑149 (MQ=255)                                                          aTCGGCGCGGTCGACAGTGTCAGCCTGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGGCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGc                                                                                           <  6:13057/148‑1 (MQ=255)                                                          aTCGGCGCGGTCGACAGTGGCAGCCTGGCAGTCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACCGTTGGTcggctt                                                                                                                                        >  1:41566/1‑102 (MQ=255)                                                          aTCGGCGCGGTCGACAGTGGCAGCCTGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACCGTTGGTcggctt                                                                                                                                        <  2:41566/103‑3 (MQ=255)                                                           tCGGCGCGGTCGACAGTGGCAGCCTGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCa                                                                                          <  4:9159/148‑1 (MQ=255)                                                           tCGGCGCGGTCGACAGTGGCAGCCTGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCa                                                                                          <  5:145745/148‑1 (MQ=255)                                                                gcgGTCGACAGTGGCAGCCTGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTgg                                                                                    <  6:390856/149‑1 (MQ=255)                                                                    tCGACAGTGGCAGCCTGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTg                                                                                >  3:161247/1‑149 (MQ=255)                                                                      gACAGTGGCAGCCTGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTgcg                                                                              >  1:165945/1‑149 (MQ=255)                                                                             gCAGCCTGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTg                                                                                                                                <  6:243839/92‑1 (MQ=255)                                                                             gCAGCCTGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAg                                                                       >  5:282922/1‑149 (MQ=255)                                                                             gCAGCCTGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAg                                                                       >  1:303043/1‑149 (MQ=255)                                                                                   tGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCACGATGAGGgcgc                                                                  >  6:296142/1‑148 (MQ=255)                                                                                               gCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAgcg                                                      >  8:127508/1‑149 (MQ=255)                                                                                               gCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAgcg                                                      >  2:329678/1‑149 (MQ=255)                                                                                                 aGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAgcgc                                                     >  8:364707/1‑148 (MQ=255)                                                                                                   ggCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACAAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAg                                                  >  3:61766/1‑149 (MQ=255)                                                                                                      tGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTg                                               >  2:231546/1‑149 (MQ=255)                                                                                                      tGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAAATGGTTCGTCGTCTGGGCAACAGCGCTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTg                                               >  1:314561/1‑149 (MQ=255)                                                                                                         ccGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAg                                            <  3:138749/149‑1 (MQ=255)                                                                                                         ccGCGGGTCAGGAAATTGTCTCTGTGCATCGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGGCTGGCCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAg                                            <  7:300672/149‑1 (MQ=255)                                                                                                           gcgGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGt                                          >  5:298033/1‑149 (MQ=255)                                                                                                           gcgGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGACAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGt                                          <  3:25815/149‑1 (MQ=255)                                                                                                                tCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGc                                     >  1:163225/1‑149 (MQ=255)                                                                                                                 cAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGagag                                                                                                      >  5:255961/1‑83 (MQ=255)                                                                                                                 cAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGagag                                                                                                      <  6:255961/83‑1 (MQ=255)                                                                                                                   ggAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTg                                   >  7:309810/1‑148 (MQ=255)                                                                                                                    gAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGcc                                                    <  6:175065/130‑1 (MQ=255)                                                                                                                    gAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGcc                                                    >  5:175065/1‑130 (MQ=255)                                                                                                                    gAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGa                                 <  4:161247/149‑1 (MQ=255)                                                                                                                     aaaTTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGAc                                <  7:269429/149‑1 (MQ=255)                                                                                                                      aaTTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGa                                 >  5:17052/1‑147 (MQ=255)                                                                                                                        ttGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGt                             >  8:38180/1‑149 (MQ=255)                                                                                                                          gTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGtt                            <  5:346878/148‑1 (MQ=255)                                                                                                                                 tgGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCAAg                    <  1:33320/149‑1 (MQ=255)                                                                                                                                 tgGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCAAg                    <  6:35506/149‑1 (MQ=255)                                                                                                                                   gATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCAAgg                   <  7:276699/148‑1 (MQ=255)                                                                                                                                    aTGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCAAGGgc                 <  6:298033/149‑1 (MQ=255)                                                                                                                                    aTGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCAAGGgc                 <  1:231546/149‑1 (MQ=255)                                                                                                                                    aTGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCAAGGgc                 >  6:45307/1‑149 (MQ=255)                                                                                                                                     tGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCAAGGgcg                <  2:303043/149‑1 (MQ=255)                                                                                                                                     tGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCAAGGgcg                <  6:282922/149‑1 (MQ=255)                                                                                                                                      ggCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCAAGGgcgc               >  1:8428/1‑149 (MQ=255)                                                                                                                                       gCAAGCCGCCCATCGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCAAGGgcgc               <  1:329678/148‑1 (MQ=255)                                                                                                                                         aaGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCAAGGGCGCcga            >  4:284452/1‑149 (MQ=255)                                                                                                                                         aaGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCAAGGGCGCcga            >  1:179220/1‑149 (MQ=255)                                                                                                                                          cgCCGACCAACCGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCAAGGGCGCcgac           >  5:183620/2‑149 (MQ=255)                                                                                                                                           gCCGACCAACCGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCAAGGGCGCcgac           >  1:149646/1‑148 (MQ=255)                                                                                                                                              gACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCAAGGGCGCCGACGAGc       >  6:350101/1‑149 (MQ=255)                                                                                                                                              gACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCAAGGGCGCCGACGAGc       <  4:366644/149‑1 (MQ=255)                                                                                                                                                   aCGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCAAGGGCGCCGACGAGCCGGAc  >  6:367561/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GCTATCCTGTTTTTCTGGGTGCTGGCGATGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAG‑CATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCAAGGGCGCCGACGAGCCGGAC  >  NC_002947/1791687‑1791980

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