Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,486,478 A→G 6.1% Y413C (TAC→TGC)  purD → phosphoribosylamine‑‑glycine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,486,4780AG 6.1% 236.5 / 3.0 82Y413C (TAC→TGC) purDphosphoribosylamine‑‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base A (34/41);  new base G (2/3);  total (38/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCCATGCCGGTACTGCGCTCAAGGATGGCAAGGTAACCACCAACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGCTCGTCGACAGGGCCCGGAATGGGCCTTGCCTTCGACTTGGCGCGC  >  NC_002947/5486330‑5486625
                                                                                                                                                    |                                                                                                                                                   
ttCAATGCCGGTACTGCGCTCAAGGATGGCAAGGTAACCACCAACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTa                                                                                                                                                     >  2:192505/1‑149 (MQ=255)
      gCCGGTACTGCGCTCAAGGATGGCAAGGTAACCACCAACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTa                                                                                                                                                     >  5:238332/1‑143 (MQ=255)
      gCCGGTACTGCGCTCAAGGATGGCAAGGTAACCACCAACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTa                                                                                                                                                     <  6:238332/143‑1 (MQ=255)
       cctgTACTGCGCTCAAGGATGGCAAGGTAACCACCAACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCTGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTcc                                                                                                                                              <  7:188050/146‑1 (MQ=255)
       ccGGTACTGCGCTCAAGGATGGCAAGGTAACCACCAACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTc                                                                                                                                               <  4:112750/148‑1 (MQ=255)
       ccGGTACTGCGCTCAAGGATGGCAAGGTAACCACCAACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTc                                                                                                                                               <  6:230766/148‑1 (MQ=255)
        cGGTACTGCGCTCAAGGATGGCAAGGTAACCACCAACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGTAGAAGCTGCCGTTcca                                                                                                                                             <  8:70588/149‑2 (MQ=255)
         ggTACTGCGCTCAAGGATGGCAAGGTAACCACCAACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTAc                                                                                                                                                    >  8:63284/1‑141 (MQ=255)
         ggTACTGCGCTCAAGGATGGCAAGGTAACCACCAACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTAc                                                                                                                                                    <  7:63284/141‑1 (MQ=255)
         ggTACGGCTCTCCAGGATGTCAACGTAACCACCAACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGTAGAAGCTGCCGTTccag                                                                                                                                            <  8:289402/149‑3 (MQ=255)
              tGCGCTCAAGGATGGCAAGGTAACCACCAACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCg                                                                                                                                       >  3:214135/1‑149 (MQ=255)
                cgcTCAAGGATGGCAAGGTAACCACCAACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCgg                                                                                                                                      <  7:221529/148‑1 (MQ=255)
                 gcTCAAGGATGGCAAGGTAACCACCAACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGc                                                                                                                                     >  3:98932/1‑148 (MQ=255)
                  cTCAAGGATGGCAAGGTAACCACCAACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCgg                                                                                                                                      >  2:29846/1‑146 (MQ=255)
                  cTCAAGGATGGCAAGGTAACCACCAACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCgg                                                                                                                                      <  1:29846/146‑1 (MQ=255)
                   tCAAGGATGGCAAGGTAACCACCAACGGCGGGCGCGTGCTGTGTGCCCCCGCCATGGGCAGCACCGTTGCCGACGCGCAGCCGCAGGCGTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCGACCGTTCCGGCCTCGGCTAt                                                                                                                                  >  8:28067/1‑149 (MQ=255)
                   tCAAGGATGGCAAGGTAACCACCAACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTa                                                                                                                                                     >  6:164265/1‑130 (MQ=255)
                   tCAAGGATGGCAAGGTAACCACCAACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTa                                                                                                                                                     <  5:164265/130‑1 (MQ=255)
                   tCAAGGATGGCAAGGTAACCACCAACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGa                                                                                                                                            <  4:245537/139‑1 (MQ=255)
                   tCAAGGATGGCAAGGTAACCACCAACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGa                                                                                                                                            >  3:245537/1‑139 (MQ=255)
                   tCAAGGATGGCAAGGTAACCACCAACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCCGGAACGGTAGAAGCTGCCGTTccagctgac                                                                                                                                       <  4:265621/144‑8 (MQ=255)
                   tCAAGGATGGCAAGGTAACCACCAACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCCGGAACGGTAGAAGCTGCCGTTccagctgac                                                                                                                                       >  3:265621/1‑137 (MQ=255)
                      aGGATGTCAAGGTAACCACCAACGGCGGGCGCGGGCTGTGTGCCACCGCCATTGGCAGCACCGTTGCCGACGCGCATCAGCATGGTTATCGCCTGGCCAATGAAGTCAGCTGGATCGGCAGGTTCTACCGTGCCTACATCGGCTATCgg                                                                                                                               <  8:274984/149‑1 (MQ=255)
                      aGGATGGCAAGGTAACCACCAACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCgg                                                                                                                               <  4:214135/149‑1 (MQ=255)
                       gaaTGGCAAGGTAACCACCAACGGCGGGCGCGTGCTGTGTGACACCGCCATGGGCAGCACCGTTTCCGACGCGCAGCCTCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTCCCGTTCAGACATCGGCTATCggg                                                                                                                              >  7:199493/3‑149 (MQ=255)
                            gCAAGGTAACCACCTAGGGCGGGCGCGTGCTGTCTGTCACCGCCATCGGCAGCACCGTTGCCGACGCGAAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACCGCAGCTTCTACCGTTCCGACATCGGCTACCCGGCCATc                                                                                                                         <  8:199493/149‑1 (MQ=255)
                              aaGGTAACCACCAACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACAGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTGCCGTGCAGACATCGGCTCTCGGGCCATCGc                                                                                                                       >  4:21027/1‑149 (MQ=255)
                                ggTAACCACCAACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGccc                                                                                                                     <  4:98932/149‑1 (MQ=255)
                                  tAACCACCAACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCcgc                                                                                                                   <  6:277619/149‑1 (MQ=255)
                                   aaCCACCAACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGAAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCcgc                                                                                                                   >  5:134735/1‑148 (MQ=255)
                                    accaccAACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGa                                                                                                                 <  1:192505/149‑1 (MQ=255)
                                       accaACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAgcg                                                                                                              >  4:8966/1‑149 (MQ=255)
                                        ccaACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGccc                                                                                                                     >  7:234280/1‑141 (MQ=255)
                                        ccaACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGccc                                                                                                                     <  8:234280/141‑1 (MQ=255)
                                        ccaACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAgcgc                                                                                                             <  4:130965/149‑1 (MQ=255)
                                        ccaACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAgcgc                                                                                                             >  3:130965/1‑149 (MQ=255)
                                          aaCGGCGGGCGCGTGCTGTGTGCCACCGCTATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCgg                                                                                                           <  3:21027/149‑1 (MQ=255)
                                          aaCGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCcg                                                                                                                    <  6:111323/140‑1 (MQ=255)
                                          aaCGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCcg                                                                                                                    >  5:111323/1‑140 (MQ=255)
                                           aCGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGt                                                                                                          <  6:42772/149‑1 (MQ=255)
                                           aCGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGt                                                                                                          <  7:28067/149‑1 (MQ=255)
                                              gcgggcgCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAg                                                                                                       <  5:152809/149‑1 (MQ=255)
                                                 ggcgCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCAc                                                                                                    >  2:267568/1‑149 (MQ=255)
                                                 ggcgCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCAc                                                                                                    <  1:267568/149‑1 (MQ=255)
                                                    gcgTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCATTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACcag                                                                                                 >  1:81140/1‑149 (MQ=255)
                                                      gTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACcag                                                                                                 >  5:4031/1‑147 (MQ=255)
                                                        gCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCagt                                                                                             <  4:61468/149‑1 (MQ=255)
                                                         cTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCAc                                                                                                    >  8:13256/1‑141 (MQ=255)
                                                         cTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCAc                                                                                                    <  7:13256/141‑1 (MQ=255)
                                                           gtgtgCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCagtaag                                                                                          >  2:2970/1‑149 (MQ=255)
                                                            tgtgCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCagtaa                                                                                           <  5:19423/147‑1 (MQ=255)
                                                             gtgCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCagtaagta                                                                                        >  1:178216/1‑149 (MQ=255)
                                                              tgCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCagtaagta                                                                                        <  3:189950/148‑1 (MQ=255)
                                                                ccACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCCGCTGGACCGGCAGCTTCTACCGTACCCACATCGGCTATCGGGCCAACGCCCGCGAGCGCGGTCAGCACCCGCAGAAAGaactc                                                                                     >  8:45429/1‑148 (MQ=255)
                                                                 cACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGtactac                                                                                    <  1:2970/149‑1 (MQ=255)
                                                                      ccATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTa                                                                               >  1:255032/1‑149 (MQ=255)
                                                                         tGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGc                                                                            <  7:209819/149‑1 (MQ=255)
                                                                            gcagcaCCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCg                                                                          >  2:193232/1‑148 (MQ=255)
                                                                                     ttGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTc                                                                <  2:255032/149‑1 (MQ=255)
                                                                                        ccGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGc                                                             <  3:868/149‑1 (MQ=255)
                                                                                             gcgcAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCg                                                        <  3:49841/149‑1 (MQ=255)
                                                                                             gcgcAGCAACAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCg                                                        >  2:76537/1‑149 (MQ=255)
                                                                                               gcagcagcaGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCACCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCgg                                                       >  6:169410/1‑148 (MQ=255)
                                                                                                cagcagcagGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTCCCGTTCCGCCATCGGCTATCGCGCCATCGCCCGCTCGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTGTAAGCCTGCGGCa                                                     >  5:232802/1‑149 (MQ=255)
                                                                                                   cagcagGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATcc                                                  >  8:69160/1‑149 (MQ=255)
                                                                                                   cagcagGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATcc                                                  <  6:134735/149‑1 (MQ=255)
                                                                                                    agcagGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCg                                                 >  2:151238/1‑149 (MQ=255)
                                                                                                      gagGGTTTTCGGCTGGCCAAGGAAGTCTGCGGGCACGGCAGCCTCTACCGTTCCGTCATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGc                                               <  7:45429/148‑1 (MQ=255)
                                                                                                      cagGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGttt                                                                 <  7:154144/131‑1 (MQ=255)
                                                                                                      cagGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGttt                                                                 >  8:154144/1‑131 (MQ=255)
                                                                                                      cagGCTTATCGCCTGGCCAAGGAAGGCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCATCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGc                                               <  1:79751/149‑1 (MQ=255)
                                                                                                       agGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGCt                                              <  8:95075/149‑1 (MQ=255)
                                                                                                           ttATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGCtcgtc                                          >  4:34326/1‑149 (MQ=255)
                                                                                                           ttATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGCtcgtc                                          >  6:150189/1‑149 (MQ=255)
                                                                                                            tATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGCtcgtc                                          >  5:24182/1‑148 (MQ=255)
                                                                                                                gCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGCTCGTCGACAg                                     >  2:56078/1‑149 (MQ=255)
                                                                                                                      ccAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGCTCGTCGACAGGGCCCg                               <  1:76537/149‑1 (MQ=255)
                                                                                                                        aaGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGCTCGTCGACAGGGCCCgg                              >  1:50816/1‑148 (MQ=255)
                                                                                                                           gAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGCTCGTCGACAGGGCCCGGAATg                          >  5:197649/1‑149 (MQ=255)
                                                                                                                            aaGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGCTCGTCGACAGGGCCCGGAATgg                         <  2:178216/149‑1 (MQ=255)
                                                                                                                             aGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGCTCGTCGACAGGGCCCGGAATggg                        <  2:50816/149‑1 (MQ=255)
                                                                                                                              gTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGCTCGTCGACAGGGCCCGGAATGGgc                       <  5:150189/149‑1 (MQ=255)
                                                                                                                                 aGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGCTCGTCGACAGGGCCCGGAATGGgcctt                    >  4:275697/1‑149 (MQ=255)
                                                                                                                                    tGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGCTCGTCGACAGGGCCCGGAATGGgccttgcc                 >  2:264056/1‑149 (MQ=255)
                                                                                                                                             aGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGCTCGTCGACAGGGCCCGGAATGGGCCTTGCCTTCGACTTg        <  3:34326/149‑1 (MQ=255)
                                                                                                                                                  cTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGCTCGTCGACAGGGCCCGGAATGGGCCTTGCCTTCGACTTGgcgc    <  3:8966/148‑1 (MQ=255)
                                                                                                                                                   tACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGCTCGTCGACAGGGCCCGGAATGGGCCTTGCCTTCGACTTGgcgcgc  >  1:37893/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTCCATGCCGGTACTGCGCTCAAGGATGGCAAGGTAACCACCAACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGCTCGTCGACAGGGCCCGGAATGGGCCTTGCCTTCGACTTGGCGCGC  >  NC_002947/5486330‑5486625

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: