Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,559,745 G→C 7.2% A251G (GCT→GGT) ‡ hisZ ← ATP phosphoribosyltransferase regulatory subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,559,7450GC 7.2% 224.2 / 5.6 69A251G (GCT→GGT) ‡hisZATP phosphoribosyltransferase regulatory subunit
Reads supporting (aligned to +/- strand):  ref base G (34/30);  new base C (2/3);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTCGGCAGGTCGGCCGT  >  NC_002947/5559603‑5559893
                                                                                                                                              |                                                                                                                                                    
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                                                                                                                      <  4:287170/147‑1 (MQ=255)
                            ccTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGc                                                                                                                      >  3:287170/1‑147 (MQ=255)
                             cTTCGCCGGGCACGAACACGGCGAACACCACCCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACCGCGGCACGTGCGGGTATCGCGACCCCACGCGCTCCGCGCTGGCCCGCAGGTCATACACCGCCGCCCGCACg                                                                                                                   >  8:302553/1‑149 (MQ=255)
                                 gCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTg                                                                                                                <  6:62296/148‑1 (MQ=255)
                                  ccTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGc                                                                                                                  <  2:231212/145‑1 (MQ=255)
                                  ccTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGc                                                                                                                  >  1:231212/1‑145 (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                                                                                                          <  4:20251/144‑1 (MQ=255)
                                           aaCACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCggg                                                                                                          >  3:20251/1‑144 (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:282677/149‑1 (MQ=255)
                                                  cGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAgg                                                                                              <  3:164180/149‑1 (MQ=255)
                                                  cGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAgg                                                                                              <  4:208590/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                                                               <  7:302553/149‑1 (MQ=255)
                                                                                 aCGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGac                                                               >  1:119359/1‑149 (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)
                                                                                                                                              gCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCTCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACACCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTCGGCAGGTCGGCCGt  >  2:64715/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GTAGCGACCGCCCTGGGCGATCGATTGGCCTTCGCCTGGCACGAACACGGCGAACACCACGCCCGTGTGGTAGTGGTAACCACGCAGCTCGCCGAGGTCGAAGTACAGCGGCAGGTCCGGGTAGCGCGACGCCAGGCGATCAGCGATCGCCAGCAGGTCATCCAGCGCCGCCAGCACGCTGGCCGGGGCGCGACCGAGGCGCACGCGGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTCGGCAGGTCGGCCGT  >  NC_002947/5559603‑5559893

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: