Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 390,093 A→C 8.1% M123L (ATG→CTG)  soxA → sarcosine oxidase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947390,0930AC 8.1% 216.3 / 9.1 74M123L (ATG→CTG) soxAsarcosine oxidase subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (35/32);  new base C (4/2);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

GATCGGCTCCAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTACATGAACCGGCACTGCGA  >  NC_002947/389945‑390241
                                                                                                                                                    |                                                                                                                                                    
gATCGGCTCCAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCatga                                                                                                                                                      >  1:3022/1‑149 (MQ=255)
 aTCGGCTCCAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCatgat                                                                                                                                                     <  7:221243/149‑1 (MQ=255)
  tCGGCTCCAGCGAAGTTACCCAGATCCCCAACGTGCGCCCCACCCAACAGGCACTGTACGAGGGCGTTGTCGCCACCCGCCCCAACGGTGGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGTCAAGGTTGGCGGCAGCatgatg                                                                                                                                                    >  6:216756/1‑149 (MQ=255)
  tCGGCTCCAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGAGGGCCTTGTCGCCACCAGCACCAACGCTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGCCGGCAGCatcatg                                                                                                                                                    >  5:111680/1‑149 (MQ=255)
     gCTCCAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATgccg                                                                                                                                                 >  5:208260/1‑149 (MQ=255)
      cTCCAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATgccgc                                                                                                                                                <  1:148752/149‑1 (MQ=255)
      cTCAAGCGAAGCTACGCAGATCACCAACGTTCTCGCCACCCAACAGGCACTGTACGCGGGCGTTGTCGACACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATgccgc                                                                                                                                                <  8:239950/149‑1 (MQ=255)
       tCCAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCATCATGCTgccgc                                                                                                                                                >  8:8504/1‑148 (MQ=255)
         cAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCATCATGCTgccgccaa                                                                                                                                             <  7:8504/149‑3 (MQ=255)
          aGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCATCATGCTgccgcc                                                                                                                                               >  1:69310/1‑146 (MQ=255)
          aGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCATCATGCTgccgcc                                                                                                                                               <  2:69310/146‑1 (MQ=255)
           gCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGc                                                                                                                                           <  6:200156/149‑1 (MQ=255)
             gAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCtt                                                                                                                                         >  2:318054/1‑149 (MQ=255)
                   aCCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTctacta                                                                                                                                   >  4:273740/1‑149 (MQ=255)
                   aCCCAGATCCCCAACGTGCGCGCCACCAAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTctacta                                                                                                                                   >  4:344684/1‑149 (MQ=255)
                    cccAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTAGTCGCCACCAGCACCAACAGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTctacta                                                                                                                                   >  6:179225/1‑148 (MQ=255)
                        gATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACaaaa                                                                                                                              <  4:106611/149‑1 (MQ=255)
                          tCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCTCCACCAGCACCATCGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAAcc                                                                                                                            <  6:331762/149‑1 (MQ=255)
                           ccccAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCt                                                                                                                           <  1:118021/149‑1 (MQ=255)
                             ccAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCtt                                                                                                                          >  1:323956/1‑148 (MQ=255)
                             ccAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTc                                                                                                                         <  8:69625/149‑1 (MQ=255)
                             ccAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTc                                                                                                                         <  5:229494/149‑1 (MQ=255)
                                         ccACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCaaa                                                                                                             <  2:2677/149‑1 (MQ=255)
                                            cccAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCg                                                                                                          >  2:263472/1‑149 (MQ=255)
                                               aaaaGGCACTGTACGCGGGCCTTGTAGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGTCGGCAGCATGATTCAGCCGGGCTTCTACTACAAAACATTCATTTACCCCAAATCttta                                                                                                       >  5:146760/4‑148 (MQ=255)
                                               aaCAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTc                                                                                                       <  5:9143/149‑1 (MQ=255)
                                               aaCAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTc                                                                                                       <  5:147183/149‑1 (MQ=255)
                                                aCAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCt                                                                                                      >  6:152151/1‑149 (MQ=255)
                                                aCAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCt                                                                                                      <  3:89005/149‑1 (MQ=255)
                                                  aGGCACTGTACGCGGGCCTTGTCGCCAGCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTgg                                                                                                    <  6:111680/149‑1 (MQ=255)
                                                  aGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTgg                                                                                                    <  7:172371/149‑1 (MQ=255)
                                                   ggCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGa                                                                                                   >  5:237509/1‑149 (MQ=255)
                                                    gCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTc                                                                                                                         <  6:288274/126‑1 (MQ=255)
                                                    gCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTc                                                                                                                         >  5:288274/1‑126 (MQ=255)
                                                      aCTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGa                                                                                                >  1:291193/1‑149 (MQ=255)
                                                       cTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGAc                                                                                               >  5:268241/1‑149 (MQ=255)
                                                           aCGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTc                                                                                                                         >  2:93085/1‑119 (MQ=255)
                                                           aCGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTc                                                                                                                         <  1:93085/119‑1 (MQ=255)
                                                            cgcgGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTAc                                                                                           <  7:34167/148‑1 (MQ=255)
                                                             gcgGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACga                                                                                         <  1:318054/149‑1 (MQ=255)
                                                             gcgGGCCTTGTCGCCACCAGCACCAACGGTTCGCCGAACGACAACAACGACGTCATGGGAATCATCGGCAAGGTTCGCGGCAGCATCATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTAc                                                                                           >  8:18143/1‑147 (MQ=255)
                                                              ggggCCTGGTCTCGACCAGCAACAACTGTTGGACGAACGTCAACAACGACGTCATGGGATTCCTCGGCAAGGTTGGCGGAAGCATGTTGCCGACGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACgag                                                                                        <  5:179225/148‑1 (MQ=255)
                                                               gggCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACATAACCTGCATGTACCCCAAATGGTTGTGGATGAGTTAGgag                                                                                        >  2:334985/1‑148 (MQ=255)
                                                                ggCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGaa                                                                                      <  8:20551/149‑1 (MQ=255)
                                                                    ttGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTa                                                                                   >  4:195379/1‑148 (MQ=255)
                                                                      gTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACa                                                                                 <  3:273740/148‑1 (MQ=255)
                                                                       tCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCGAAATCGTTCTGGATGAGTTACGAGAAGTACAt                                                                                <  2:323956/148‑1 (MQ=255)
                                                                       tCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATc                                                                               >  4:246123/1‑149 (MQ=255)
                                                                        cGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATc                                                                               <  5:152151/148‑1 (MQ=255)
                                                                           caccaGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTc                                                                                                       <  1:248403/121‑1 (MQ=255)
                                                                           caccaGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTc                                                                                                       >  2:248403/1‑121 (MQ=255)
                                                                             ccaGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCGTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTaaa                                                                         <  2:3022/149‑1 (MQ=255)
                                                                                    cAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGccgccg                                                                  >  1:243011/1‑149 (MQ=255)
                                                                                             gCCTAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGc                                                          <  3:195379/148‑1 (MQ=255)
                                                                                              ccGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACAt                                                                                >  7:154725/1‑125 (MQ=255)
                                                                                              ccGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACAt                                                                                <  8:154725/125‑1 (MQ=255)
                                                                                               cGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTTCATCCGTAAAGCCGCCGGCCTGGGCCGt                                                       >  8:130767/1‑149 (MQ=255)
                                                                                                 caCGTCAACAACGACGTCCTCGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGAACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTgc                                                     <  8:328402/148‑1 (MQ=255)
                                                                                                    gTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGcgcc                                                  >  2:316137/1‑149 (MQ=255)
                                                                                                      caacaaCGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGcgccgc                                                >  6:28835/1‑149 (MQ=255)
                                                                                                          aacGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCGTCATGTACCCCAAATCGGTGTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTg                                              >  4:119402/1‑147 (MQ=255)
                                                                                                               cgTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAAc                                        >  6:207122/1‑148 (MQ=255)
                                                                                                                 tCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTCCCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGACGCTGCAGAACGaa                                     >  6:261950/1‑148 (MQ=255)
                                                                                                                    tGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCt                                  <  2:243011/149‑1 (MQ=255)
                                                                                                                          tCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGCCGTGCGCCGCTGCAGAACGATCCTGACAGc                            >  1:90004/1‑148 (MQ=255)
                                                                                                                           ccTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGc                            >  5:310358/1‑148 (MQ=255)
                                                                                                                            ctgggCATGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTgc                                                     <  3:191670/119‑1 (MQ=255)
                                                                                                                             tCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTAc                         >  3:285895/1‑149 (MQ=255)
                                                                                                                              cGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTAc                         >  3:193657/1‑148 (MQ=255)
                                                                                                                               ggCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGa                       >  4:279812/1‑149 (MQ=255)
                                                                                                                                  aaGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATTACTTACGAGAAGTAAATCCGTAAAGCCGCCGGGCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTa                    <  5:216756/149‑1 (MQ=255)
                                                                                                                                   aGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGTATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTAc                   <  3:344684/149‑1 (MQ=255)
                                                                                                                                     ggtGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTACAt                 <  6:208260/147‑1 (MQ=255)
                                                                                                                                      ttGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTACAt                 >  5:359880/1‑148 (MQ=255)
                                                                                                                                      ccggcggcATCATGCTGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTACAt                 >  7:71314/3‑148 (MQ=255)
                                                                                                                                       cggcggcATCATGCTGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTACATg                >  4:308161/2‑148 (MQ=255)
                                                                                                                                          cggcAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTACATGAAcc            >  8:350081/1‑149 (MQ=255)
                                                                                                                                             cagcaTGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTACATGAACCgg          >  3:314960/1‑148 (MQ=255)
                                                                                                                                                 atgatgCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTACATGAACCGGCACTg     >  2:15020/1‑149 (MQ=255)
                                                                                                                                                  cgatgCCGCCGGGCTTCTACTACAAAACCTACATGTACCCCAAAGCGTTCTGGATGACTTACGAGAAGTACATCCGTAACGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTACATGAACCGGCACTGc    <  6:310358/148‑1 (MQ=255)
                                                                                                                                                    atgCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTACATGAACCGGCACTGCGa  <  2:90004/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GATCGGCTCCAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTACATGAACCGGCACTGCGA  >  NC_002947/389945‑390241

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: