Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,004,997 A→T 6.4% intergenic (+35/‑141) 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,004,9970AT 6.4% 233.0 / 5.0 92intergenic (+35/‑141)PP_5259/ydcJLysR family transcriptional regulator/metalloprotein
Reads supporting (aligned to +/- strand):  ref base A (39/46);  new base T (4/2);  total (44/48)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.75e-01

CCATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCG  >  NC_002947/6004851‑6005144
                                                                                                                                                   |                                                                                                                                                    
ccATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGAcc                                                                                                                                                    >  5:99140/1‑149 (MQ=255)
 cATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACcg                                                                                                                                                   >  8:290743/1‑149 (MQ=255)
 cATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCACGGCTAAAGCCGCTCCTACGGGGACcg                                                                                                                                                   >  8:201021/1‑149 (MQ=255)
  aTGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGa                                                                                                                                                      >  1:233663/1‑145 (MQ=255)
  aTGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGa                                                                                                                                                      <  2:233663/145‑1 (MQ=255)
  aTGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACcgc                                                                                                                                                  >  5:295349/1‑149 (MQ=255)
   tGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACcgcg                                                                                                                                                 >  8:197216/1‑149 (MQ=255)
   tGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACcgcg                                                                                                                                                 >  6:209890/1‑149 (MQ=255)
     cAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATCAGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGt                                                                                                                                                <  1:298709/149‑1 (MQ=255)
          cgcAGTGCCGGAGCGCATTACCTGGCGGATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGgcgc                                                                                                                                          <  6:54971/149‑1 (MQ=255)
           gcAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGtgcg                                                                                                                                         >  7:353990/1‑149 (MQ=255)
             aGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGc                                                                                                                                       <  5:138626/149‑1 (MQ=255)
              gTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTGCAACGGGGGCCGCGTAGGCGCGGCt                                                                                                                                     >  4:20660/1‑149 (MQ=255)
              gTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGcg                                                                                                                                      >  4:168877/1‑148 (MQ=255)
               tGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCtt                                                                                                                                    <  4:358924/149‑1 (MQ=255)
                  cGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAg                                                                                                                                 <  4:313230/149‑1 (MQ=255)
                  cGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAg                                                                                                                                 <  1:147916/149‑1 (MQ=255)
                    gAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGcc                                                                                                                               <  1:222548/149‑1 (MQ=255)
                    gAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGAGCGGCTTTAGcc                                                                                                                               >  6:171339/1‑149 (MQ=255)
                     aGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGGTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGAGCGCGTAGGCGCGGCTTTAGCCg                                                                                                                              <  6:99140/149‑1 (MQ=255)
                      gcgcATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCgg                                                                                                                             >  4:335128/1‑149 (MQ=255)
                      gcgcATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCgg                                                                                                                             <  2:91784/149‑1 (MQ=255)
                      gagcATTACTTGGAGTATGCCGAGGATGCGGCAGAGGTGACAAAGGTGCGTGCGCTGGTGGAATTGATTGGCGAGCAGTTGCAGGTTTGACGGACTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCg                                                                                                                              <  8:317350/146‑1 (MQ=255)
                         ccttACCTTGCGTATGCCGTTCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGgcgc                                                                                                                                          <  7:228205/132‑1 (MQ=255)
                          aTTACCTGGTGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGaaa                                                                                                                         <  8:326295/149‑1 (MQ=255)
                          aTTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGaaa                                                                                                                         <  5:348659/149‑1 (MQ=255)
                              ccTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGc                                                                                                                     >  6:297475/1‑149 (MQ=255)
                                         ccGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCg                                                                                                          <  7:290743/149‑1 (MQ=255)
                                          cGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGa                                                                                                         <  3:335128/149‑1 (MQ=255)
                                            aGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATc                                                                                                       >  5:119254/1‑148 (MQ=255)
                                                 gCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCaaaa                                                                                                   >  4:72181/1‑148 (MQ=255)
                                                 gCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCaaaa                                                                                                   >  3:284502/1‑148 (MQ=255)
                                                      agagGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGaa                                                                                             <  5:29253/149‑1 (MQ=255)
                                                       gagGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGGTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAAt                                                                                            <  2:302727/149‑1 (MQ=255)
                                                       gagGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTGGCAGGTGTGACGGCCTCTTCACGGCTACAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAAt                                                                                            <  7:201021/149‑1 (MQ=255)
                                                        agGTGCCAAAGTTGCGGGCGCTGGTGGCAGGGAGTGGCGAGCAGTTTCTGGTTTGCCGGCCTGTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGTCTTTAGCCGTGAAATGGCCGTCGAGGTCGATCAAAAAAGGAAt                                                                                            <  4:213613/148‑1 (MQ=255)
                                                          gTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGa                                                                                          >  3:318228/1‑148 (MQ=255)
                                                           tGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGaca                                                                                        <  3:169041/149‑1 (MQ=255)
                                                               aaaGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGacacac                                                                                     >  5:159707/1‑148 (MQ=255)
                                                                    tGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAACCCGCTCCTACGGGGACCGCCTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACtt                                                                               >  4:82303/1‑149 (MQ=255)
                                                                       gggCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACttttt                                                                            <  5:101832/149‑1 (MQ=255)
                                                                        ggCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAACCCCAACAACGGCGACCGGGAAAGCGCCGGATTATCCGGGAAAGGGCCCGAGAGGTACAAAAAAAAAAGAAAGGAACACGCACattttt                                                                           >  2:132088/1‑149 (MQ=255)
                                                                         gcgcTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACttttttt                                                                          <  4:284502/149‑1 (MQ=255)
                                                                         gcgcTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGAGCAAAAAAGGAATGACACACACACtgttttt                                                                          >  4:223686/1‑149 (MQ=255)
                                                                            cTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACcgc                                                                                                                                                  <  5:236592/75‑1 (MQ=255)
                                                                            cTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACcgc                                                                                                                                                  >  6:236592/1‑75 (MQ=255)
                                                                              ggtggAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCtt                                                                     >  3:157969/1‑149 (MQ=255)
                                                                              ggtggAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCtt                                                                     >  3:163715/1‑149 (MQ=255)
                                                                               tttgAATGGATTGGCGAGCAGTTGGAGGTTTGACGGCGTCTTCTCGGCTAAAGCCGCTGCTACGGTGAGCGCGTAGGCGCTGCTTTAGCCGGGAAAGGGCGGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTtg                                                                    <  1:132088/146‑1 (MQ=255)
                                                                               gtggAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTtg                                                                    <  4:63289/149‑1 (MQ=255)
                                                                               gtggAATGGATTGGCGAGCAGTTGAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTtgt                                                                   <  6:119254/149‑1 (MQ=255)
                                                                                tggAATGTATTGGCGTGCGTTTGCAGGTTTGACGGCTTCTTGGCGGATAAAGTTGCTCCTACGGGGACCGCGTTGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAATAGGAATTACACACTCACTTTTTTTCGCTtgt                                                                   <  7:234309/149‑1 (MQ=255)
                                                                                tggAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGTGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTtgt                                                                   <  8:353990/149‑1 (MQ=255)
                                                                                 ggAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTtgtg                                                                  <  3:20660/149‑1 (MQ=255)
                                                                                         ttGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGt                                                          <  7:197216/149‑1 (MQ=255)
                                                                                             ggcaGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAATGCCGCTCCTCGGGGGTCCGGGTAGGGGCGGCTTTAGCGGGGATAGTGCCGTCGAGGTCGATCAAAAAAGTAATGacaca                                                                                      <  2:272565/115‑1 (MQ=255)
                                                                                             gCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGaccca                                                                                      >  1:272565/1‑118 (MQ=255)
                                                                                                  cAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGTGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTccc                                                  <  5:297475/149‑1 (MQ=255)
                                                                                                  cAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTccc                                                  <  4:318228/149‑1 (MQ=255)
                                                                                                   aGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCa                                                 <  6:130092/149‑1 (MQ=255)
                                                                                                    gTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTTGGCGCGGCTTTAGCCGGGAAAGGGCCGTGGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAg                                                >  6:1962/1‑149 (MQ=255)
                                                                                                     ttGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACATTTTTTCGCTTGTGAGCAAAGTGGATTACCAGc                                               >  3:20900/1‑149 (MQ=255)
                                                                                                       tcAGGTTTTCCGGCCTTTTTTGTGCTTAAGCCTTTTCTACGTGGACCGCGTATGCGTGGTTTTATCCGGGAAATGGCCGTCGAGTTCGATTAAAAAATGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCt                                              <  8:235461/147‑1 (MQ=255)
                                                                                                          ggTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCa                                          <  3:72181/149‑1 (MQ=255)
                                                                                                           gTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTGGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCAt                                         >  4:31124/1‑149 (MQ=255)
                                                                                                                aCGGCCTCTTCGCGGCTAATGCCGCTCCTAGGGGGACCGCTTAGGCGCGGCTTTAGCGGGGAATGGGCGGTCGAGGTCTATCAAAAAAGTAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCTTGTaa                                     <  6:159707/148‑1 (MQ=255)
                                                                                                                aCGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGTCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTaa                                     <  5:209890/148‑1 (MQ=255)
                                                                                                                aCGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCGTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGt                                       >  1:280591/1‑147 (MQ=255)
                                                                                                                aCGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAg                                                                >  5:382785/1‑121 (MQ=255)
                                                                                                                aCGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAg                                                                <  6:382785/121‑1 (MQ=255)
                                                                                                                      ctttCCCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTcc                               >  6:131292/3‑148 (MQ=37)
                                                                                                                            cgGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTTAGCGTCCGAAccc                         >  1:138265/1‑148 (MQ=255)
                                                                                                                                tAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACATTTTTTAGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAAc                      >  6:363679/1‑147 (MQ=255)
                                                                                                                                 aaaGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAAc                      >  8:7367/1‑146 (MQ=255)
                                                                                                                                 aaaGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAAc                      <  7:7367/146‑1 (MQ=255)
                                                                                                                                 aaaGCCGCTCCGACGGGGACCGCGTAGGCGCGTCTATAGCCGGGAAAGGGCCGACGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCAAGTAAGCGTCCGAACCCAACCAg                   >  6:98717/1‑149 (MQ=255)
                                                                                                                                    gCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGaa                <  3:31124/149‑1 (MQ=255)
                                                                                                                                    gCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGaa                <  4:157969/149‑1 (MQ=255)
                                                                                                                                    gCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGaa                <  4:163715/149‑1 (MQ=255)
                                                                                                                                    gCCGCGCCTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAgg                  >  3:325514/1‑147 (MQ=38)
                                                                                                                                     ccGCTCCTACGTGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAg               <  1:369751/149‑1 (MQ=255)
                                                                                                                                     ccGCGCCTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGCGAAGAGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAg               >  6:254406/1‑149 (MQ=37)
                                                                                                                                       gcgccTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAg                                                >  8:87533/4‑114 (MQ=37)
                                                                                                                                       gcgccTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAg                                                <  7:87533/111‑1 (MQ=37)
                                                                                                                                           cTACGGGTACCGCGTAGGCGCGGCTTTAGCCGTGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGc          <  5:171339/148‑1 (MQ=255)
                                                                                                                                           cTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGc          >  5:64847/1‑148 (MQ=255)
                                                                                                                                            tACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGa        >  4:220250/1‑149 (MQ=255)
                                                                                                                                            tACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCAAACCAGGAAGCTAGCGa        >  4:225063/1‑149 (MQ=255)
                                                                                                                                             aCGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCaa                                                                                                     >  4:132186/1‑55 (MQ=255)
                                                                                                                                             aCGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCaa                                                                                                     <  3:132186/55‑1 (MQ=255)
                                                                                                                                               ggggACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCAGGTAAGCGTCCGAACCCCACCAGGCAGCTAGCGATGc     >  8:206787/1‑149 (MQ=255)
                                                                                                                                               ggggACAGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGc     <  3:223686/149‑1 (MQ=255)
                                                                                                                                                  gACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAGTGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCg  >  5:199592/1‑149 (MQ=255)
                                                                                                                                                  gACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCg  >  7:241075/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCG  >  NC_002947/6004851‑6005144

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