Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 392,295 C→G 8.5% L857V (CTG→GTG) ‡ soxA → sarcosine oxidase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947392,2950CG 8.5% 277.8 / 12.1 92L857V (CTG→GTG) ‡soxAsarcosine oxidase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (41/42);  new base G (3/5);  total (45/47)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.42e-01

GGCCCGCGTGTTCCGTATCTCGTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCATGAACCGCGAAGACTGCGTGC  >  NC_002947/392147‑392442
                                                                                                                                                      |                                                                                                                                                   
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                                                             >  7:194492/1‑149 (MQ=255)
                                                                                        atcatcGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTgg                                                             <  3:156801/149‑1 (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                                                             <  6:33980/132‑1 (MQ=255)
                                                                                                         gaagTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTgg                                                             >  5:33980/1‑132 (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)
                                                                                                                                                     aCTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCATGAACCGCGAAGACtgcgtgc  >  7:119510/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
GGCCCGCGTGTTCCGTATCTCGTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCATGAACCGCGAAGACTGCGTGC  >  NC_002947/392147‑392442

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: