Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,734,025 A→C 5.2% L537R (CTG→CGG)  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,0250AC 5.2% 246.4 / 3.8 77L537R (CTG→CGG) sucA2‑oxoglutarate dehydrogenase subunit E1
Reads supporting (aligned to +/- strand):  ref base A (38/35);  new base C (2/2);  total (40/37)
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.94e-01

CTTCGACACCTGACGCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCGGCATACAGCTCACGGGTAG  >  NC_002947/4733883‑4734172
                                                                                                                                               |                                                                                                                                                   
cTTCGACACCTGACGCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCGTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTcgc                                                                                                                                               >  3:101371/1‑149 (MQ=255)
 ttCGACACCTGACGCTGGACCACGGAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTcgcg                                                                                                                                              >  8:167454/1‑149 (MQ=255)
  tCGACACCTGACGCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCgg                                                                                                                                             >  7:204847/1‑149 (MQ=255)
     acacCTGACGCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTc                                                                                                                                          <  2:289083/149‑1 (MQ=255)
        ccTGACGCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGc                                                                                                                                       <  6:164364/149‑1 (MQ=255)
             cGCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCtt                                                                                                                                  <  7:88844/149‑1 (MQ=255)
              gCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCtt                                                                                                                                  >  1:170213/1‑148 (MQ=255)
              gCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCtt                                                                                                                                  >  7:17803/1‑148 (MQ=255)
                  gACCACGGAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGAc                                                                                                                               <  7:167454/147‑1 (MQ=255)
                    ccACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACCGCGAGCTGTTCgc                                                                                                                                        <  1:103842/136‑2 (MQ=255)
                    ccACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACCGCGAGCTGTTCgc                                                                                                                                        >  2:103842/1‑135 (MQ=255)
                    ccACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAgg                                                                                                                           >  3:66210/1‑149 (MQ=255)
                     cACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTcgc                                                                                                                                               >  8:141296/1‑128 (MQ=255)
                      aCGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCt                                                                                                                         >  3:56108/1‑149 (MQ=255)
                       cGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCtt                                                                                                                        >  5:413002/1‑149 (MQ=255)
                             cTTCCGGCAGCTCGAGCAGCTTGGCCGAGCAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTAc                                                                                                                   >  1:24083/1‑149 (MQ=255)
                                ccGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCAc                                                                                                                <  2:170213/148‑1 (MQ=255)
                                ccGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCAc                                                                                                                <  8:17803/148‑1 (MQ=255)
                                 cGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGt                                                                                                              >  1:321149/1‑149 (MQ=255)
                                       cTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTcg                                                                                                                                                >  2:160644/1‑109 (MQ=255)
                                       cTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTcg                                                                                                                                                <  1:160644/109‑1 (MQ=255)
                                            gcagcTTGGCCGACAGATCCTGCAGGGTCTGGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTg                                                                                                   <  4:138399/149‑1 (MQ=255)
                                              agcTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTc                                                                                                 <  4:65392/149‑1 (MQ=255)
                                              agcTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTc                                                                                                 <  8:92842/149‑1 (MQ=255)
                                               gcTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTc                                                                                                 <  3:95033/148‑1 (MQ=255)
                                                cTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTc                                                                                                 >  4:95033/1‑147 (MQ=255)
                                                     ccGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGcgcg                                                                                          >  5:53418/1‑149 (MQ=255)
                                                      cGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGAc                                                                                                                               <  3:178546/111‑1 (MQ=255)
                                                      cGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGAc                                                                                                                               >  4:178546/1‑111 (MQ=255)
                                                      cGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGt                                                                                         >  5:187946/1‑149 (MQ=255)
                                                       gACAGATGATGCAGTGTCTTGTGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGACAGTCGACGAAAAGCTCGCGGTTCGGCTACTTTACCAGGCTTTTTACCAGGTGCAGGCCTTTGTCCAGAGCGtt                                                                                        <  4:101371/149‑1 (MQ=255)
                                                       gACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGtt                                                                                        <  2:24083/149‑1 (MQ=255)
                                                           ataTCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTc                                                                                                                                          <  8:362203/94‑1 (MQ=255)
                                                           aGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCg                                                                                     >  4:299882/1‑149 (MQ=255)
                                                             aTCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCgg                                                                                    >  7:131098/1‑148 (MQ=255)
                                                                cTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACt                                                                                >  2:21281/1‑149 (MQ=255)
                                                                   cAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACtcgt                                                                             <  7:209099/149‑1 (MQ=255)
                                                                   cAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACtcgt                                                                             <  6:53418/149‑1 (MQ=255)
                                                                   cAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACtcgt                                                                             <  4:56108/149‑1 (MQ=255)
                                                                   cAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACtcgt                                                                             <  6:187946/149‑1 (MQ=255)
                                                                     gggTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACtcgtcg                                                                           <  8:131098/149‑1 (MQ=255)
                                                                       gTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGAt                                                                         <  8:204847/149‑1 (MQ=255)
                                                                        tCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTAACCACGTGCAGGCCATTGTCCAGCGCGTTGCTGTACTCGTCGAt                                                                         >  4:16195/1‑148 (MQ=255)
                                                                              ggTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGcc                                                                  >  3:310596/1‑149 (MQ=255)
                                                                              ggTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGc                                                                   <  3:299882/148‑1 (MQ=255)
                                                                              ggTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCAAGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCAATGTCCAGCGCGTTGCGGTACTCGTCGATCTTGTcc                                                                  >  3:305057/1‑149 (MQ=255)
                                                                               gTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCt                                                                 >  6:199087/1‑149 (MQ=255)
                                                                                 cGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCt                                                                 >  1:294840/1‑147 (MQ=255)
                                                                                    aaCGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTAc                                                                                 <  4:252391/128‑1 (MQ=255)
                                                                                    aaCGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTAc                                                                                 >  3:252391/1‑128 (MQ=255)
                                                                                           gtCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCgg                                                     <  4:66210/149‑1 (MQ=255)
                                                                                              gTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCCGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGAGGTACTCGTCGATCTTGGCCTGAGAGCTCTCGGCGt                                                  >  8:294986/1‑149 (MQ=255)
                                                                                              gTGACGCGCGGTCCAGGCATGGCCCAGGTACGGAAGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTAATCGTCGATCTTGGCCTGAGCGCGCTCGGCGt                                                  >  7:72359/1‑149 (MQ=255)
                                                                                                  cgcgcgGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGAt                                              <  8:302502/149‑1 (MQ=255)
                                                                                                   gcgcgGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACtcgt                                                                             >  4:342862/1‑117 (MQ=255)
                                                                                                   gcgcgGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATg                                             >  1:266554/1‑149 (MQ=255)
                                                                                                    cgcgGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGc                                            <  6:413002/149‑1 (MQ=255)
                                                                                                         tCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCGATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGAGCGCTCGGCGTCGATGCGACCg                                       >  5:301773/1‑149 (MQ=255)
                                                                                                            aGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGcc                                    <  5:199087/149‑1 (MQ=255)
                                                                                                                    cccAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAg                            >  1:371361/1‑149 (MQ=255)
                                                                                                                    cccAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAg                            <  2:321149/149‑1 (MQ=255)
                                                                                                                     ccAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAgc                           >  6:420018/1‑149 (MQ=255)
                                                                                                                       aGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCt                                                                 >  1:78092/1‑109 (MQ=255)
                                                                                                                       aGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCt                                                                 <  2:78092/109‑1 (MQ=255)
                                                                                                                             ggACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGc                                                                   >  4:15015/1‑101 (MQ=255)
                                                                                                                             ggACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGc                                                                   <  3:15015/101‑1 (MQ=255)
                                                                                                                             ggACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCGGc                   <  7:252609/149‑1 (MQ=255)
                                                                                                                                        gccgaACCGCGAGCTGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCGGCATACAGCt           >  1:374351/3‑146 (MQ=255)
                                                                                                                                        gccgaACCGCGAGCTGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCGGCATACAGCt           <  2:374351/144‑1 (MQ=255)
                                                                                                                                          cgaACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCGGCATACAGCTCACgg      <  2:266554/149‑1 (MQ=255)
                                                                                                                                          cgaACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCGGCATACAGCTCACgg      <  2:371361/149‑1 (MQ=255)
                                                                                                                                           gaACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCt                                   <  6:127735/119‑1 (MQ=255)
                                                                                                                                           gaACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCt                                   >  5:127735/1‑119 (MQ=255)
                                                                                                                                           gaACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCGGCATAc               <  6:377562/139‑1 (MQ=255)
                                                                                                                                           gaACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCGGCATAc               >  5:377562/1‑139 (MQ=255)
                                                                                                                                           gaACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTAGGCATACAGCTCACgg      >  3:48983/1‑148 (MQ=255)
                                                                                                                                              cAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCGGCATACAGCTCACGGGTAg  >  1:234036/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CTTCGACACCTGACGCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCGGCATACAGCTCACGGGTAG  >  NC_002947/4733883‑4734172

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: