Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,239,889	T→C	19.7%	P142P (CCA→CCG)	ldh ←	leucine dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,239,889	0	T	C	19.7%	152.5 / 31.8	77	P142P (CCA→CCG)	ldh	leucine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (30/31);  new base C (8/7);  total (38/39)
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.33e-01

CGCGCCCACCGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCAAAGGCCTC  >  NC_002947/5239742‑5240032                                                                                                                                                       |                                                                                                                                                cgcgCCCACCGCATGTAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCAtgtg                                                                                                                                                >  8:263370/1‑149 (MQ=255)  gcgcCCACCGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCACATGCgg                                                                                                                                               <  5:94407/149‑1 (MQ=255)      ccACCGCATGCAATTGCTCCGCCAGGGCTTAACCAACATTACCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCCATCGCCACATGCGGc                                                                                                                                              >  1:210884/1‑147 (MQ=255)      ccACCGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGACCA‑CCCCCATCGCCGCATGTGGCGAt                                                                                                                                           >  5:343504/1‑149 (MQ=255)        aCCGGATTCTATTGCTGCGCGAGGTCATAACCTACATTGCCCATCCCGTGGACGGCTACCGGCAAGCCTCTCTGTTCATC‑‑GCTGCCCAGACGTAACGAGGTTGTGGCACGTATGCCCGCTAACA‑CCCCCTTCGCCAGATGCGGCGATgg                                                                                                                                         <  4:105031/149‑1 (MQ=255)         ccGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATggg                                                                                                                                        >  4:219793/1‑149 (MQ=255)         ccGCATGCAATTGCTCCGCCAGGGCATAACCAACATTACCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACACCCCCCATCGCCACATGCGGCGAggg                                                                                                                                         <  2:210884/149‑1 (MQ=255)          cGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATCCGGCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATgg                                                                                                                                         <  1:115887/149‑1 (MQ=255)            cATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CACCCATCGCCACATGCGGCGATgggggt                                                                                                                                     <  3:338360/149‑4 (MQ=255)                 aaTTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGcc                                                                                                                                   <  4:75307/146‑1 (MQ=255)                 aaTTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGcc                                                                                                                                   >  3:75307/1‑146 (MQ=255)                       tCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAgt                                                                                                                          >  1:154755/1‑149 (MQ=255)                         cGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAgt                                                                                                                          <  8:103438/147‑1 (MQ=255)                             aGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAACCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCt                                                                                                                    <  3:147298/149‑1 (MQ=255)                               ggCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTgg                                                                                                                  >  7:34692/1‑149 (MQ=255)                                gCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGt                                                                                                                 <  3:348390/149‑1 (MQ=255)                                  aTAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGa                                                                                                                                <  3:298725/132‑1 (MQ=255)                                  aTAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGa                                                                                                                                >  4:298725/1‑132 (MQ=255)                                    aaCCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCAACGCCGCATGTGGCGATGGGTCCCCGGATGCAGTGGTGCTGGACACg                                                                                                             >  1:522/1‑149 (MQ=255)                                       cAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAg                                                                                                          >  3:324026/1‑149 (MQ=255)                                        aaCATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAgg                                                                                                         <  3:219793/149‑1 (MQ=255)                                                ccAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGcc                                                                                                                                   >  7:299831/1‑115 (MQ=255)                                                ccAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGcc                                                                                                                                   <  8:299831/115‑1 (MQ=255)                                                 cAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCtgt                                                                                                <  6:385311/149‑1 (MQ=255)                                                     cccTGGACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTgg                                                                                              >  6:278093/1‑147 (MQ=255)                                                     cccTGGACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTgg                                                                                              <  5:278093/147‑1 (MQ=255)                                                      ccTGGACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGc                                                                                            >  6:348281/1‑148 (MQ=255)                                                       cgggACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCCGCCGGCCCGTGGGTGTGGCACGTATGCCCGCGACCA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCaa                                                                                          <  6:343504/147‑1 (MQ=255)                                                         ggACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCACATGCGGCGATgggg                                                                                                                                       <  3:173795/102‑2 (MQ=255)                                                         ggACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCACATGCGGCGATgggg                                                                                                                                       >  4:173795/1‑101 (MQ=255)                                                          gACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATg                                                                                        <  8:34692/148‑1 (MQ=255)                                                             cGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGt                                                                                    <  5:348281/149‑1 (MQ=255)                                                              gCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTc                                                                                   <  8:194860/149‑1 (MQ=255)                                                              gCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGGGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTc                                                                                   <  6:287625/149‑1 (MQ=255)                                                                  cccGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCAt                                                                                <  5:209048/148‑1 (MQ=255)                                                                        aGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAg                                                                         <  2:154755/149‑1 (MQ=255)                                                                           ctctCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTc                                                                      <  8:302923/149‑1 (MQ=255)                                                                            tctcAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAg                                                                         >  8:290413/1‑145 (MQ=255)                                                                            tctcAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAg                                                                         <  7:290413/145‑1 (MQ=255)                                                                            tctcAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTc                                                                      >  1:79512/1‑148 (MQ=255)                                                                            tctcAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGAGGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCg                                                                     >  4:100890/1‑149 (MQ=255)                                                                             ctcAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATAGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTc                                                                      >  7:84460/1‑147 (MQ=255)                                                                                aGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCtgt                                                                                                <  6:194070/118‑1 (MQ=255)                                                                                aGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCtgt                                                                                                >  5:194070/1‑118 (MQ=255)                                                                                aGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAg                                                                   <  2:19910/147‑1 (MQ=255)                                                                                   tCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGc                                                               >  7:289207/1‑148 (MQ=255)                                                                                     aTC‑‑GGTGCCGAGCCGGAACGAGGTTGTGGCACGTATGCCCGAGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGGCGAGGTGCCa                                                             <  3:240881/148‑1 (MQ=255)                                                                                       atggCTGCCCAGCCGGAAGGAGGTTGTGGCACGTATGCGCGCGAACA‑GCCCCTTCGCGGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGGCa                                                             <  2:125461/145‑1 (MQ=255)                                                                                           cTGCCCACCCGCACCGCGGTGGGGGCACGTATCCCCGCTAACA‑CCCCCATCGCCGCATGTGGCTATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTc                                                        <  2:355514/149‑1 (MQ=255)                                                                                              cccAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCa                                                      <  7:268223/148‑1 (MQ=255)                                                                                               ccAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACg                                                    >  1:126401/1‑149 (MQ=255)                                                                                                 aGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGc                                                  >  2:266567/1‑149 (MQ=255)                                                                                                      gAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCgatga                                             >  6:131812/1‑149 (MQ=255)                                                                                                              ttGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAgcgc                                     >  2:113408/1‑149 (MQ=255)                                                                                                                 tgGCACGTATGCCCGCGCACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCCCGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCCCGGCGATGATGAAGCGCCCa                                  >  2:32028/1‑149 (MQ=255)                                                                                                                      cGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCaa                             <  7:263370/149‑1 (MQ=255)                                                                                                                            cccGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATc                       >  8:43467/1‑149 (MQ=255)                                                                                                                             ccGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTgg                                                                                              >  8:309272/1‑77 (MQ=255)                                                                                                                             ccGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTgg                                                                                              <  7:309272/77‑1 (MQ=255)                                                                                                                             ccGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATc                       >  2:71876/1‑148 (MQ=255)                                                                                                                                   gggcgacccATCGCCACATGCGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATg                   >  5:341576/7‑147 (MQ=255)                                                                                                                                     a‑cCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCg              <  5:131812/149‑1 (MQ=255)                                                                                                                                        ccccATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAg                <  7:141963/145‑1 (MQ=255)                                                                                                                                        ccccATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAg                >  8:146921/1‑145 (MQ=255)                                                                                                                                        ccccATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAg                >  8:141963/1‑145 (MQ=255)                                                                                                                                        ccccATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAg                <  7:146921/145‑1 (MQ=255)                                                                                                                                        ccccATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGc            >  5:311501/1‑149 (MQ=255)                                                                                                                                        ccccATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGc            <  2:126401/149‑1 (MQ=255)                                                                                                                                          ccATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGGAAGGTATCGATGAAGCGGcc           >  8:172571/1‑148 (MQ=255)                                                                                                                                          ccATCGCCACATGCGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAgg                                                                        <  4:248977/87‑1 (MQ=255)                                                                                                                                          ccATCGCCACATGCGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAgg                                                                        >  3:248977/1‑87 (MQ=255)                                                                                                                                           cATCGCCACATGCGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCa          >  3:41007/1‑148 (MQ=255)                                                                                                                                           cATCGCCACATGCGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCa          >  5:49977/1‑148 (MQ=255)                                                                                                                                            aTCGCCACATGCGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCa          >  5:199563/1‑147 (MQ=255)                                                                                                                                            aTCGCCACATGCGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCa          >  6:375410/1‑147 (MQ=255)                                                                                                                                                 cgcATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGTTGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCAAAGGCCt   <  2:522/149‑1 (MQ=255)                                                                                                                                                  gcATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCAAAGGCCTc  >  5:67094/1‑149 (MQ=255)                                                                                                                                                  acATGCGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCAAAGGCCTc  <  4:41007/148‑1 (MQ=255)                                                                                                                                                       tgGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCgat                                               >  8:65054/1‑99 (MQ=255)                                                                                                                                                       tgGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCgat                                               <  7:65054/99‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                CGCGCCCACCGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCAAAGGCCTC  >  NC_002947/5239742‑5240032

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