Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,005,056 G→C 12.8% intergenic (+38/‑138) PP_5259 → / → ydcJ LysR family transcriptional regulator/metalloprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,005,0560GC12.8% 193.7 / 16.9 70intergenic (+38/‑138)PP_5259/ydcJLysR family transcriptional regulator/metalloprotein
Reads supporting (aligned to +/- strand):  ref base G (28/33);  new base C (4/5);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

CATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGCCAA  >  NC_002947/6004908‑6005204
                                                                                                                                                    |                                                                                                                                                    
cATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACcgc                                                                                                                                                     >  2:947097/1‑150 (MQ=255)
  tGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTAGGGGGACCGCGt                                                                                                                                                   >  2:1168530/1‑150 (MQ=255)
  tGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGt                                                                                                                                                   >  1:1065004/1‑150 (MQ=255)
  tGCAATGCGCAGTCCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGCGGACCGCGt                                                                                                                                                   <  2:1152457/150‑1 (MQ=255)
    cAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAg                                                                                                                                                 <  2:964829/150‑1 (MQ=255)
    cAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAg                                                                                                                                                 >  1:1188037/1‑150 (MQ=255)
    cAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAg                                                                                                                                                 >  1:1187809/1‑150 (MQ=255)
    cAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCCCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAg                                                                                                                                                 <  2:297990/150‑1 (MQ=39)
    cAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCCTTTCCCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAg                                                                                                                                                 <  1:1053070/150‑1 (MQ=37)
     aaTGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAgg                                                                                                                                                >  2:12352/1‑150 (MQ=255)
          gcAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCGAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGAGCgg                                                                                                                                           >  1:388225/1‑150 (MQ=38)
          gcAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGAGCgg                                                                                                                                           <  1:1000359/150‑1 (MQ=38)
          gcAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGAGCgg                                                                                                                                           <  2:388225/150‑1 (MQ=38)
            aGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCCCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCt                                                                                                                                         <  1:12352/150‑1 (MQ=255)
                 cGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGc                                                                                                                                    <  1:391256/150‑1 (MQ=255)
                   gAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCg                                                                                                                                  >  2:12358/1‑150 (MQ=255)
                    aGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCgg                                                                                                                                 <  2:395633/150‑1 (MQ=255)
                        cATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCGAAATGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGaaa                                                                                                                             >  2:1102921/1‑150 (MQ=255)
                         aTTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAg                                                                                                                            >  2:1212569/1‑150 (MQ=255)
                         aTTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAg                                                                                                                            >  1:192557/1‑150 (MQ=255)
                          ttACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAgg                                                                                                                           <  2:545140/150‑1 (MQ=255)
                           tACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAggg                                                                                                                          <  2:1187809/150‑1 (MQ=255)
                           tACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAggg                                                                                                                          <  2:1188037/150‑1 (MQ=255)
                              cTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAgg                                                                                                                           <  1:515810/146‑1 (MQ=255)
                              cTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAgg                                                                                                                           >  2:515810/1‑146 (MQ=255)
                                   gTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGTGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCt                                                                                                                                         >  2:1291145/1‑127 (MQ=255)
                                   gTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCt                                                                                                                                         <  1:1291145/127‑1 (MQ=255)
                                   gTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAg                                                                                                                  >  2:1239921/1‑150 (MQ=255)
                                         cGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTTCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGAt                                                                                                            <  2:769450/150‑1 (MQ=255)
                                         cGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAg                                                                                                                  <  2:1417780/144‑1 (MQ=255)
                                         cGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAg                                                                                                                  >  1:1417780/1‑144 (MQ=255)
                                         cGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGAt                                                                                                            <  1:129525/150‑1 (MQ=255)
                                          gAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATc                                                                                                           <  1:996775/150‑1 (MQ=255)
                                             cATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCaaa                                                                                                        >  1:1166878/1‑150 (MQ=255)
                                                gCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCaaacaa                                                                                                     <  2:1065004/150‑1 (MQ=255)
                                                gCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCaaaaaa                                                                                                     <  1:947097/150‑1 (MQ=255)
                                                 cGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAg                                                                                                    <  2:1166878/150‑1 (MQ=255)
                                                  ggCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAgg                                                                                                   >  2:29205/1‑150 (MQ=255)
                                                  ggCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCAGTCCGACTGGGACGGCGTGGGCGGGGCTTTAGCCGGGAAGGGGCAGTCGAGGTCGATCAAAAAAgg                                                                                                   >  2:1229627/1‑150 (MQ=255)
                                                         gTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGaca                                                                                            <  1:29205/150‑1 (MQ=255)
                                                           gCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGacaca                                                                                          <  1:954935/150‑1 (MQ=255)
                                                              aaaGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGCGGTCCCCGTAgga                                                                                                                                               <  1:1183297/94‑2 (MQ=38)
                                                              aaaGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGCGGTCCCCGTAgga                                                                                                                                               >  2:1183297/1‑93 (MQ=38)
                                                                aGGTGCGGCCGCTGCTGGAATGGATTGGCGACCAGTTTCAGCTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACTGCGACCGCGTACGAGCGTTTTTAGCCGGGAAAGGGCCGTAACGAGCTATCAAAAAAGGAATGACAcactcac                                                                                     <  1:398816/150‑1 (MQ=255)
                                                                  gTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGTCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACtt                                                                                   <  1:405721/150‑1 (MQ=255)
                                                                          gcTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGc                                                                           <  2:1130744/150‑1 (MQ=255)
                                                                           cTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCt                                                                          >  1:368968/1‑150 (MQ=255)
                                                                             ggtggAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAATCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGCAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTtg                                                                        <  1:1229627/150‑1 (MQ=255)
                                                                               tggAATAGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGTGGCTAAAGCCGCTCCTACGGGGACCGGGTAGGCGCGGGTTTATCCGGGAAAGGGGCGTCGGGGGGGATCAAAAAAGGAATGAGATAATAACTTTTTTTTGTTtgtg                                                                      >  2:573662/1‑150 (MQ=255)
                                                                                     ggATTGGCAGGCAGGTGCAGGTTTGACGGCCTCTTCGCGGCTAAAACCGTACTAACGGGCGCCGTGTGGGCCGGGCTTCTCCAGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCaaa                                                                <  1:1168530/150‑1 (MQ=255)
                                                                                           gCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCAt                                                          <  1:1212569/150‑1 (MQ=255)
                                                                                            cGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCAtt                                                         >  1:539945/1‑150 (MQ=255)
                                                                                               gCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTccc                                                      <  1:1102921/150‑1 (MQ=255)
                                                                                                     gCACGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCTGCTTTAGCCGGGAAAGCGCCGTCGAGGTCGATCAAAAAAGGAATTACACACTCACTTTTTTGCGCTTGTGAGCAAAGTGCATTCCCAGCttt                                                <  1:573662/150‑1 (MQ=255)
                                                                                                         gTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGaca                                                                                            >  1:692500/1‑102 (MQ=255)
                                                                                                         gTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGaca                                                                                            <  2:692500/102‑1 (MQ=255)
                                                                                                         gTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATg                                            >  1:1218328/1‑150 (MQ=255)
                                                                                                           tcgACGGCCCTTTCCCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGCGAAGAGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTa                                          >  2:1325805/3‑150 (MQ=17)
                                                                                                            tGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTaa                                         <  2:1157309/150‑1 (MQ=255)
                                                                                                            tGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTAAGTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTaa                                         >  1:726417/1‑150 (MQ=255)
                                                                                                                  ccTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTcc                                   >  1:191791/1‑150 (MQ=255)
                                                                                                                        cgcgGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAAccc                             >  1:860310/1‑150 (MQ=255)
                                                                                                                         ccgGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGGGAAAGGGCGGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTGGCTTTTGAGCAAAGTACATTCCCAGCTTTCATGTTAGCGTCCGAACCCa                            >  1:276754/2‑150 (MQ=37)
                                                                                                                          cgGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCaa                           <  1:1363913/150‑1 (MQ=255)
                                                                                                                              tAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAg                       >  1:51577/1‑150 (MQ=255)
                                                                                                                                 aGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGaa                    <  2:191791/150‑1 (MQ=255)
                                                                                                                                   ccGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGc                  <  1:48468/150‑1 (MQ=255)
                                                                                                                                     gCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTa                >  2:1092200/1‑150 (MQ=255)
                                                                                                                                     gCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTa                >  1:973618/1‑150 (MQ=255)
                                                                                                                                          tACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGaa                    >  2:1154350/1‑141 (MQ=255)
                                                                                                                                          tACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGaa                    <  1:1154350/141‑1 (MQ=255)
                                                                                                                                                   cgcgTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATTTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGCCaa  >  1:1015551/1‑150 (MQ=255)
                                                                                                                                                   cgcgTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGCCaa  <  1:12358/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGCCAA  >  NC_002947/6004908‑6005204

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