Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,906,508	A→T	7.7%	L112Q (CTG→CAG)	PP_4315 ←	epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,906,508	0	A	T	7.7%	226.8 / 7.3	75	L112Q (CTG→CAG)	PP_4315	epimerase
Reads supporting (aligned to +/- strand):  ref base A (38/31);  new base T (3/3);  total (41/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 = 1.00e+00

TCAACGTTTCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTGGGGGT  >  NC_002947/4906359‑4906656                                                                                                                                                      |                                                                                                                                                     tCAACGTTTCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCca                                                                                                                                                      <  1:724021/150‑1 (MQ=255) tCAACGTTTCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCca                                                                                                                                                      >  2:505851/1‑150 (MQ=255)   aaCGTTTCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCcagc                                                                                                                                                    >  2:884887/1‑150 (MQ=255)       tttCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGc                                                                                                                                                >  2:508624/1‑150 (MQ=255)           gAGCAGCATGGCGTCATCGATCACCCCGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCTGCTGGCCAtg                                                                                                                                            <  2:1244077/150‑2 (MQ=255)           gAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATc                                                                                                                                            >  1:224151/1‑150 (MQ=255)           gAGCAGCATGGCGTCATCGATCACCACGACATAATCGTAGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATc                                                                                                                                            >  2:286869/1‑150 (MQ=255)            agcagcATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCa                                                                                                                                           <  2:224151/150‑1 (MQ=255)              cagcaTGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTTCGGGTTGTTTGGCGGGGAAGTCCATGGCCTGCTGGCCAtggac                                                                                                                                         <  2:1033777/150‑5 (MQ=255)                  aTGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCTGCTGGCCAtgga                                                                                                                                          <  2:905048/145‑4 (MQ=255)                  aTGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCTGCTGGCCAtgga                                                                                                                                          >  1:905048/1‑142 (MQ=255)                     gCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACgg                                                                                                                                  <  1:508624/150‑1 (MQ=255)                        tCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTg                                                                                                                               >  2:417586/1‑150 (MQ=255)                             gATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGAcgcg                                                                                                                          >  1:245745/1‑150 (MQ=255)                                 accacGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACgg                                                                                                                                  >  1:1367499/1‑138 (MQ=255)                                 accacGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACgg                                                                                                                                  <  2:1367499/138‑1 (MQ=255)                                  ccacGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGACACGGCTGGCGCGCAACt                                                                                                                     >  2:1170202/1‑150 (MQ=255)                                   cacGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTc                                                                                                                    >  2:1104497/1‑150 (MQ=255)                                         aTAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCa                                                                                                                                           <  1:140900/121‑1 (MQ=255)                                         aTAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCa                                                                                                                                           >  2:140900/1‑121 (MQ=255)                                          tAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTg                                                                                                             <  2:1336611/150‑1 (MQ=255)                                           aaTCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGACGTGCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGc                                                                                                            <  1:1170202/150‑1 (MQ=255)                                           aaTCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGAcgc                                                                                                                           >  1:1318441/1‑135 (MQ=255)                                           aaTCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGAcgc                                                                                                                           <  2:1318441/135‑1 (MQ=255)                                             tcgtcgCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACgcgc                                                                                                                         <  2:522990/135‑1 (MQ=255)                                             tcgtcgCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACgcgc                                                                                                                         >  1:522990/1‑135 (MQ=255)                                             tcgtcgCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCgg                                                                                                          <  1:1104497/150‑1 (MQ=255)                                                   cTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTg                                                                                                    >  2:776197/1‑150 (MQ=255)                                                   cTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTg                                                                                                    <  1:108300/150‑1 (MQ=255)                                                    tAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACgcgc                                                                                                                         <  2:1411684/128‑1 (MQ=255)                                                    tAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACgcgc                                                                                                                         >  1:1411684/1‑128 (MQ=255)                                                        gTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCa                                                                                                                                              >  2:310495/1‑103 (MQ=255)                                                        gTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCa                                                                                                                                              <  1:310495/103‑1 (MQ=255)                                                             ccGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCa                                                                                          <  1:776197/150‑1 (MQ=255)                                                              cGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCAc                                                                                         >  1:688340/1‑150 (MQ=255)                                                                ccTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGTAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACtt                                                                                       >  1:987505/1‑150 (MQ=255)                                                                 cTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTAt                                                                                                      <  1:1045140/134‑1 (MQ=255)                                                                 cTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTAt                                                                                                      >  2:1045140/1‑134 (MQ=255)                                                                  ttCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTgc                                                                                     <  1:402114/150‑1 (MQ=255)                                                                     gCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGcc                                                                                  >  1:296504/1‑150 (MQ=255)                                                                            tCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCAc                                                                           >  2:359043/1‑150 (MQ=255)                                                                            tCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCAc                                                                           >  2:627877/1‑150 (MQ=255)                                                                            tCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGATTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCAc                                                                           >  1:1389947/1‑150 (MQ=255)                                                                                gCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAg                                                                       >  1:511534/1‑150 (MQ=255)                                                                                   cAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTg                                                                    >  2:844879/1‑150 (MQ=255)                                                                                        cttgcctaaGTCCCGCTGTGCTGCCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCACCCCCACGTCTGACGCGCAACTCGCCACTCCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGa                                                               <  2:559934/143‑1 (MQ=255)                                                                                          tGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAAc                                                             >  1:1281273/1‑150 (MQ=255)                                                                                          tGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAAc                                                             >  1:775828/1‑150 (MQ=255)                                                                                           gCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACt                                                                                        >  2:403941/1‑122 (MQ=255)                                                                                           gCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACt                                                                                        <  1:403941/122‑1 (MQ=255)                                                                                           gCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACa                                                            >  1:1121031/1‑150 (MQ=255)                                                                                                   cccGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTaa                                                    <  2:775828/150‑1 (MQ=255)                                                                                                     cGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGc                                                  >  2:524468/1‑150 (MQ=255)                                                                                                       cTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTAAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCg                                                >  2:75274/1‑150 (MQ=255)                                                                                                        tGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGa                                               <  2:245745/150‑1 (MQ=255)                                                                                                          gCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAgg                                             <  1:417586/150‑1 (MQ=255)                                                                                                           cATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGc                                            <  1:844879/150‑1 (MQ=255)                                                                                                            aTGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGcc                                           >  1:222485/1‑150 (MQ=255)                                                                                                             tGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCa                                                         <  2:326277/135‑1 (MQ=255)                                                                                                             tGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCa                                                         >  1:326277/1‑135 (MQ=255)                                                                                                                   ccGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCaa                                         <  2:1409225/145‑1 (MQ=255)                                                                                                                   ccGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCaa                                         >  1:1409225/1‑145 (MQ=255)                                                                                                                     gCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCa                                  <  2:688340/150‑1 (MQ=255)                                                                                                                      cTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCAt                                 <  1:524468/150‑1 (MQ=255)                                                                                                                       tACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATg                                >  1:1208176/1‑150 (MQ=255)                                                                                                                       tACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATg                                <  1:855439/150‑1 (MQ=255)                                                                                                                            gttgttTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGc                           >  2:545424/1‑150 (MQ=255)                                                                                                                             ttgttTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGca                          <  2:1281273/150‑1 (MQ=255)                                                                                                                                     cGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGa                  >  1:1300169/1‑150 (MQ=255)                                                                                                                                          aaGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTcc             >  1:1120865/1‑150 (MQ=255)                                                                                                                                           aGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCa            <  1:545424/150‑1 (MQ=255)                                                                                                                                            gctgaTGGCCTGCTGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCAc           >  1:1278147/5‑150 (MQ=255)                                                                                                                                               cATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGt        >  2:1311322/1‑150 (MQ=255)                                                                                                                                               cATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGt        <  2:1389947/150‑1 (MQ=255)                                                                                                                                                aTGGCCTGCTGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTg       >  1:140368/1‑150 (MQ=255)                                                                                                                                                aTGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTg       <  2:296504/150‑1 (MQ=255)                                                                                                                                                 tGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTgg      <  1:359043/150‑1 (MQ=255)                                                                                                                                                 tGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTgg      >  1:640878/1‑150 (MQ=255)                                                                                                                                                     cagcagGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTGGGGGt  <  1:404336/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     TCAACGTTTCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTGGGGGT  >  NC_002947/4906359‑4906656

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