Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,734,025 A→C 6.9% 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 6.9% 212.9 / 6.6 73L537R (CTG→CGG) sucA2‑oxoglutarate dehydrogenase subunit E1
Reads supporting (aligned to +/- strand):  ref base A (33/35);  new base C (2/3);  total (35/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 = 8.18e-01

AGATCTTCGACACCTGACGCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCGGCATACAGCTCACGGGTA  >  NC_002947/4733879‑4734171
                                                                                                                                                  |                                                                                                                                                  
aGATCTTCGACACCTGACGCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAg                                                                                                                                                   >  4:34132/1‑148 (MQ=255)
aGATCTTCGACACCTGACGCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGc                                                                                                                                                  >  4:38813/1‑149 (MQ=255)
    cTTCGACACCTGACGCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTcgc                                                                                                                                              >  4:303771/1‑149 (MQ=255)
             cTGACGCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCt                                                                                                                                     >  4:150230/1‑149 (MQ=255)
               gACGCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTcc                                                                                                                                   >  3:334422/1‑149 (MQ=255)
                aCGCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACCGCGAGCTGTTCgtcga                                                                                                                                    <  4:186408/147‑5 (MQ=255)
                aCGCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTcc                                                                                                                                   <  3:232652/148‑1 (MQ=255)
                aCGCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTcc                                                                                                                                   <  5:345248/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                                                                                                                                      <  8:201170/131‑3 (MQ=255)
                              aaCCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACCGCGAGCTGTTCgtc                                                                                                                                      >  7:201170/1‑129 (MQ=255)
                                ccTTCCGGCAGCTCGAGCAGCTTGGCCGACAGTTCGTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGATGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTAc                                                                                                                  <  8:192188/149‑1 (MQ=255)
                                 cTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTAcc                                                                                                                 <  3:34132/149‑1 (MQ=255)
                                 cTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTAcc                                                                                                                 <  3:38813/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                                                                                                                <  6:319961/147‑1 (MQ=255)
                                    ccGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCa                                                                                                                <  4:334422/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                                                                            <  6:233359/149‑1 (MQ=255)
                                                                      cAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACtcgt                                                                            >  5:328339/1‑149 (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:295426/1‑149 (MQ=255)
                                                                                                 gTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGt                                                 <  2:27279/149‑1 (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)
                                                                                                                                                  |                                                                                                                                                  
AGATCTTCGACACCTGACGCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCGGCATACAGCTCACGGGTA  >  NC_002947/4733879‑4734171

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: