Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,734,032 G→T 7.2% R535S (CGC→AGC)  sucA ← 2‑oxoglutarate dehydrogenase subunit E1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,734,0320GT 7.2% 208.6 / 5.9 70R535S (CGC→AGC) sucA2‑oxoglutarate dehydrogenase subunit E1
Reads supporting (aligned to +/- strand):  ref base G (30/35);  new base T (2/3);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-01

CTGACGCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCGGCATACAGCTCACGGGTAGTGCGCT  >  NC_002947/4733892‑4734178
                                                                                                                                            |                                                                                                                                                  
cTGACGCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCt                                                                                                                                            >  4:150230/1‑149 (MQ=255)
  gACGCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTcc                                                                                                                                          >  3:334422/1‑149 (MQ=255)
   aCGCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACCGCGAGCTGTTCgtcga                                                                                                                                           <  4:186408/147‑5 (MQ=255)
   aCGCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTcc                                                                                                                                          <  5:345248/148‑1 (MQ=255)
   aCGCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTcc                                                                                                                                          <  3:232652/148‑1 (MQ=255)
   aCGCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCt                                                                                                                                         <  8:78086/149‑1 (MQ=255)
    cGCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCtt                                                                                                                                        <  5:329541/149‑1 (MQ=255)
       tGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGAc                                                                                                                                     >  5:233359/1‑149 (MQ=255)
            cACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACCAACAGCTCGCGCTTCGCCTCCTTGACCAGGc                                                                                                                                >  4:113611/1‑149 (MQ=255)
                agctcTTCCGGCAGCTGGAGCTGCTTGGCCGACAGATCGTGCTGTGTCTTGAGGTCGAAACGGGTGTCGTGTCGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCtttt                                                                                                                            <  6:234101/145‑1 (MQ=255)
                aaaCCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCtttt                                                                                                                            <  3:194365/149‑1 (MQ=255)
                 atccTTCCGGCCGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCttttt                                                                                                                           <  2:294463/147‑1 (MQ=255)
                 aagcTTCCGGCAGCTCGAGCTGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGTGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCttttt                                                                                                                           <  7:127924/146‑1 (MQ=255)
                 aaCCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACCGCGAGCTGTTCgtc                                                                                                                                             >  7:201170/1‑129 (MQ=255)
                 aaCCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACCGCGAGCTGTTCgtc                                                                                                                                             <  8:201170/131‑3 (MQ=255)
                   ccTTCCGGCAGCTCGAGCAGCTTGGCCGACAGTTCGTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGATGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTAc                                                                                                                         <  8:192188/149‑1 (MQ=255)
                    cTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTAcc                                                                                                                        <  3:38813/149‑1 (MQ=255)
                    cTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTAcc                                                                                                                        <  3:34132/149‑1 (MQ=255)
                    cTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTAc                                                                                                                         >  2:83290/1‑148 (MQ=255)
                     ttCCGGCAGCTCGAGCAGCTTGGCCGACACATCCTGCAGGGTCTTGAGGTCGAAACGGGTGACGTGACGCGAGGTCCAGGCATGGCCCAGGTACGGACGCAAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCtttt                                                                                                                            >  5:234101/1‑144 (MQ=255)
                       ccGGCAGCTCGCCCAGCTTGGCCGAACGATCCTTCAGGGTCTTGAGGTCGCAACGGGTGGCGTGACGCGCGGTCCAGGCATGGCCCATGGACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTTACCAGGCTTTTTACCa                                                                                                                       <  4:62101/147‑1 (MQ=255)
                       ccGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCa                                                                                                                       <  4:334422/147‑1 (MQ=255)
                       ccGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCa                                                                                                                       <  6:319961/147‑1 (MQ=255)
                          gCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCCGGTACGGACGCCAGTCGCCGAACAGGTCGCGGTGCGGCTCCTGGCCGAGGCGTTTTGCCACGTGc                                                                                                                  >  4:8164/1‑149 (MQ=255)
                          gCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCTTTGACCAGGCTTTTTACCACGTg                                                                                                                   >  4:57773/1‑148 (MQ=255)
                                            ccGACAGATCCTGCAGGGTCTTGAGGTCGATACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTa                                                                                                                          >  1:35759/1‑123 (MQ=255)
                                            ccGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTa                                                                                                                          <  2:35759/123‑1 (MQ=255)
                                              gACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCttttt                                                                                                                           <  5:58866/120‑1 (MQ=255)
                                              gACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCttttt                                                                                                                           >  6:58866/1‑120 (MQ=255)
                                              gACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGtt                                                                                              >  1:198238/1‑149 (MQ=255)
                                                 aGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCg                                                                                           <  2:198743/149‑1 (MQ=255)
                                                   aTCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCgg                                                                                          >  8:130272/1‑148 (MQ=255)
                                                      cTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACt                                                                                      <  1:222335/149‑1 (MQ=255)
                                                         cAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACtcgt                                                                                   >  5:328339/1‑149 (MQ=255)
                                                         cAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACtcgt                                                                                   <  6:233359/149‑1 (MQ=255)
                                                            ggTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGa                                                                                >  6:79060/1‑149 (MQ=255)
                                                              tCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATc                                                                              <  3:303771/149‑1 (MQ=255)
                                                                ttGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCtt                                                                            <  3:150230/149‑1 (MQ=255)
                                                                   aGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTgg                                                                          <  2:198238/148‑1 (MQ=255)
                                                                        gAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTcc                                                                                                                                          >  4:41115/1‑79 (MQ=255)
                                                                        gAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTcc                                                                                                                                          <  3:41115/79‑1 (MQ=255)
                                                                        gAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAg                                                                    >  1:27279/1‑149 (MQ=255)
                                                                        gAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAg                                                                    >  1:140849/1‑149 (MQ=255)
                                                                            cGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAgcgcg                                                                <  2:140849/149‑1 (MQ=255)
                                                                                  tCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGc                                                          >  7:325423/1‑149 (MQ=255)
                                                                                  gcGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGc                                                          <  5:33437/148‑1 (MQ=255)
                                                                                    gTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGt                                                        <  2:27279/149‑1 (MQ=255)
                                                                                    gTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGt                                                        >  2:295426/1‑149 (MQ=255)
                                                                                     tGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTc                                                       >  5:366494/1‑149 (MQ=255)
                                                                                        cgcgcgGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGAt                                                    <  7:130272/149‑1 (MQ=255)
                                                                                         gcgcgGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATg                                                   >  2:140856/1‑149 (MQ=255)
                                                                                          cgcgGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGc                                                  <  6:328339/149‑1 (MQ=255)
                                                                                           gcgGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGTCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTg                                                                                             <  1:340475/105‑1 (MQ=255)
                                                                                           gcgGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGACTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTg                                                                                             >  2:340475/1‑105 (MQ=255)
                                                                                            cgGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCTTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGa                                                <  3:57773/149‑1 (MQ=255)
                                                                                                 cAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGc                                           >  5:38552/1‑149 (MQ=255)
                                                                                                  aGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGcc                                          >  4:214676/1‑149 (MQ=255)
                                                                                                        ggCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGAt                                     >  7:30522/1‑148 (MQ=255)
                                                                                                         gCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTg                                                                                                         <  5:44438/79‑1 (MQ=255)
                                                                                                         gCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTg                                                                                                         >  6:44438/1‑79 (MQ=255)
                                                                                                         gCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCa                                   <  3:113611/149‑1 (MQ=255)
                                                                                                           ccAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGAGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAgc                                 >  3:285961/1‑149 (MQ=255)
                                                                                                                  cGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCg                           <  5:79060/148‑1 (MQ=255)
                                                                                                                    gACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCGGCa                        <  1:295426/149‑1 (MQ=255)
                                                                                                                           gTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCGGCATACAGCt                 >  1:41045/1‑149 (MQ=255)
                                                                                                                            tCGACCAACAGCTCGAGGTTCGGCTCCTTGACCAGGCTTTTTACAACCTGCAGGCCATTCTCCAGAGCGTTGCCGTACTCGTCGATCTTGGCCTGAGAGCGCTCGGCGTCGATGCGACCGGACTGGATCAGATGCTCGGCATACAGCt                 >  8:82218/1‑148 (MQ=255)
                                                                                                                            gagccgaACCGCGAGCTGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCGGCATACAGCt                 >  4:28077/5‑148 (MQ=255)
                                                                                                                             agccgaACCGCGAGCTGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCGGCATACAGCTCa               <  3:28077/146‑1 (MQ=255)
                                                                                                                                 gaACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCGGCATACAGCTCACggg           >  4:317811/1‑149 (MQ=255)
                                                                                                                                   aCAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCGGCATACAGCTCACGGGTa         <  3:8164/149‑1 (MQ=255)
                                                                                                                                          gcgGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCGGCATACAGCTCACGGGTAGTgcgc   <  1:140856/148‑1 (MQ=255)
                                                                                                                                          gcgGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCGGCATACAGCTCACGGGTAGTTCgct  >  6:186098/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
CTGACGCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCGGCATACAGCTCACGGGTAGTGCGCT  >  NC_002947/4733892‑4734178

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