Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 392,293 T→C 8.7% V856A (GTA→GCA)  soxA → sarcosine oxidase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947392,2930TC 8.7% 258.0 / 11.0 92V856A (GTA→GCA) soxAsarcosine oxidase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (41/43);  new base C (3/5);  total (44/48)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCCGCGTGTTCCGTATCTCGTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCATGAACCGCGAAGACTGCGT  >  NC_002947/392147‑392440
                                                                                                                                                    |                                                                                                                                                   
ggCCCGCGTGTTCCGTATCTCGTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTAc                                                                                                                                                   <  3:10498/149‑1 (MQ=255)
   ccGCGTGTTCCGTATCTCGTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACACCGTACTGCACCGAGACCATGCACGTACTGc                                                                                                                                                <  4:56057/149‑1 (MQ=255)
   acgcgTGTTCCGTATCTCGTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGTCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGc                                                                                                                                                <  3:221077/148‑1 (MQ=255)
    cgcgTGTTCCGTATCTCGTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCg                                                                                                                                               >  3:116985/1‑149 (MQ=255)
       gtgtTCCGTATCTCGTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGc                                                                                                                                            <  7:174470/149‑1 (MQ=255)
        tgtTCCGTATCTCGTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCa                                                                                                                                           <  4:116985/149‑1 (MQ=255)
          ttCCGTATCTCGTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCTTCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGCAGTACGTGCatg                                                                                                                                         <  1:275148/149‑3 (MQ=255)
          ttCCGTATCTCGTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCATGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGCAGTACGTGCatg                                                                                                                                         <  1:259592/149‑3 (MQ=255)
          ttCCGTATCTCGTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACt                                                                                                                                                  <  4:36156/140‑1 (MQ=255)
          ttCCGTATCTCGTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACt                                                                                                                                                  >  3:36156/1‑140 (MQ=255)
             cGTATCTCGTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCagag                                                                                                                                        >  6:256032/1‑147 (MQ=255)
             cGTATCTCGTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGCAGTACGTGCatggtc                                                                                                                                      <  3:227109/149‑6 (MQ=255)
               tATCTCGTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGAGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCAGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAgg                                                                                                                                    >  8:267661/1‑149 (MQ=255)
               tATCTCGTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGACAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAgg                                                                                                                                    <  8:198727/149‑1 (MQ=255)
                 tctcGTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGc                                                                                                                                  >  5:55664/1‑149 (MQ=255)
                     gTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTtca                                                                                                                              >  3:2078/1‑149 (MQ=255)
                     gTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTtc                                                                                                                               <  3:272646/148‑1 (MQ=255)
                        cACCGGCGAGCTGTCGTACGATGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTtcatc                                                                                                                            <  4:34116/148‑1 (MQ=255)
                         aCCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTcatcat                                                                                                                          >  2:106204/1‑149 (MQ=255)
                          ccGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCatcatc                                                                                                                         <  6:159983/149‑1 (MQ=255)
                           cGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCg                                                                                                                        <  5:335373/149‑1 (MQ=255)
                             gCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTg                                                                                                                      <  5:256032/149‑1 (MQ=255)
                              cGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTgg                                                                                                                     <  5:132179/149‑1 (MQ=255)
                                aGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGc                                                                                                                   <  5:122176/149‑1 (MQ=255)
                                          cGAAGTGAACGTGCAGGCCAACTACACCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACAcc                                                                                                          <  7:258793/148‑1 (MQ=255)
                                           gAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCatcatc                                                                                                                         >  1:91342/1‑132 (MQ=255)
                                           gAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCatcatc                                                                                                                         <  2:91342/132‑1 (MQ=255)
                                           gAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGa                                                                                                        >  3:57309/1‑149 (MQ=255)
                                           gAAGTGAACGTCCAGGCCAACTACGCCATGGGAGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACACCGTACGGCACAGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGacacaca                                                                                                        >  2:11940/1‑146 (MQ=255)
                                              gTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATgg                                                                                                     <  7:267661/149‑1 (MQ=255)
                                                 aaCGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTc                                                                                                  >  4:90414/1‑149 (MQ=255)
                                                   cGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCTAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGTAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGa                                                                                                <  5:285862/149‑1 (MQ=255)
                                                    gtccacgCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAAACAGACCCCGTACGGCACAGAGACCATGCACGTAAAGAGTGCAGAGAAGGGCTTCATCATAGTGGTCCAGGACAACGAAGGTACGa                                                                                                >  8:117231/7‑148 (MQ=255)
                                                       cAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACTCTGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGAt                                                                                               <  1:106204/148‑1 (MQ=255)
                                                          gCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAAccc                                                                                         >  6:204200/1‑149 (MQ=255)
                                                             aaCTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGga                                                                                      >  4:156801/1‑149 (MQ=255)
                                                             aaCTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGga                                                                                      >  3:208496/1‑149 (MQ=255)
                                                              aCTACGACATGGTCGTGCTGGAACAGATCATCGAGGCCTTCTAGAAGTACAACGTGAACCCGTACGTCAACGAGACCATGCACGTACTGCTTGAAGAGAAGGGCTTCATCATCGTGGTCCAGGACAAGGATGGTTCGATGAACCCGgac                                                                                     <  7:117231/149‑1 (MQ=255)
                                                                    ccATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCt                                                                                >  2:113650/1‑148 (MQ=255)
                                                                     cATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGc                                                                                                                                  <  7:278155/97‑1 (MQ=255)
                                                                     cATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGc                                                                                                                                  >  8:278155/1‑97 (MQ=255)
                                                                     cATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGTCACCGATGGTTCGATGAACCCGGACGACCTGa                                                                              <  1:212250/149‑1 (MQ=255)
                                                                     cATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGgac                                                                                     <  7:3912/142‑1 (MQ=255)
                                                                     cATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGgac                                                                                     >  8:3912/1‑142 (MQ=255)
                                                                      aTGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGaa                                                                             >  2:8341/1‑149 (MQ=255)
                                                                       tGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAAc                                                                            <  4:2078/149‑1 (MQ=255)
                                                                       tGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAAc                                                                            <  2:190864/149‑1 (MQ=255)
                                                                          gCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACACCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACAt                                                                          >  7:307356/1‑148 (MQ=255)
                                                                               cTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCGCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGAt                                                                                                       >  1:65075/1‑114 (MQ=255)
                                                                               cTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGAt                                                                                                       <  2:65075/114‑1 (MQ=255)
                                                                                tGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACAGCCTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACAACGATGGTTCGATGAACCCGGACGACCTGAACATGATCTgg                                                                   >  6:59757/1‑149 (MQ=255)
                                                                                   aaCAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGc                                                                <  4:57309/149‑1 (MQ=255)
                                                                                   aaCAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGc                                                                <  3:90414/149‑1 (MQ=255)
                                                                                        atcatcGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTgg                                                           >  2:190083/1‑149 (MQ=255)
                                                                                        atcatcGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTgg                                                           <  3:156801/149‑1 (MQ=255)
                                                                                        atcatcGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTgg                                                           >  7:194492/1‑149 (MQ=255)
                                                                                           atcGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCg                                                        >  1:248859/1‑149 (MQ=255)
                                                                                            tcGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACAt                                                                          >  6:169905/1‑130 (MQ=255)
                                                                                            tcGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACAt                                                                          <  5:169905/130‑1 (MQ=255)
                                                                                            tcGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGAGCAGAGAAGGGCTTCATCATAGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCg                                                        >  1:268739/1‑148 (MQ=255)
                                                                                                ggCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTg                                                                                                                      >  2:15904/1‑82 (MQ=255)
                                                                                                ggCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTg                                                                                                                      <  1:15904/82‑1 (MQ=255)
                                                                                                    ggCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGtt                                                             <  4:295726/135‑1 (MQ=255)
                                                                                                    ggCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGtt                                                             >  3:295726/1‑135 (MQ=255)
                                                                                                     gCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCg                                              >  5:122659/1‑149 (MQ=255)
                                                                                                         gaagTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTgg                                                           >  5:33980/1‑132 (MQ=255)
                                                                                                         gaagTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTgg                                                           <  6:33980/132‑1 (MQ=255)
                                                                                                                  ccTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGacac                                                                                                           <  3:178814/75‑1 (MQ=255)
                                                                                                                  ccTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGacac                                                                                                           >  4:178814/1‑75 (MQ=255)
                                                                                                                     gACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCTCAACAAGACGTTCTCGTGGATCGGcc                              >  2:89362/1‑149 (MQ=255)
                                                                                                                     gACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGcc                              >  8:109127/1‑149 (MQ=255)
                                                                                                                          cGTACGGGACCGAGACCATGCACGTACTGCGTGCAGAGAAGTGCTTCATCATAGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAAAATGAGCTGTTGCGTTGGCCGCAACAAGCCGTTCTCTAGCATCGGACTCCGt                         >  6:345223/1‑149 (MQ=255)
                                                                                                                          cGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGt                         >  1:201142/1‑149 (MQ=255)
                                                                                                                           gTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTg                        >  5:166128/1‑149 (MQ=255)
                                                                                                                           gTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTg                        >  5:166124/1‑149 (MQ=255)
                                                                                                                              cGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCa                     <  2:248859/149‑1 (MQ=255)
                                                                                                                              cGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCa                     <  6:166128/149‑1 (MQ=255)
                                                                                                                              cGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCa                     <  6:166124/149‑1 (MQ=255)
                                                                                                                               ggCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACTAGCCGTTCTCGTGGATCGGCCTGCGTGGCAt                    <  1:34504/149‑1 (MQ=255)
                                                                                                                                     cGAGACCATGCACGTACTGGGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGGTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCATGAAcc               <  1:11940/149‑1 (MQ=255)
                                                                                                                                     cGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCATGAAcc               <  4:208496/149‑1 (MQ=255)
                                                                                                                                       gagaCCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCATGAACcg              <  7:109127/149‑1 (MQ=255)
                                                                                                                                          aCCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCATGAACcgcg            >  7:215230/1‑148 (MQ=255)
                                                                                                                                           ccATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCATGAACCGCGaa          <  6:122659/149‑1 (MQ=255)
                                                                                                                                            tctGCACGCAGTACGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCATGAACcgcg            >  2:201861/3‑146 (MQ=255)
                                                                                                                                             ctGCACGCAGTACGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCTTGGATCg                                 <  8:274040/123‑1 (MQ=255)
                                                                                                                                             ctGCACGCAGTACGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCg                                 >  7:274040/2‑124 (MQ=255)
                                                                                                                                             ctGCACGCAGTACGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCg                                 <  6:164925/123‑1 (MQ=255)
                                                                                                                                             ctGCACGCAGTACGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCg                                 >  5:164925/2‑124 (MQ=255)
                                                                                                                                                cACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCATGAACCGCGAAGACtg     >  8:198990/1‑149 (MQ=255)
                                                                                                                                                  cGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCATGAACCGCGAAGACtgcg   <  3:142482/149‑1 (MQ=255)
                                                                                                                                                   gTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCATGAACCGCGAAGACtgcgt  >  3:303954/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGCCCGCGTGTTCCGTATCTCGTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCATGAACCGCGAAGACTGCGT  >  NC_002947/392147‑392440

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: