Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,135,883 G→C 11.3% D229E (GAC→GAG)  PP_4521 ← aerotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,135,8830GC11.3% 244.8 / 20.3 89D229E (GAC→GAG) PP_4521aerotaxis receptor
Reads supporting (aligned to +/- strand):  ref base G (38/41);  new base C (5/5);  total (43/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.50e-01

AGTCTCGCTGCGCTGCCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCCAGCAGCAGACTTGCCA  >  NC_002947/5135736‑5136031
                                                                                                                                                   |                                                                                                                                                    
aGTCTCGCTGCGCTGCCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGt                                                                                                                                                     <  7:141546/149‑1 (MQ=255)
aGTCTCGCTGCGCTGCCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGt                                                                                                                                                     >  1:228571/1‑149 (MQ=255)
aGTCTCGCTGCGCTGCCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTGTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGt                                                                                                                                                     >  5:221159/1‑149 (MQ=255)
   ctcGCTGCGCTGCCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGGATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTc                                                                                                                                                    >  7:199008/1‑147 (MQ=255)
    tcGCTGCGCTGCCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAg                                                                                                                                                 >  6:358023/1‑149 (MQ=255)
    tcGCTGCGCTGCCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGCTCGac                                                                                                                                                 >  7:25363/1‑148 (MQ=255)
        tgcgctgcCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTgg                                                                                                                                             >  8:168227/1‑149 (MQ=255)
          cgctgcCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGGCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTgg                                                                                                                                             <  4:357393/147‑1 (MQ=255)
          cgctgcCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGCTCGACCTGGcc                                                                                                                                           <  3:349978/149‑1 (MQ=255)
           gctgcCGGTCGAGTATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGccc                                                                                                                                          <  5:299024/149‑1 (MQ=255)
           gctgcCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGCTCGACCTGGcca                                                                                                                                          <  4:318813/149‑2 (MQ=255)
           gctgcCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGCTCGACCTGGcca                                                                                                                                          <  8:25363/149‑2 (MQ=255)
            ctgcCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCt                                                                                                                                         >  2:158742/1‑149 (MQ=255)
            ctgcCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCt                                                                                                                                         <  6:283907/149‑1 (MQ=255)
            ctgcCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGCTCGACCTGGccat                                                                                                                                         <  1:366039/149‑3 (MQ=255)
                     gAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAg                                                                                                                                     >  6:47550/1‑144 (MQ=255)
                     gAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAg                                                                                                                                     <  5:47550/144‑1 (MQ=255)
                       gCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTggg                                                                                                                              <  1:100912/149‑1 (MQ=255)
                       gCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTggg                                                                                                                              >  4:13020/1‑149 (MQ=255)
                        cATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGt                                                                                                                             <  5:251622/149‑1 (MQ=255)
                        cATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGCTCGACCTGGcca                                                                                                                                          >  7:339789/1‑135 (MQ=255)
                        cATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGCTCGACCTGGcca                                                                                                                                          <  8:339789/136‑2 (MQ=255)
                           cTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGt                                                                                                                          >  4:77843/1‑149 (MQ=255)
                              cgcCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTc                                                                                                                         <  3:304565/147‑1 (MQ=255)
                              cgcATGGGCATCGGAGAAGTCGGTAGATTGAGCCGCCCGCGACTTTACCCCGACCGCGGCATCCACCAGGGGTGTCAGGGGCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCTAGCTGGCCCTTCTGGCCTGGGTTGTc                                                                                                                         <  8:199008/147‑1 (MQ=255)
                                ccTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTg                                                                                                                     >  8:87787/1‑149 (MQ=255)
                                ccTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTAGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTg                                                                                                                     >  4:29758/1‑149 (MQ=255)
                                 cTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTg                                                                                                                     >  8:69424/1‑148 (MQ=255)
                                   gggCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGCCCCCGGCAGCAACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTGCTTCGCGGATGATGCCCCGGGGGAGCTGGCCGTGCGGGCGTGGGTTggggatac                                                                                                                    >  4:104192/1‑140 (MQ=255)
                                    ggCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAg                                                                                                                  <  3:124891/148‑1 (MQ=255)
                                     gCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGt                                                                                                                <  7:168227/149‑1 (MQ=255)
                                       aTCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGt                                                                                                                >  3:241776/1‑147 (MQ=255)
                                       aTCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCa                                                                                                              <  3:77843/149‑1 (MQ=255)
                                           gcccaGTCGGCAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCCGGCGGGTCAGGGCCGTTGTCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGGTGTCGCTGTAGGTCAgcgc                                                                                                          <  2:259461/145‑1 (MQ=255)
                                              aaGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGccac                                                                                                       >  5:67868/1‑149 (MQ=255)
                                                gTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGt                                                                                                                                                     >  4:19077/1‑101 (MQ=255)
                                                gTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGt                                                                                                                                                     <  3:19077/101‑1 (MQ=255)
                                                gTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGc                                                                                                                       >  4:205658/1‑131 (MQ=255)
                                                gTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGc                                                                                                                       <  3:205658/131‑1 (MQ=255)
                                                gTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGccacc                                                                                                      >  1:278593/1‑148 (MQ=255)
                                                   ggAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACg                                                                                                  >  1:85733/1‑149 (MQ=255)
                                                   ggAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACg                                                                                                  >  8:271290/1‑149 (MQ=255)
                                                      aGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGAc                                                                                               <  3:13020/149‑1 (MQ=255)
                                                         ttGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTgg                                                                                            <  6:290058/149‑1 (MQ=255)
                                                            cgccgcCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTa                                                                                          >  4:130449/1‑148 (MQ=255)
                                                            cgccgcCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTa                                                                                          <  2:48838/148‑1 (MQ=255)
                                                             gccgccCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTaaa                                                                                        <  7:271290/149‑1 (MQ=255)
                                                             gccgccCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTa                                                                                          >  8:66406/1‑147 (MQ=255)
                                                              ccgcccgcGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAg                                                                                       <  7:192982/149‑1 (MQ=255)
                                                               cgcccgcGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAgg                                                                                      >  1:160629/1‑149 (MQ=255)
                                                                  ccgcGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTAAGCGCCACCAACGGCCTGGTAACGGCCt                                                                                   >  8:327736/1‑149 (MQ=255)
                                                                   cgcgACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCtt                                                                                  <  2:58665/149‑1 (MQ=255)
                                                                      gACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTggg                                                                               <  4:325498/149‑1 (MQ=255)
                                                                          ttACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGc                                                                           >  3:186118/1‑149 (MQ=255)
                                                                              cccgaccccgGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTa                                                                                                                   >  3:259558/1‑105 (MQ=255)
                                                                              cccgaccccgGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTa                                                                                                                   <  4:259558/105‑1 (MQ=255)
                                                                              cccgaccccgGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAg                                                                       >  7:205048/1‑149 (MQ=255)
                                                                              cccgaccccgGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAg                                                                       >  1:333222/1‑149 (MQ=255)
                                                                                cgaccccgGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGACGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGCCTGGTAAAGGCCTTGGGATGCTCAGc                                                                      >  8:342039/1‑148 (MQ=255)
                                                                                   ccccgGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAAc                                                                  <  2:316170/149‑1 (MQ=255)
                                                                                      cgGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTg                                                               <  3:29758/149‑1 (MQ=255)
                                                                                       ggCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGCTGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGc                                                              <  5:358023/149‑1 (MQ=255)
                                                                                        gCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGcg                                                             <  6:67868/149‑1 (MQ=255)
                                                                                          aTCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGcgac                                                           <  1:158742/149‑1 (MQ=255)
                                                                                             caccaGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGaa                                                        >  4:22469/1‑149 (MQ=255)
                                                                                              accaGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGc                                                              >  7:117593/1‑142 (MQ=255)
                                                                                              accaGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGc                                                              <  8:117593/142‑1 (MQ=255)
                                                                                              accaGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAAt                                                       >  3:185687/1‑149 (MQ=255)
                                                                                               ccaGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATg                                                      <  6:165160/149‑1 (MQ=255)
                                                                                                   gCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCt                                                  <  7:66406/149‑1 (MQ=255)
                                                                                                           gggCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGc                                          <  4:185687/149‑1 (MQ=255)
                                                                                                              ccGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGGCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCgg                                       <  7:327736/149‑1 (MQ=255)
                                                                                                                 tttGCCGACGGGCTTCTTCGCGTATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGc                                    <  8:73813/149‑1 (MQ=255)
                                                                                                                 tttGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGc                                    >  5:68634/1‑149 (MQ=255)
                                                                                                                    gCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTc                                  <  2:333222/148‑1 (MQ=255)
                                                                                                                       gACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGAcc                              <  7:87787/149‑1 (MQ=255)
                                                                                                                         cGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCa                            <  3:130449/149‑1 (MQ=255)
                                                                                                                            gCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTa                         >  6:120382/1‑149 (MQ=255)
                                                                                                                             cttcttCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTAc                        >  6:254391/1‑149 (MQ=255)
                                                                                                                                 ttcGATGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCAGCGTTGACGCTGTAGCTCAGCGCCACCAACGGACTGGTAAAGTCCTTGGGATGCTCAGCCAACGTGc                                                              >  7:315429/1‑107 (MQ=255)
                                                                                                                                     cTGATGACGGCCAGGTCGAGCTGGCCCTTCATGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCATCTGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCcag                <  7:342039/149‑1 (MQ=255)
                                                                                                                                        tgaagggCCAGCTCGACCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCcagcag             >  7:204792/6‑149 (MQ=255)
                                                                                                                                        tgaagggCCAGCTCGACCTGGCCCTTCAGGCCTGGGTTGACGCTGTAGGTCAGCGCCACCACCGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCcaccag             >  7:139386/6‑149 (MQ=255)
                                                                                                                                        atgatgGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCcagcag             <  4:241776/149‑1 (MQ=255)
                                                                                                                                        atgatgGCCAGGGCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCcagcag             <  2:278593/149‑1 (MQ=255)
                                                                                                                                         gaagggCCAGCTCGACCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCcagcag             >  1:20645/5‑148 (MQ=255)
                                                                                                                                              gCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCGCGACCCAGTACCAAGCCAGCAGCAGa         <  2:85733/147‑1 (MQ=255)
                                                                                                                                               ccAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCCAGCAGCAGACtt      >  5:221862/1‑149 (MQ=255)
                                                                                                                                                cAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTAGATGTAGGTCAGCGACACCAACGGACTGGTAAAGGCCTTGGGATGCTAAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAATCCAGCAGCAGACTTg     >  3:71213/1‑149 (MQ=255)
                                                                                                                                                  ggTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCCAGCAGCAGACTTGcc   <  5:254391/149‑1 (MQ=255)
                                                                                                                                                   gTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCCAGCAGCAGACTTGCCa  >  1:105130/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGTCTCGCTGCGCTGCCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCCAGCAGCAGACTTGCCA  >  NC_002947/5135736‑5136031

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