Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,559,744 A→G 7.2% A251A (GCT→GCC) ‡ hisZ ← ATP phosphoribosyltransferase regulatory subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,559,7440AG 7.2% 193.4 / 3.7 69A251A (GCT→GCC) ‡hisZATP phosphoribosyltransferase regulatory subunit
Reads supporting (aligned to +/- strand):  ref base A (33/30);  new base G (2/3);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTCGGCAG  >  NC_002947/5559596‑5559884
                                                                                                                                                    |                                                                                                                                            
tatCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCa                                                                                                                                              <  6:111717/149‑1 (MQ=255)
       gtAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATc                                                                                                                                        <  7:118306/148‑1 (MQ=255)
       caatcGATCGCCCAGGGCGGTCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATGAGCGATCg                                                                                                                                       >  4:140894/5‑149 (MQ=255)
        tAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGGCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGc                                                                                                                                      <  7:323347/149‑1 (MQ=255)
         atcGATCGCCCAGGGCGGTCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGc                                                                                                                                      >  2:187451/3‑148 (MQ=255)
         aGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGc                                                                                                                                      >  8:16792/1‑148 (MQ=255)
          gCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCGCTGATCGcct                                                                                                                                    <  1:205010/149‑2 (MQ=255)
          gCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCca                                                                                                                                    >  4:282677/1‑149 (MQ=255)
           cGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCGCTGATCGcctg                                                                                                                                   <  1:65509/149‑3 (MQ=255)
                gCCCAGGGCGGTCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAgg                                                                                                                               <  1:187451/148‑1 (MQ=255)
                 cccTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGt                                                                                                                              >  4:120223/1‑148 (MQ=255)
                  ccTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCa                                                                                                                            >  5:62296/1‑149 (MQ=255)
                         gatcgatTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAgcg                                                                                                                     >  5:321495/1‑149 (MQ=255)
                           tcgatTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCGGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGcgcc                                                                                                                   <  1:228536/149‑1 (MQ=255)
                              atTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCgccgcc                                                                                                                >  6:301035/1‑149 (MQ=255)
                               ttGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGcgccgc                                                                                                                 >  1:214095/1‑147 (MQ=255)
                                  gCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATGAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCa                                                                                                            <  3:140894/149‑1 (MQ=255)
                                   ccTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGc                                                                                                             >  3:287170/1‑147 (MQ=255)
                                   ccTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGc                                                                                                             <  4:287170/147‑1 (MQ=255)
                                    cTTCGCCGGGCACGAACACGGCGAACACCACCCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACCGCGGCACGTGCGGGTATCGCGACCCCACGCGCTCCGCGCTGGCCCGCAGGTCATACACCGCCGCCCGCACg                                                                                                          >  8:302553/1‑149 (MQ=255)
                                        gCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTg                                                                                                       <  6:62296/148‑1 (MQ=255)
                                         ccTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGc                                                                                                         >  1:231212/1‑145 (MQ=255)
                                         ccTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGc                                                                                                         <  2:231212/145‑1 (MQ=255)
                                         ccTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGc                                                                                                     >  4:232013/1‑149 (MQ=255)
                                             gCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCggg                                                                                                 <  2:18888/149‑1 (MQ=255)
                                             gCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCg                                                                                                   <  3:232013/147‑1 (MQ=255)
                                               aCGAGCACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGgc                                                                                               >  3:208590/1‑149 (MQ=255)
                                               aCGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGgc                                                                                               >  3:139934/1‑149 (MQ=255)
                                               aCGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGgc                                                                                               <  5:301035/149‑1 (MQ=255)
                                                 gAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGgcgc                                                                                             >  2:303595/1‑149 (MQ=255)
                                                 gAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGgcgc                                                                                             >  6:141822/1‑149 (MQ=255)
                                                  aaCACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCggg                                                                                                 >  3:20251/1‑144 (MQ=255)
                                                  aaCACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCggg                                                                                                 <  4:20251/144‑1 (MQ=255)
                                                       ggCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGAcc                                                                                         >  3:21122/1‑147 (MQ=255)
                                                       ggCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGa                                                                                       <  3:120223/149‑1 (MQ=255)
                                                        gCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATc                                                                                                                          >  6:108650/1‑113 (MQ=255)
                                                        gCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATc                                                                                                                          <  5:108650/113‑1 (MQ=255)
                                                        gCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAg                                                                                      >  2:176644/1‑149 (MQ=255)
                                                         cGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAgg                                                                                     <  3:164180/149‑1 (MQ=255)
                                                         cGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAgg                                                                                     <  4:208590/149‑1 (MQ=255)
                                                         cGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAgg                                                                                     <  3:282677/149‑1 (MQ=255)
                                                          gAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGgc                                                                                    <  1:303595/149‑1 (MQ=255)
                                                            acacCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGgcgc                                                                                  <  4:21122/149‑1 (MQ=255)
                                                                cacGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCAcgc                                                                              >  2:7863/1‑149 (MQ=255)
                                                                        tgtgGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCa                                                                      <  1:7863/149‑1 (MQ=255)
                                                                            gtagtggtaACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGc                                                                    >  1:17060/1‑147 (MQ=255)
                                                                             tagtggtaACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAg                                                                 <  5:141822/149‑1 (MQ=255)
                                                                              agtggtaACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAg                                                                 >  4:153354/1‑148 (MQ=255)
                                                                                 ggtaACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACt                                                             >  6:368993/1‑149 (MQ=255)
                                                                                        aCGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGac                                                      >  1:119359/1‑149 (MQ=255)
                                                                                        aCGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGac                                                      <  7:302553/149‑1 (MQ=255)
                                                                                           cAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGaccac                                                   <  1:176644/149‑1 (MQ=255)
                                                                                            aGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCa                                                                                                               >  8:127943/1‑88 (MQ=255)
                                                                                            aGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCa                                                                                                               <  7:127943/88‑1 (MQ=255)
                                                                                            aGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGaccacc                                                  >  6:348772/1‑149 (MQ=255)
                                                                                                 gCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACcacac                                              <  4:139934/148‑1 (MQ=255)
                                                                                                   cGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGc                                           >  5:369359/1‑149 (MQ=255)
                                                                                                     aGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTc                                         <  2:214095/149‑1 (MQ=255)
                                                                                                         cGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGAc                                      >  6:72714/1‑148 (MQ=255)
                                                                                                                   ctgCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCgca                           <  8:165285/147‑1 (MQ=255)
                                                                                                                               gTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCa               <  6:369359/149‑1 (MQ=255)
                                                                                                                                   cgcgACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGgcgc                              <  2:269635/130‑1 (MQ=255)
                                                                                                                                   cgcgACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGgcgc                              >  1:269635/1‑130 (MQ=255)
                                                                                                                                   cgcgACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTc           >  3:164296/1‑149 (MQ=255)
                                                                                                                                    gcgACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCt          >  2:326516/1‑149 (MQ=255)
                                                                                                                                       aCGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTc       <  7:16792/149‑1 (MQ=255)
                                                                                                                                         tgctggCGATCGCTGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAgg                                  >  8:341024/5‑120 (MQ=255)
                                                                                                                                         tgctggCGATCGCTGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAgg                                  <  7:341024/116‑1 (MQ=255)
                                                                                                                                           cAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTCgg     >  7:160359/1‑147 (MQ=255)
                                                                                                                                            tggCGATCGCTGATCGCCAGCAGGTCATCCAGCGACGCCAGCACGCTCGCCGCGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTCCCCAGGTCATACGGCAg  >  4:90810/2‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
TATCGTCGTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTCGGCAG  >  NC_002947/5559596‑5559884

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