Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 390,088 G→T 7.2% S121I (AGC→ATC)  soxA → sarcosine oxidase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947390,0880GT 7.2% 243.3 / 4.4 77S121I (AGC→ATC) soxAsarcosine oxidase subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (35/36);  new base T (3/3);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

ATCGGCTCCAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTACAT  >  NC_002947/389946‑390226
                                                                                                                                              |                                                                                                                                           
aTCGGCTCCAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACGAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCATCatgct                                                                                                                                       >  1:115346/1‑147 (MQ=255)
      tCCAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCatgatcgggcg                                                                                                                                 >  1:308324/1‑143 (MQ=255)
      tCCAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATgccgcc                                                                                                                                 >  8:191693/1‑149 (MQ=255)
        cAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAATAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGccgccg                                                                                                                                <  8:296723/148‑1 (MQ=255)
        cAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCATCATGCTgccgccaa                                                                                                                               <  6:105597/149‑3 (MQ=255)
        cAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCATCATGCTgccgcc                                                                                                                                 <  2:115346/147‑1 (MQ=255)
         aGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCggg                                                                                                                              >  4:5237/1‑149 (MQ=255)
         aGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCggg                                                                                                                              >  3:116365/1‑149 (MQ=255)
          gCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGc                                                                                                                             <  4:116365/149‑1 (MQ=255)
              aGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCACCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCAGGATGCCGCCGGGCTTct                                                                                                                         <  1:281167/149‑1 (MQ=255)
               gCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAGCGACGTCAGGGGCATCGTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTcta                                                                                                                        >  4:298942/1‑149 (MQ=255)
                  aCCCAGATACCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCAACCAGCCCCAACGGTTGGCCGAACGTCAACACCGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCCCCGGGCGTgtacta                                                                                                                     >  2:161121/1‑149 (MQ=255)
                   cccAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTctacta                                                                                                                     >  4:326708/1‑148 (MQ=255)
                    ccAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCtactac                                                                                                                   >  7:57122/1‑148 (MQ=255)
                         tCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAAcc                                                                                                             >  7:26880/1‑149 (MQ=255)
                           cccAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCtt                                                                                                           >  5:185843/1‑149 (MQ=255)
                            ccAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTTCTACAAAACCtt                                                                                                           >  1:207810/1‑148 (MQ=255)
                            ccAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCtt                                                                                                           >  3:280538/1‑148 (MQ=255)
                            ccAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTc                                                                                                          <  4:280538/149‑1 (MQ=255)
                            ccAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTc                                                                                                          <  3:5237/149‑1 (MQ=255)
                             cAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACaaa                                                                                                                >  6:290630/1‑142 (MQ=255)
                             cAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACaaa                                                                                                                <  5:290630/142‑1 (MQ=255)
                              aaCGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCAt                                                                                                        <  8:39407/149‑1 (MQ=255)
                                   gcgcgcCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGTTTCTACTACAAAACCTTc                                                                                                          <  8:57408/142‑1 (MQ=255)
                                   gcgcgcCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGTTTCTACTACAAAACCTTc                                                                                                          >  7:57408/1‑142 (MQ=255)
                                        ccACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCaaa                                                                                              <  4:31783/149‑1 (MQ=255)
                                         cACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAAt                                                                                             >  8:179015/1‑149 (MQ=255)
                                          aCCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTTCTACAAAACCTTCATGTACCCCAAATc                                                                                            <  2:207810/149‑1 (MQ=255)
                                          aCCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATc                                                                                            >  5:130177/1‑149 (MQ=255)
                                             cAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGtt                                                                                         >  2:221806/1‑149 (MQ=255)
                                              aaCAGTCACTGTACGCGGGCCTTGTCGCCACCAGCACCATCTGTGGGCCGAACGTCAACAACGACGTGATGGGCATCCTCGGCAAGGTTGGCGTCATCATGATGCCGCCGTGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTc                                                                                        <  4:317554/149‑1 (MQ=255)
                                                 aGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTgg                                                                                     <  7:287113/149‑1 (MQ=255)
                                                  ggCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGa                                                                                    >  2:352170/1‑149 (MQ=255)
                                                   gCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGAt                                                                                   <  7:191693/149‑1 (MQ=255)
                                                    cACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGAt                                                                                   >  5:360906/1‑148 (MQ=255)
                                                           cgcgGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTAc                                                                            <  8:26880/148‑1 (MQ=255)
                                                           cgcgGGCCGTGTCGCCACCATCACCAAAGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGGTGGCGGCCGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTAc                                                                            <  6:352156/148‑1 (MQ=255)
                                                              gggCCTTGTCGCCACCAGCACCATCGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCAGGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACgaga                                                                        <  3:298942/149‑1 (MQ=255)
                                                              gggCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTc                                                                                                          >  4:283455/1‑115 (MQ=255)
                                                              gggCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTc                                                                                                          <  3:283455/115‑1 (MQ=255)
                                                               ggCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGaa                                                                       <  6:130177/149‑1 (MQ=255)
                                                                 ccTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGt                                                                     <  4:307380/149‑1 (MQ=255)
                                                                         ccaccaGCACCAACGGTTGGCCGAACGTCATCAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCg                                                              >  1:30840/1‑148 (MQ=255)
                                                                         ccaccaGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCg                                                              >  2:164119/1‑148 (MQ=255)
                                                                         ccaccaGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCg                                                              >  3:147923/1‑148 (MQ=255)
                                                                          caccaGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGt                                                             >  7:220197/1‑148 (MQ=255)
                                                                            ccaGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTaaa                                                          <  3:326190/149‑1 (MQ=255)
                                                                                  ccAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAgccgcc                                                    >  8:40533/1‑149 (MQ=255)
                                                                                   cAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGccgccg                                                   <  8:57122/149‑1 (MQ=255)
                                                                                   cAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGccgccg                                                   <  1:164119/149‑1 (MQ=255)
                                                                                           ggCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTggg                                            <  1:221806/148‑1 (MQ=255)
                                                                                           ggCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGAGGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTggg                                            <  1:161121/148‑1 (MQ=255)
                                                                                             ccGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAACGCCGCCGGCCTGGGCCg                                         >  2:52511/1‑149 (MQ=255)
                                                                                              cGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGt                                        <  4:164480/149‑1 (MQ=255)
                                                                                               gAACGTCAATAACGAAGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTg                                       <  6:157493/149‑1 (MQ=255)
                                                                                               gAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAATCCTTCATGTACCCCAAATAGTTCTGGAAGACTTACGAGAAGTACATCCGTAAAGCAGCCGGCCTGCGCCGTg                                       >  8:200278/1‑149 (MQ=255)
                                                                                                   gTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAACCCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGcgcc                                   >  8:265656/1‑149 (MQ=255)
                                                                                                   gTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGcgcc                                   >  6:144764/1‑149 (MQ=255)
                                                                                                    tCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGcgccg                                  <  7:179015/149‑1 (MQ=255)
                                                                                                      aacaacGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAg                                                                      <  3:101534/111‑1 (MQ=255)
                                                                                                      aacaacGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAg                                                                      >  4:101534/1‑111 (MQ=255)
                                                                                                         aacGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTc                                                                                        >  3:245195/1‑90 (MQ=255)
                                                                                                         aacGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTc                                                                                        <  4:245195/90‑1 (MQ=255)
                                                                                                         aacGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCa                             <  6:360906/149‑1 (MQ=255)
                                                                                                               gTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGa                       >  3:56373/1‑149 (MQ=255)
                                                                                                               gTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGa                       <  4:244234/149‑1 (MQ=255)
                                                                                                                  aTGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATc                     <  2:308324/148‑1 (MQ=255)
                                                                                                                      gCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGt                                        >  4:215683/1‑125 (MQ=255)
                                                                                                                      gCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGt                                        <  3:215683/125‑1 (MQ=255)
                                                                                                                      gCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTgcgc                                    >  8:258026/1‑129 (MQ=255)
                                                                                                                      gCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTgcgc                                    <  7:258026/129‑1 (MQ=255)
                                                                                                                              ggCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTAc          <  8:48739/147‑1 (MQ=255)
                                                                                                                              ggCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGa        <  7:302356/149‑1 (MQ=255)
                                                                                                                              ggCAAGGTTGCCGGCAGCATGATGCCGGCGGGCTTCTACTACAAAACGTTCATGTACCCCATATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGa        <  4:310456/149‑1 (MQ=255)
                                                                                                                                 aaGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTa     >  5:118015/1‑149 (MQ=255)
                                                                                                                                     ccggcggcATCATGCTGCCGACGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTACAt  >  6:246848/3‑148 (MQ=255)
                                                                                                                                      cggcggcATCATGCTGCCGCCGGGCTTCTACTATCAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTACAt  >  1:43109/2‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                           
ATCGGCTCCAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTACAT  >  NC_002947/389946‑390226

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: