Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 760,485 G→A G594S (GGC→AGC)  sucA → subunit of E1(0) component of 2‑oxoglutarate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913760,4850GA97.0% 170.0 / NA 67G594S (GGC→AGC) sucAsubunit of E1(0) component of 2‑oxoglutarate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (29/35);  minor base C (2/0);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAAGTTGAGATGAAGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGGGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCAGTCTAACGGTTCCACTTACAC  >  NC_000913/760346‑760624
                                                                                                                                           |                                                                                                                                           
cAAAGTTGAGATGAAGCGCCTGCAGGAGCTGGCGAAACTCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCg                                                                                                                                          <  2:326919/143‑1 (MQ=255)
  aaGTTGTGTTGTTGCGCCTGCAGGAGCTGGCGAAACGCATGAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGt                                                                                                                                        <  2:62979/143‑1 (MQ=255)
  aaGTTGAGATGAAGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGt                                                                                                                                        <  1:79598/143‑1 (MQ=255)
     ttGAGATGAAGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCg                                                                                                                                     >  2:202477/1‑143 (MQ=255)
     ttGAGATGAAGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGAAACGCGTCTTTTCAATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCg                                                                                                                                     >  1:259611/1‑143 (MQ=255)
      ggagaTGAAGCGCCTGCAGGAGCGGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCgg                                                                                                                                    <  1:256432/142‑1 (MQ=255)
          aggaggcgcCGGCAGGAGCTGGCGAAACGCGTCAGCGCGGTGCCGGAAGCAGTTGAAATGCAGGCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGaaaa                                                                                                                                <  2:355740/138‑1 (MQ=255)
            gAAGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGTAGCTGACCGACTGGAGCGGTGCGGACAAcc                                                                                                                              >  2:229879/1‑143 (MQ=255)
            gAAGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGCGCGGTGCGGAAAAcc                                                                                                                              >  2:713122/1‑143 (MQ=255)
            gAAGCGCCTGCAGGAACTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGATGCGGaaaaac                                                                                                                              >  2:959251/1‑141 (MQ=255)
             aaGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGACAAACTGTTCGACTGGAGCGGTGCGGCAAACCt                                                                                                                             >  2:817758/1‑143 (MQ=255)
              aGCGCCTGTGGCAGCTGTCAAAACGCATCACCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAATATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTc                                                                                                                            <  2:728816/143‑1 (MQ=255)
              aGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGGAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAGCCTc                                                                                                                            >  2:1019910/1‑143 (MQ=255)
               gcgcCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCg                                                                                                                           >  1:974570/1‑143 (MQ=255)
                     gCAGGGGCTGGGGAGACGCACCGACACGGGGCCGGAAGCAGTTGAAATGCATTCTCGCGTTGCCAAGATTCATAGCGATCGCCAGGCGATGGCTGCCGGTGAGAGACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACg                                                                                                                     <  1:489633/143‑1 (MQ=255)
                     gCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGCAAACCTCGCTTACg                                                                                                                     >  1:115938/1‑143 (MQ=255)
                     gCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCCCGCTTACg                                                                                                                     >  2:380250/1‑143 (MQ=255)
                          aGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGAGTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGGGCGGTGCGGAACACCTCGCGTACGCCACg                                                                                                                >  2:1034409/1‑143 (MQ=255)
                            cTGGCGAAACGCATCAGCACGGCGCTGGAAGCAGTTGAAATGGAGGCTCGCGTTGCCAATGTTTAGGGCGACCGCCAGGCGATGGCTGGTGGTGAGAAAGTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCt                                                                                                              <  1:428712/143‑1 (MQ=255)
                               gCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGt                                                                                                           >  1:363129/1‑143 (MQ=255)
                                cggaaCGCCTCAGCACGGTGCCGGAAGCTGTTGCCCTGCCGCCTCGCGTTGCCAAGATTGGTGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGtt                                                                                                          <  2:498575/140‑1 (MQ=255)
                                cGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGACTCGCGTTGCCAAGATTTATGGCGATCGCGAGGCGATGGCTGCCGGTGAGAAACAGTTCGCCTGGAGCGGTGGGGGACACCGCGCTTACGCCACGCCGGtt                                                                                                          >  2:49649/1‑143 (MQ=255)
                                 gAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGAGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTGAGCCACGATGGTtg                                                                                                         >  2:621298/1‑143 (MQ=255)
                                 gAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGAGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGGGCGGAGAGCCTCGCTTACGCCACGCTCGTtg                                                                                                         >  2:1014700/1‑143 (MQ=255)
                                    aCGCACCAGCACGGTGCCGGAAGAAGTTGAAATGCAGTCTCGCGTTGCCAAGAATTATGGCGATCGACAGCCGATGGCTGCCGGTGAGAAACTGTTAGACTGGAGCGGGCCGGAAAACCCAGCTTACGCAACGCCGGCagatg                                                                                                      >  2:322859/1‑143 (MQ=255)
                                     cGCATCGGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCCGGTTTGATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTtgatga                                                                                                     <  1:1126948/143‑1 (MQ=255)
                                      gCATCAGCACGGCGGCGGAAGCAGTTGAAATGCAGAACGTCATGGCTAAGATTTATGGCGATCGCGAGCCGATGGCTGCCGGTGAGAACCTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGaa                                                                                                    <  2:235286/143‑1 (MQ=255)
                                       cATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGAGTGGAGCGGTGGGGAAAACCTCGCTTACGCCACACTGATtgatgacg                                                                                                   >  2:784136/1‑141 (MQ=255)
                                           aGCACGGTGCCGGAAGCAGATGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGCTGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGACACGCTGGTTGAATGAGGCAt                                                                                               >  2:996715/1‑143 (MQ=255)
                                              aCGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTcc                                                                                            >  1:531163/1‑143 (MQ=255)
                                               cGGTGCCGGAAGCAGTTGGGATGCAGTCTCGCGTTGCGAAGATTTATGGCGCTCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCg                                                                                           <  2:331292/143‑1 (MQ=255)
                                               cGGTGCCGGAAGCAGTTGAGATGCAGTCTCGCGTGGCCAAGATTTATGGCGAGCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCg                                                                                           <  1:866190/143‑1 (MQ=255)
                                                ggTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCgg                                                                                          <  2:986921/143‑1 (MQ=255)
                                                            ggTTGAAATGCAGTCTCGCGTCGCCAAGATGTATGGCGATCGCCAGGCGATGGCGGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCgg                                                                              <  1:223362/142‑1 (MQ=255)
                                                             gTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCTACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCggg                                                                             >  1:228010/1‑143 (MQ=255)
                                                              ttGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACACAAAGCTGGTTGTTGTTGGCATTCCGGTTCGCCTGTCGGGt                                                                            >  1:880980/1‑143 (MQ=255)
                                                                 aagtGCGGTCTCGGGGTGCCAAGATGTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCTGGTTCGCCTGTCGGGTGaa                                                                         <  2:933046/140‑1 (MQ=255)
                                                                 aaaTGCGGCCTCGCGTGGCCAAGATTCATGGCGACCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGaa                                                                         <  2:914799/143‑1 (MQ=255)
                                                                 aaaTGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGGTGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGaa                                                                         <  2:835231/143‑1 (MQ=255)
                                                                    tGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGCTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGCTGAACGCATTCCGGTTCGCCTGACGGGTGCAGAc                                                                      >  2:1126060/1‑143 (MQ=255)
                                                                              cgTTGCCAAGGTTTATGGCGATCGCCATGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTcgcg                                                            <  2:578825/143‑1 (MQ=255)
                                                                                  gCCAAGCTTTATGGCGATCGCCAGGCGACGGCTCCCGGTGAGAAACTGTTCGACTGGCGCGGTGCGGAAAACCTCGCCTACGCCACCCAGGTTGAAGAAATCATGCCCGTCCCCCCTTCCGGTGAAGCCTCCGGACGCCGTAc                                                        >  2:851958/1‑143 (MQ=255)
                                                                                   ccAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAACCCTCGCTTACGCCACGCGGGTTGATGAAGGCATTCCGGTGCGCCTGTCGGGTGCAGCATCCGGTCGCggcagc                                                       >  2:48522/1‑139 (MQ=255)
                                                                                    cAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGGCTCCGGGCGCGGTACct                                                      >  1:825195/1‑143 (MQ=255)
                                                                                    cAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACct                                                      >  1:658328/1‑143 (MQ=255)
                                                                                     aaGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAACACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACctt                                                     >  1:553121/1‑143 (MQ=255)
                                                                                     aaGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGCTGCCGGCAATCCGGTTCGCCTTTCGGGTGAAGACTCCGGTCGCGGTACctt                                                     >  2:1043656/1‑143 (MQ=255)
                                                                                         tttATGGCGATCGCCCGGCGATCGCTGCCGTTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCttcttc                                                 <  1:202477/143‑1 (MQ=255)
                                                                                            aTGGCGATCTCCAGGCGATGGCTTCCGGTTAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCAc                                              <  1:266991/143‑1 (MQ=255)
                                                                                             tGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCAcc                                             >  1:68248/1‑143 (MQ=255)
                                                                                              ggCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCg                                            <  1:1111825/143‑1 (MQ=255)
                                                                                                       cAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAACCCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGt                                   <  1:622190/143‑1 (MQ=255)
                                                                                                        agccGATGGCTGCCGGTGAGAAGCTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTg                                  <  2:411769/140‑1 (MQ=255)
                                                                                                         ggCGAGGGCTGCCGGTGAGAATCTGTGCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGa                                 <  1:621253/143‑1 (MQ=255)
                                                                                                          gagATGGCTGCCGGTGAGAAACTGTTTGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGAt                                <  1:983148/141‑1 (MQ=255)
                                                                                                                 cTGCCGGTGAGAAACCGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAAc                         <  1:462297/143‑1 (MQ=255)
                                                                                                                   gacGGTTAGAATCTTTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTGGATGAAGGCATGCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCa                       <  1:531949/141‑1 (MQ=255)
                                                                                                                         gagaATCTGCTCGACTGGAGCGGTCCGGAAATCCCCGCCTACGCCACGCTGGTGGATGACGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCAGTCTaa                 <  1:983432/143‑1 (MQ=255)
                                                                                                                           gaAACTGTTCGTCTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCAGTCTAACg               <  2:974570/143‑1 (MQ=255)
                                                                                                                             aaCTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCAGTCTAACGGt             <  1:1020995/143‑1 (MQ=255)
                                                                                                                               cTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCAGTCTAACGGTTc           <  2:531163/143‑1 (MQ=255)
                                                                                                                                  ttCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGGTGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCAGTCTCACGGTTCCAc        >  2:1061374/1‑143 (MQ=255)
                                                                                                                                  ttCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGCTCCACAACCAGTCTAACGGTTCCAc        >  2:1071142/1‑143 (MQ=255)
                                                                                                                                  gtCGACTGGATCGGTGCGGAAAACCTCGCTGACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCAGTCTAACGGTTCCAc        <  1:412962/142‑1 (MQ=255)
                                                                                                                                     gACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCAGTCTAACGGTTCCACTTa     >  2:231297/1‑143 (MQ=255)
                                                                                                                                       cTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCAGTCTAACGGTTCCACTTaca   <  1:1019910/143‑1 (MQ=255)
                                                                                                                                        tGGAGCCGTGCCGAAAACCTCGCTTACGCCACGCTGGTTCATAAAGTCCTTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCAGTCTAACGGTTCCACTTacac  <  1:318255/143‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                           
CAAAGTTGAGATGAAGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGGGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCAGTCTAACGGTTCCACTTACAC  >  NC_000913/760346‑760624

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