Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 392,297 G→A 8.8% L857L (CTG→CTA) ‡ soxA → sarcosine oxidase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947392,2970GA 8.8% 255.4 / 10.5 91L857L (CTG→CTA) ‡soxAsarcosine oxidase subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (42/41);  new base A (3/5);  total (45/46)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.77e-01

CCGCGTGTTCCGTATCTCGTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCATGAACCGCGAAGACTGCGTGCGT  >  NC_002947/392150‑392444
                                                                                                                                                     |                                                                                                                                                   
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)
          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                                                                                          >  3:208496/1‑149 (MQ=255)
                                                          aaCTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGga                                                                                          >  4:156801/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                                                                    <  3:90414/149‑1 (MQ=255)
                                                                                aaCAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGc                                                                    <  4:57309/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                                                                              <  5:169905/130‑1 (MQ=255)
                                                                                         tcGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACAt                                                                              >  6:169905/1‑130 (MQ=255)
                                                                                         tcGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGAGCAGAGAAGGGCTTCATCATAGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCg                                                            >  1:268739/1‑148 (MQ=255)
                                                                                             ggCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTg                                                                                                                          <  1:15904/82‑1 (MQ=255)
                                                                                             ggCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTg                                                                                                                          >  2:15904/1‑82 (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                                                                                                               >  4:178814/1‑75 (MQ=255)
                                                                                                               ccTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGacac                                                                                                               <  3:178814/75‑1 (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)
                                                                                                                                                    tGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCATGAACCGCGAAGACTgcgtgcgt  >  2:229721/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCGCGTGTTCCGTATCTCGTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCATGAACCGCGAAGACTGCGTGCGT  >  NC_002947/392150‑392444

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: