Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,239,884	G→A	18.4%	A144V (GCG→GTG)	ldh ←	leucine dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,239,884	0	G	A	18.4%	165.2 / 27.2	76	A144V (GCG→GTG)	ldh	leucine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (31/31);  new base A (8/6);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 7.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.42e-01

CTCCGCGCCCACCGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCAAAGGCCTC  >  NC_002947/5239739‑5240032                                                                                                                                                     |                                                                                                                                                     cTCCGCGCCCACCGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCAt                                                                                                                                                   >  7:302923/1‑149 (MQ=255)    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                                                                                                                                   >  3:75307/1‑146 (MQ=255)                    aaTTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGcc                                                                                                                                   <  4:75307/146‑1 (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                                                                                                                                >  4:298725/1‑132 (MQ=255)                                     aTAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGa                                                                                                                                <  3:298725/132‑1 (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                                                                                                                                   <  8:299831/115‑1 (MQ=255)                                                   ccAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGcc                                                                                                                                   >  7:299831/1‑115 (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                                                                                                                                       >  4:173795/1‑101 (MQ=255)                                                            ggACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCACATGCGGCGATgggg                                                                                                                                       <  3:173795/102‑2 (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                                                                         <  7:290413/145‑1 (MQ=255)                                                                               tctcAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAg                                                                         >  8:290413/1‑145 (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                                                                                                >  5:194070/1‑118 (MQ=255)                                                                                   aGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCtgt                                                                                                <  6:194070/118‑1 (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:146921/145‑1 (MQ=255)                                                                                                                                           ccccATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAg                >  8:141963/1‑145 (MQ=255)                                                                                                                                           ccccATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAg                >  8:146921/1‑145 (MQ=255)                                                                                                                                           ccccATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAg                <  7:141963/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)                                                                                                                                                     |                                                                                                                                                     CTCCGCGCCCACCGCATGCAATTGCTCCGCCAGGGCATAACCAACATTGCCCAACCCCTGGACCGCTACCCGCAAGCCTCTCAGGTCATC‑‑GCTGCCCAGCCGGAACGAGGTTGTGGCACGTATGCCCGCGAACA‑CCCCCATCGCCGCATGTGGCGATGGGTCGCCCGATGCAGTGGTGCTGGTCACGTAGGGCGTGCTGTGGGCAATGCAGTCCATGTCCAAGGTCGAGGTGCCACTGTCCACGGCGATGATGAAGCGCCCATGCAAGGTATCGATGAAGCGGCCAAAGGCCTC  >  NC_002947/5239739‑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: ‑