Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,237,681 C→A 6.9% E205* (GAG→TAG) ‡ catB ← muconate cycloisomerase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,237,6810CA 6.9% 213.7 / 3.0 69E205* (GAG→TAG) ‡catBmuconate cycloisomerase 1
Reads supporting (aligned to +/- strand):  ref base C (38/26);  new base A (3/2);  total (41/28)
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.80e-01

TTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCTTGAGGTCCTGCGCCACCGGGTTGGCGCCGATTTTCAGCTTGAACACGCGGTGCCGGCGAATGTCCAGCAT  >  NC_002947/4237534‑4237825
                                                                                                                                                   |                                                                                                                                                
tttCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGt                                                                                                                                                >  1:198449/1‑150 (MQ=255)
     atcgaCTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACGAGTCCCAg                                                                                                                                           <  1:2451058/150‑1 (MQ=255)
      tcgaCTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCAATCAGGTCGATGCGATTATCGGCGAATTCCCAGCAGTTGCGGATGTCCTGGGACTCGTCCCAga                                                                                                                                          >  2:1205032/1‑149 (MQ=255)
             aTCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGtt                                                                                                                                   >  2:2252759/1‑150 (MQ=255)
             aTCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGtt                                                                                                                                   >  1:1936691/1‑150 (MQ=255)
             aTCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAAGCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACAAGGCATGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGtt                                                                                                                                   >  2:1379936/1‑150 (MQ=255)
                   cATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTc                                                                                                                             >  1:716357/1‑150 (MQ=255)
                    aTGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCg                                                                                                                            >  1:241015/1‑150 (MQ=255)
                                  aCTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGgc                                                                                                              <  2:723000/150‑1 (MQ=255)
                                  aCTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGGCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGCTTGACGTCGACCCGCACGCTGgc                                                                                                              <  2:1193934/150‑1 (MQ=255)
                                    tGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCAGAACGCTGgcgc                                                                                                            >  1:833119/1‑150 (MQ=255)
                                      cgcTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTg                                                                                                          >  1:732923/1‑150 (MQ=255)
                                       tcTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTCCTCGATCAGGTCGATGCCGTTGTCGCCCCATACCTGGCAGGCGCGGATGCTGCGGCGTTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGt                                                                                                         <  1:2274598/149‑1 (MQ=255)
                                       gcTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTGCTGGTTGACGTCGACCCGCACGCTGGCGCTGt                                                                                                         >  2:1975448/1‑150 (MQ=255)
                                        cTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTc                                                                                                        >  2:951922/1‑150 (MQ=255)
                                         tGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCCACCCGCAGACTTGCGCTGTCa                                                                                                       >  2:2274598/1‑150 (MQ=255)
                                           gTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCAcc                                                                                                     <  2:241015/150‑1 (MQ=255)
                                             tCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCa                                                                                                   <  2:1992132/150‑1 (MQ=255)
                                                 gcgcACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGCTCCGGTACCTGGTCATGTCGCCCAATAACTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTc                                                                                               >  2:1797351/1‑150 (MQ=255)
                                                 gcgcACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGCCAAGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTc                                                                                               >  1:1331436/1‑150 (MQ=255)
                                                 gcgcACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTc                                                                                               <  2:1942599/150‑1 (MQ=255)
                                                 gcgcACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACGAGTCCCAgg                                                                                                                                          >  2:488550/1‑106 (MQ=255)
                                                 gcgcACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACGAGTCCCAgg                                                                                                                                          <  1:488550/107‑2 (MQ=255)
                                                    cACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGAACCGCACGCTAGCGCTGTCACCCAGCTcgcg                                                                                            >  1:603338/1‑150 (MQ=255)
                                                     aCCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCgcgc                                                                                           <  1:1677474/150‑1 (MQ=255)
                                                     aCCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGACAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCgcgc                                                                                           <  1:2252759/150‑1 (MQ=255)
                                                       cTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCtt                                                                                         >  2:1384298/1‑150 (MQ=255)
                                                              gcgcGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTc                                                                                                                                                  >  2:589207/1‑86 (MQ=255)
                                                              gcgcGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTc                                                                                                                                                  <  1:589207/86‑1 (MQ=255)
                                                              gcgcGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATcgcg                                                                                  >  2:883438/1‑150 (MQ=255)
                                                              gcgcGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCAAAACGGACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATcgcg                                                                                  >  2:1796737/1‑150 (MQ=255)
                                                               cgcgGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCATCTCGCGCTTGATCGCGa                                                                                 <  1:618980/150‑1 (MQ=255)
                                                                gcgGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGAc                                                                                >  1:295268/1‑150 (MQ=255)
                                                                           gcgAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCTTGa                                                                     <  1:1384298/150‑1 (MQ=255)
                                                                            cgAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCTTGAg                                                                    <  2:297750/150‑1 (MQ=255)
                                                                            cgAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCTTGAg                                                                    >  2:58735/1‑150 (MQ=255)
                                                                                 tCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCTTGAGGTCCt                                                               >  1:1219436/1‑150 (MQ=255)
                                                                                         tCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGc                                                                         <  1:360560/132‑1 (MQ=255)
                                                                                         tCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGc                                                                         >  2:360560/1‑132 (MQ=255)
                                                                                         tCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCTTGAg                                                                    <  2:1165043/137‑1 (MQ=255)
                                                                                         tCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCTTGAg                                                                    >  1:1165043/1‑137 (MQ=255)
                                                                                           gATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCTTGAGGTCCTGCGCCACCgg                                                     >  1:1249490/1‑150 (MQ=255)
                                                                                           gATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGCTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCt                                                                        <  1:536244/131‑1 (MQ=255)
                                                                                           gATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGCTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCt                                                                        >  2:536244/1‑131 (MQ=255)
                                                                                                 gTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCTTGAGGTCCTGCGCCACCGGGTTGgc                                               >  2:304005/1‑150 (MQ=255)
                                                                                                 gTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCTTGAGGTCCTGCGCCACCGGGTTGgc                                               >  2:303984/1‑150 (MQ=255)
                                                                                                    gATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCCCGCGCTTGTTGAAGACCGCGTGCAGCATGTCTTGCGCCACCGTGTTGGCGcc                                            >  2:857020/1‑150 (MQ=255)
                                                                                                         cGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCTTGa                                                                     <  1:1824353/120‑1 (MQ=255)
                                                                                                         cGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCTTGa                                                                     >  2:1824353/1‑120 (MQ=255)
                                                                                                          gTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCTTGAGGTCCTGCGCCACCGGGTTGGCGCCGAtttt                                      >  1:2128889/1‑150 (MQ=255)
                                                                                                            tGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCTTGAGGTCCTGCGCCACCGGGTTGGCGCCGATTTTCa                                    <  1:304005/150‑1 (MQ=255)
                                                                                                            tGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCTTGAGGTCCTGCGCCACCGGGTTGGCGCCGATTTTCa                                    <  1:303984/150‑1 (MQ=255)
                                                                                                             gTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCTTGAGGTCCTGCGCCACCGGGTTGGCGCCGATTTTCAg                                   <  1:883438/150‑1 (MQ=255)
                                                                                                                gCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGt                                                                                                                              <  1:13408/56‑1 (MQ=255)
                                                                                                                gCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGt                                                                                                                              >  2:13408/1‑56 (MQ=255)
                                                                                                                gCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCTTGAGGTCCTGCGCCACCGGGTTGGCGCCGATTTTCAGCtt                                >  2:986033/1‑150 (MQ=255)
                                                                                                                 cccAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCTTGAGGTCCTGCGCCACCGGGTTGGCGCCGATTTTCAGCTTg                               >  1:1768870/1‑150 (MQ=255)
                                                                                                                    aaTACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCTTGAGGTCCTGCGCCACCGGGTTGGCGCCGATTTTCAGCTTGAac                            >  1:1271213/1‑150 (MQ=255)
                                                                                                                       aCCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCTTGAGGTCCTGCGCCACCGGGTTGGCGCCGATTTTCAGCTTGAACAcg                         >  1:979855/1‑150 (MQ=255)
                                                                                                                        ccTGGCAGGCGCGGATTGCATGGGACCCGCCACAATACTGGTTGACACCGACCCGCACGCTGGCGCTGTCACCCAGCCCGCGCTTGATCGCGACCACGTGCTTGAGGTCCTGCGCCACCGGGTTGGCGCCGATTTTCAGCTTGAACAcgc                        <  2:833119/150‑1 (MQ=255)
                                                                                                                              aGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCTTGAGGTCCTGCGCCACCGGGTTGGCGCCGATTTTCAGCTTGAACACGCGGTGcc                  >  2:519887/1‑150 (MQ=255)
                                                                                                                                gcgcGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCTTGAGGTCCTGCGCCACCGGGTTGGCGCCGATTTTCAGCTTGAACACGCGGTGCCgg                <  1:1797351/150‑1 (MQ=255)
                                                                                                                                gcgcGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCTTGAGGTCCTGCGCCACCGGGTTGGCGCCGATTTTCAGCTTGAACACGCGGTGCCgg                <  2:1936691/150‑1 (MQ=255)
                                                                                                                                         tacTGGGACGAGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGTCGCGGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCTTGAGGTCCTGCGCCACCGGGTTGGCGCCGATTTTCAGCTTGAACACGCGGTGCCGGCGAATGTcc       >  1:1232876/3‑150 (MQ=255)
                                                                                                                                         gCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCTTGAGGTCCTGCGCCACCGGGTTGGCGCCGAGGTGCAGCTTGAACACGCGGTGCCGGCGAATGTcc       <  1:1379936/150‑1 (MQ=255)
                                                                                                                                          ccTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCTTGAGGTCCTGCGCCACCGGGTTGGCGCCGATTTTCAGCTTGAACACGCGGTGCCGGCGAATGTCCa      <  2:1219436/150‑1 (MQ=255)
                                                                                                                                            tGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCTTGAGGTCCTGCGCCACCGGGTTGGCGCCGATTTTCAGCTTGAACACGCGGTGCCGGCGAATGTCCAGc    >  1:1321093/1‑150 (MQ=255)
                                                                                                                                              ggACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCTTGAGGTCCTGCGCCACCGGGTTGGCGCCGATTTTCAGCTTGAACACGCGGTGCCGGCGAATGTCCAGCAt  <  1:1796737/150‑1 (MQ=255)
                                                                                                                                              ggACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCTTGAGGTCCTGCGCCACCGGGTTGGCGCCGATTTTCAGCTTGAACACGCGGTGCCGGCGAATGTCCAGCAt  >  2:146200/1‑150 (MQ=255)
                                                                                                                                                tcTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCTTGAGGTCCTGCGCCACCGGGTTGGCGCCGATTTTCAGCTTGAACACGCGGTGCCGGCGAAt           <  2:14804/138‑1 (MQ=255)
                                                                                                                                                gcTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCTTGAGGTCCTGCGCCACCGGGTTGGCGCCGATTTTCAGCTTGAACACGCGGTGCCGGCGAAt           >  1:14804/2‑139 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGACCACGTGCTTGAGGTCCTGCGCCACCGGGTTGGCGCCGATTTTCAGCTTGAACACGCGGTGCCGGCGAATGTCCAGCAT  >  NC_002947/4237534‑4237825

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