Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,520,860 C→T 7.2% G308S (GGT→AGT)  icd ← NADP(+)‑specific isocitrate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,520,8600CT 7.2% 234.6 / 6.1 84G308S (GGT→AGT) icdNADP(+)‑specific isocitrate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (39/38);  new base T (3/3);  total (43/41)
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.70e-01

GAGCCCGGGTTGACCTTGTCCTGCCCGGCATACTTGGGCGCAGTGCCGTGAGTGGCCTCGAACATGGCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGT  >  NC_002947/4520712‑4521008
                                                                                                                                                    |                                                                                                                                                    
gAGCCCGGGTTGACCTTGTCCTGCCCGGCATACTTGGGCGCAGTGCCGTGAGTGGCCTCGAACATGGCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCg                                                                                                                                                     <  2:1767915/150‑1 (MQ=255)
gAGCCCGGGTTGACCTTGTCCTGCCCGGCATACTTGGGCGCAGTGCCGTGAGTGGCCTCGAACATGGCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCg                                                                                                                                                     <  1:1657445/150‑1 (MQ=255)
gAGCCCGGGTTGACCTTGTCCTGCCCGCTATACTTGGGCGCAGTGCCGTGAGTGGCCTCGAACATGGCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCg                                                                                                                                                     <  1:481111/150‑1 (MQ=255)
 aGCCCGGGTTGACCTTGTCCTGCCCGGCATACTTGGGCGCAGTGCCGTGAGTGGCCTCGAACATGGCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCgt                                                                                                                                                    >  1:257518/1‑150 (MQ=255)
 aGCCCGGGTTGACCTTGTCCTGCCCGGCATACTTGGGCGCAGTGCCGTGAGTGGCCTCGAACATGGCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCgt                                                                                                                                                    >  2:265171/1‑150 (MQ=255)
   cccGGGTTGACCTTGTCCTGCCCGGCATACTTGGGCGCAGTGCCGTGAGTGGCCTCGAACATGGCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCgttg                                                                                                                                                  >  1:519934/1‑150 (MQ=255)
    ccGGGTTGACCTTGTCCTGCCCGGCACACTTGGGCGCAGTGCCGTGAGTGGCCTCGAACATGGCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGACGATAACGCCCACTTCTGCGGCCCGGGCATCGGACAGGTAGTCCCCgtaga                                                                                                                                                 >  2:391374/1‑150 (MQ=255)
     cGGGTTGACCTTGTCCTGCCCGGCATACTTGGGCGCAGTGCCGTGAGTGGCCTCGAACATGGCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCgttgag                                                                                                                                                <  2:257518/150‑1 (MQ=255)
            aCCTTGTCCTGCCCGGCATACTTGGGCGCAGTGCCGTGAGTGGCCTCGAACATGGCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGttgaggttgagg                                                                                                                                         >  1:503212/1‑150 (MQ=255)
            aCCTTGTCCTGCCCGGCATACTTGGGCGCAGTGCCGTGAGTGGCATAGATGGAGTCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTNCTGCGGCCAGGGCATCGGACAGGTAGTCACCGttgaggttgagg                                                                                                                                         <  1:593805/150‑1 (MQ=255)
              cTTGTCCTGCCCTGCATACTTGGGCGCAGTGCCGTGAGTGGCCTCGAACATGGCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGt                                                                                                                                       >  1:217997/1‑150 (MQ=255)
              cTTGTCCTGCCCGGCATACTTGGGCGCAGTGCCGTGAGTGGCCTCGAACATGGCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGt                                                                                                                                       >  1:481283/1‑150 (MQ=255)
               ttGTCCTGCCCGGCATACTTGGGCGCAGTGCCGTGAGTGGCCTCGAACATGGCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTg                                                                                                                                      >  2:560211/1‑150 (MQ=255)
               ttGTCCTGCCCGGCATACTTGGGCGCAGTGCCGTGAGTGGCCTCGAACATGGCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTg                                                                                                                                      >  1:411770/1‑150 (MQ=255)
                  tCCTGCCCGGCATACTTGGGCGCAGTGCCGTGAGTGGCCTCGAACATGGCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCg                                                                                                                                   >  1:734493/1‑150 (MQ=255)
                       cccGGCATACTTGGGCGCAGTGCCGTGAGTGGCCTCGAACATGGCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCAc                                                                                                                              >  1:1677331/1‑150 (MQ=255)
                          ggCATACTTGGGCGCAGTGCCGTGAGTGGCCTCGAACATGGCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCAcgtc                                                                                                                           <  2:41625/150‑1 (MQ=255)
                                cTTGGGCGCAGTGCCGTGAGTGGCCTCGAACATGGCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTc                                                                                                                     >  2:1444793/1‑150 (MQ=255)
                                   gggCGCAGTGCCGTGAGTGGCCTCGAACATGGCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGc                                                                                                                  <  2:1677331/150‑1 (MQ=255)
                                          gTGCCGTGAGTGGCCTCGAACATGGCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGt                                                                                                           >  2:931850/1‑150 (MQ=255)
                                             ccGTGAGTGGCCTCGAACATGGCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCAAAGTCGAGGCTGAAGGCGCGGCGCAGGTTGTTCGCGGCCGGGCGTAgg                                                                                                        >  2:958270/1‑149 (MQ=255)
                                               gTGAGTGGCCTCGAACATGGCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAg                                                                                                      <  2:503212/150‑1 (MQ=255)
                                                tGAGTGGCCTCGAACATGGCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAgg                                                                                                     <  2:217997/150‑1 (MQ=255)
                                                  aagggCCTCGAACTTGGCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGCTACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGAt                                                                                                   <  1:645657/147‑1 (MQ=255)
                                                       ccTCGAACATGGCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTg                                                                                                                                      >  2:1415073/1‑110 (MQ=255)
                                                       ccTCGAACATGGCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTg                                                                                                                                      <  1:1415073/110‑1 (MQ=255)
                                                         tCGAACATGGCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCtgctgc                                                                                            >  1:1695134/1‑150 (MQ=255)
                                                         tCGAACATGGCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCtgctgc                                                                                            >  1:64480/1‑150 (MQ=255)
                                                            aaCATGGCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTgcagc                                                                                         >  2:1583384/1‑150 (MQ=255)
                                                             aCATGGCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTgcagca                                                                                        <  2:1224973/150‑1 (MQ=255)
                                                               aTGGCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATg                                                                                      >  1:242342/1‑150 (MQ=255)
                                                                 ggCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGc                                                                                    <  1:1394123/150‑1 (MQ=255)
                                                                   ccACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGt                                                                                  >  1:1028286/1‑150 (MQ=255)
                                                                     aCGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATg                                                                                      >  1:59230/1‑144 (MQ=255)
                                                                     aCGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATg                                                                                      <  2:59230/144‑1 (MQ=255)
                                                                         tgtCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGa                                                                            <  2:411770/150‑1 (MQ=255)
                                                                         tgtCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGa                                                                            <  2:899856/150‑1 (MQ=255)
                                                                           gcGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATg                                                                          <  2:38088/149‑1 (MQ=255)
                                                                             ggACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGt                                                                                                                                       >  1:1303104/1‑87 (MQ=255)
                                                                             ggACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGt                                                                                                                                       <  2:1303104/87‑1 (MQ=255)
                                                                               aCAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATAGGAAAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATggagg                                                                      >  2:596215/1‑147 (MQ=255)
                                                                                   gTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCtt                                                                  >  1:238760/1‑150 (MQ=255)
                                                                                   gTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCtt                                                                  >  1:1695388/1‑150 (MQ=255)
                                                                                     tGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGCTCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGa                                                                <  1:1660111/150‑1 (MQ=255)
                                                                                     tGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGa                                                                <  2:1332419/150‑1 (MQ=255)
                                                                                                gCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTcg                                                     >  1:1268101/1‑150 (MQ=255)
                                                                                                gCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGAGACATGGCGTCGGCGATGGCGTCCTTGACGATGACCTcg                                                     >  1:99896/1‑150 (MQ=255)
                                                                                                  aaTGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACTGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTcgcg                                                   >  1:471330/1‑150 (MQ=255)
                                                                                                  aaTGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACTGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTcgcg                                                   >  1:471169/1‑150 (MQ=255)
                                                                                                  aaTGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACTGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTcgcg                                                   >  1:470769/1‑150 (MQ=255)
                                                                                                    tGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGAc                                                 <  2:238760/150‑1 (MQ=255)
                                                                                                             cGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTgg                                        >  2:1599532/1‑150 (MQ=255)
                                                                                                             cGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTgg                                        <  1:1599532/150‑1 (MQ=255)
                                                                                                              gCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTggg                                       <  2:1579124/150‑1 (MQ=255)
                                                                                                                 cACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTggg                                       >  2:498961/1‑147 (MQ=255)
                                                                                                                 cACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTggg                                       <  1:498961/147‑1 (MQ=255)
                                                                                                                  aCTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCCATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGTAGCCTGCCTGCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTc                                   <  1:1788672/150‑1 (MQ=255)
                                                                                                                   cTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCt                                  <  2:1695388/150‑1 (MQ=255)
                                                                                                                     tCTGCGGCCAGGGCATCGGACAGGTAGTCACTGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTg                                <  2:470769/150‑1 (MQ=255)
                                                                                                                     tCTGCGGCCAGGGCATCGGACAGGTAGTCACTGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTg                                <  2:471330/150‑1 (MQ=255)
                                                                                                                     tCTGCGGCCAGGGCATCGGACAGGTAGTCACTGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTGCTTg                                <  2:471169/150‑1 (MQ=255)
                                                                                                                      cTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGa                               <  1:931850/150‑1 (MQ=255)
                                                                                                                           gCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTcgc                                                    >  2:685426/1‑124 (MQ=255)
                                                                                                                           gCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTcgc                                                    <  1:685426/124‑1 (MQ=255)
                                                                                                                            ccAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAg                                                                                                      >  1:48902/1‑73 (MQ=255)
                                                                                                                            ccAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAg                                                                                                      <  2:48902/73‑1 (MQ=255)
                                                                                                                             cAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAgg                                                                                                     <  2:1758724/73‑1 (MQ=255)
                                                                                                                             cAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAgg                                                                                                     >  1:1758724/1‑73 (MQ=255)
                                                                                                                              aGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATc                       <  1:1583384/150‑1 (MQ=255)
                                                                                                                                ggCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCa                     <  2:734493/150‑1 (MQ=255)
                                                                                                                                  cATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAg                   <  2:1917045/150‑1 (MQ=255)
                                                                                                                                    tCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGATTCTTGAACTTCATCCAAggg                 >  2:918014/1‑150 (MQ=255)
                                                                                                                                      ggACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGgcc               >  1:1365669/1‑150 (MQ=255)
                                                                                                                                      ggACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGgcc               >  2:510993/1‑150 (MQ=255)
                                                                                                                                          aGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATg                                                                          <  2:1814033/87‑1 (MQ=255)
                                                                                                                                          aGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATg                                                                          >  1:1814033/1‑87 (MQ=255)
                                                                                                                                            gTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCAt                                                                                       >  1:1485154/1‑72 (MQ=255)
                                                                                                                                            gTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCAt                                                                                       <  2:1485154/72‑1 (MQ=255)
                                                                                                                                            gTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGc                                                                                    <  2:1024179/75‑1 (MQ=255)
                                                                                                                                            gTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGc                                                                                    >  1:1024179/1‑75 (MQ=255)
                                                                                                                                            gTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATc         >  2:1766528/1‑150 (MQ=255)
                                                                                                                                               gTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAg      <  1:1681953/150‑1 (MQ=255)
                                                                                                                                                tCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTc                                   <  2:1208525/120‑1 (MQ=255)
                                                                                                                                                tCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTc                                   >  1:1208525/1‑120 (MQ=255)
                                                                                                                                                   ccGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGt  <  1:1548700/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GAGCCCGGGTTGACCTTGTCCTGCCCGGCATACTTGGGCGCAGTGCCGTGAGTGGCCTCGAACATGGCCACGGTGTCGGACAGGTTGGCGCCCGGGGCAATGCCGATACCGCCCACTTCTGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGT  >  NC_002947/4520712‑4521008

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