Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,185,122 C→T 5.2% L446L (CTG→CTA)  aceK ← isocitrate dehydrogenase kinase/phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,185,1220CT 5.2% 239.4 / 2.1 78L446L (CTG→CTA) aceKisocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base C (35/38);  new base T (2/2);  total (38/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

CGGAATGTGGCGGAAATTCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTCCAACAGTGGCGGATCAGCACCGTGTCGCCCTCCAGCGCCACCGTCGACGGCGCCACCTCCA  >  NC_002947/5184978‑5185270
                                                                                                                                                |                                                                                                                                                    
cGGAATGTGGCGGAAATTCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGc                                                                                                                                                  <  1:189704/149‑1 (MQ=255)
 ggAATGTGGCGGAAATTCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCTAGGcc                                                                                                                                                 >  4:267689/1‑149 (MQ=255)
 ggAATGTGGCGGAAATTCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCTAGGcc                                                                                                                                                 >  4:262383/1‑149 (MQ=255)
    aTGTGGCGGAAATTCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTa                                                                                                                                              <  7:229598/149‑1 (MQ=255)
    aTGTGGCGGAAATTCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTa                                                                                                                                              <  4:179992/149‑1 (MQ=255)
     tgtgGCGGAAATTCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAg                                                                                                                                             <  3:309701/149‑1 (MQ=255)
               aTTCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAgcg                                                                                                                                   <  7:349983/149‑1 (MQ=255)
               aTTCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAgcg                                                                                                                                   >  4:324009/1‑149 (MQ=255)
                ttCACCTCGGACAGGAAGCTGATTTCGTCATAGTCCTAGAACACGACCCGGCCATCACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAgcgc                                                                                                                                  <  8:172245/149‑1 (MQ=255)
                   aCCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCCCCGCCAACTGCTTGATGGCAAGGCCTTAGTCCTCCAGCGCCTc                                                                                                                               >  2:89453/1‑149 (MQ=255)
                    ccTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTc                                                                                                                               >  5:294972/1‑148 (MQ=255)
                     cTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCa                                                                                                                             <  3:174269/149‑1 (MQ=255)
                      tCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAg                                                                                                                            >  3:254669/1‑149 (MQ=255)
                          tCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCAAGCAc                                                                                                                         >  8:331443/1‑148 (MQ=255)
                            aGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCTAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTg                                                                                                                      <  3:267689/149‑1 (MQ=255)
                            aGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCTAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTg                                                                                                                      <  3:262383/149‑1 (MQ=255)
                             ggAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCAAACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCAGTAGTCCTCCAGCGCCTCCAGCACCTgg                                                                                                                     >  4:258976/1‑149 (MQ=255)
                              gAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGc                                                                                                                    >  4:240893/1‑149 (MQ=255)
                                     tattCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCACACCAAAGTTGTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCg                                                                                                               <  5:231407/145‑1 (MQ=255)
                                      tttCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGCTGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTg                                                                                                            <  1:167931/149‑1 (MQ=255)
                                           tCATAGTCGTAGAACACGACCCGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTg                                                                                                       >  5:263750/1‑149 (MQ=255)
                                            cATAGTCGTAGAACACGACCCGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTg                                                                                                       >  6:60763/1‑148 (MQ=255)
                                             aTAGTCGTAGAACACGACCCGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGc                                                                                                      <  1:334623/148‑1 (MQ=255)
                                             aTAGTCGTAGAACACGACCCGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCCGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGc                                                                                                      <  3:180873/148‑1 (MQ=255)
                                               aGTCGTAGAACACGACCCGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGcc                                                                                                         <  1:235728/143‑1 (MQ=255)
                                               aGTCGTAGAACACGACCCGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGcc                                                                                                         >  2:235728/1‑143 (MQ=255)
                                               aGTCGTAGAACACGACCCGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTc                                                                                                    >  1:46359/1‑148 (MQ=255)
                                                gTCGTAGAACACGACCCGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCa                                                                                                  <  6:294972/149‑1 (MQ=255)
                                                 tCGTAGAACACGACCCGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAg                                                                                                 >  6:231796/1‑149 (MQ=255)
                                                             aCCCGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCaa                                                                                     >  5:100020/1‑149 (MQ=255)
                                                              cccGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCCCCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAAc                                                                                    >  8:197859/1‑149 (MQ=255)
                                                               ccGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCg                                                                                                                                                <  7:291445/88‑1 (MQ=255)
                                                               ccGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCg                                                                                                                                                >  8:291445/1‑88 (MQ=255)
                                                               ccGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTc                                                                                                    >  4:106266/1‑132 (MQ=255)
                                                               ccGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTc                                                                                                    <  3:106266/132‑1 (MQ=255)
                                                               ccGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACg                                                                                   <  3:324009/149‑1 (MQ=255)
                                                                  gCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCg                                                                                >  6:61245/1‑149 (MQ=255)
                                                                    cATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTc                                                                              <  3:258976/149‑1 (MQ=255)
                                                                     aTGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCa                                                                             <  1:89453/149‑1 (MQ=255)
                                                                        aCGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATg                                                                           <  5:60763/148‑1 (MQ=255)
                                                                           ggTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGcgcc                                                                       <  6:99578/149‑1 (MQ=255)
                                                                            gTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAgg                                                                                          <  4:289231/129‑1 (MQ=255)
                                                                            gTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAgg                                                                                          >  3:289231/1‑129 (MQ=255)
                                                                               acacCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCt                                                                    <  5:231796/148‑1 (MQ=255)
                                                                                   cAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTc                                                               >  6:365531/1‑149 (MQ=255)
                                                                                   cAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTc                                                               >  3:186288/1‑149 (MQ=255)
                                                                                    aaaGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCCCCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCACCGGCGTCATGCGCCGCTCGGTc                                                               >  7:135769/1‑148 (MQ=255)
                                                                                        ttcttcAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCGGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTCCaa                                                            >  6:25650/1‑147 (MQ=255)
                                                                                             cAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTCCAACAGTGGc                                                     >  5:211956/1‑149 (MQ=255)
                                                                                              aGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTCCAACAGTGGc                                                     <  2:46359/148‑1 (MQ=255)
                                                                                              aGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTCCAACAGTGGc                                                     <  3:240893/148‑1 (MQ=255)
                                                                                              aGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTCCAACAGTGGc                                                     <  6:340554/148‑1 (MQ=255)
                                                                                                cAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTCCAACAGTGGCGGa                                                  >  6:379642/1‑149 (MQ=255)
                                                                                                cAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTCCAACAGTGGCGGa                                                  >  2:128403/1‑149 (MQ=255)
                                                                                                 aaCATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTCCAACAGTGGCGGAt                                                 <  4:254669/149‑1 (MQ=255)
                                                                                                 aaCATGTCGCCGGCGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTCCAACAGTGGCGGAt                                                 <  7:349034/149‑1 (MQ=255)
                                                                                                   cATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTCCAACAGTGGCGGATCa                                               <  4:186288/149‑1 (MQ=255)
                                                                                                     tGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTCCAACAGTGGCGGATCAGc                                             >  3:354942/1‑149 (MQ=255)
                                                                                                       tCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTCCAACAGTGGCGGATCAGCAc                                           >  6:98558/1‑149 (MQ=255)
                                                                                                               gAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTCCAACAGTGGCGGATCAGCACCgt                                        <  1:272730/144‑1 (MQ=255)
                                                                                                               gAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTCCAACAGTGGCGGATCAGCACCgt                                        >  2:272730/1‑144 (MQ=255)
                                                                                                               gAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTCCAACAGTGGCGGATCAGCACCGTGTCGc                                   >  1:50832/1‑149 (MQ=255)
                                                                                                                 aGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTCCAACAGTGGCGGATCAGCACCGTGTCGccc                                 >  2:385209/1‑149 (MQ=255)
                                                                                                                    tGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTCCAACAGTGGCGGATCAGCACCGTGTCGCCCTcc                              <  4:343913/149‑1 (MQ=255)
                                                                                                                      ttCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTCCAACAGTGGCGGATCAGCACCGTGTCGCCCTCCAg                            >  6:332923/1‑149 (MQ=255)
                                                                                                                        cgccgccgcCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTCCAACAGTGGCGGATCAGCACCGTGTCGCCCTCCAgcg                          >  7:302409/1‑149 (MQ=255)
                                                                                                                             ccgccAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTCCAACAGTGGCGGATCAGCACCGTGTCGCCCTCCAGCGCCAcc                     <  7:331443/149‑1 (MQ=255)
                                                                                                                                  aaCTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTCCAACAGTGGCGGATCAGCACCGTGTCGCCCTCCAGCGCCACCGTCGa                >  6:187957/1‑149 (MQ=255)
                                                                                                                                    cTGCTTGATGGCCAGGGCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTCCAACAGTGGCGGATCAGCACCGTGTCGCCCTCCAGCGCCa                       <  6:136672/140‑1 (MQ=255)
                                                                                                                                    cTGCTTGATGGCCAGGCCGTAGACCTCCAGCGCCTCCAGCACCTGCCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACCGCGTCATGCGCCGCTCGGTCCAACAGTGGCGGATCAGCACCGTGTCGCCCTCCAGCGCCa                       >  5:136672/1‑140 (MQ=255)
                                                                                                                                       cTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTCCAACAGTGGCGGATCAGCACCGTGTCGCCCTCCAGCGCCACCGTCGACGgc            <  5:98558/148‑1 (MQ=255)
                                                                                                                                       cTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTCCAACAGTGGCGGATCAGCACCGTGTCGCCCTCCAGCGCCACCGTCGACGgc            <  5:379642/148‑1 (MQ=255)
                                                                                                                                          gATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTCCAACAGTGGCGGATCAGCAcc                                          >  1:340268/1‑115 (MQ=255)
                                                                                                                                          gATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTCCAACAGTGGCGGATCAGCAcc                                          <  2:340268/115‑1 (MQ=255)
                                                                                                                                          gATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTCCAACAGTGGCGGATCAGCACCGTGTCGCCCTCCAGCGCCACCGTCGACGgcgc          <  2:50832/147‑1 (MQ=255)
                                                                                                                                          gATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTCCAACAGTGGCGGATCAGCACCGTGTCGCCCTCCAGCGCCACCGTCGACGGCGcc         <  6:263750/148‑1 (MQ=255)
                                                                                                                                             ggCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTCCAACAGTGGCGGATCAGCACCGTGTCGCCCTCCAGCGCCACCGTCGACGGCGCCACCt     >  6:116158/1‑149 (MQ=255)
                                                                                                                                              gCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTCCAACAGTGGCTGATCAGCACCGTGTCGCCCTCCATCGCCACCGTCGACGGCGCCACCTc    <  5:25650/149‑1 (MQ=255)
                                                                                                                                               ccAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTCCAACAGTGGCGGATCAGCACCGTGTCGCCCTCCAGCGCCACCGTCGACGGCGCCACCTc    >  3:111029/1‑148 (MQ=255)
                                                                                                                                                cAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCTTGTACAGGTTCAACGGCGACATGCGCCGCTCGGTCCAACAGTGGCGGATCAGCACCGTGTCGCCCTCCAGCGCCACCGTCGACGGCGCCACCTCCa  <  1:128403/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CGGAATGTGGCGGAAATTCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCACACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATGTTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGCGCCTCCAGCACCTGGCCTTCGGTGGCCTGCTCCAGGTACAGGTTCAACGGCGTCATGCGCCGCTCGGTCCAACAGTGGCGGATCAGCACCGTGTCGCCCTCCAGCGCCACCGTCGACGGCGCCACCTCCA  >  NC_002947/5184978‑5185270

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: ‑