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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947390,0960AC 9.6% 222.8 / 16.0 83M123L (ATG→CTG) soxAsarcosine oxidase subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (41/34);  new base C (5/3);  total (46/37)
Fisher's exact test for biased strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.35e-01

TCGGCTCCAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTACATGAACCGGCA  >  NC_002947/389950‑390238
                                                                                                                                                  |                                                                                                                                              
tCGGCTCCAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATgc                                                                                                                                             >  1:264161/1‑150 (MQ=255)
     tCCAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGccgccg                                                                                                                                        >  2:261513/1‑150 (MQ=255)
     tCCAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGccgccg                                                                                                                                        >  2:910755/1‑150 (MQ=255)
       cAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCggg                                                                                                                                      >  1:203269/1‑150 (MQ=255)
        aGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGc                                                                                                                                     >  1:1326031/1‑150 (MQ=255)
           gAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTc                                                                                                                                  <  1:1187345/150‑1 (MQ=255)
           gAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTc                                                                                                                                  <  2:1326031/150‑1 (MQ=255)
              gCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTctac                                                                                                                               <  2:264161/150‑1 (MQ=255)
              gCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAAAGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTctac                                                                                                                               <  2:800715/150‑1 (MQ=255)
                tACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTctacta                                                                                                                             >  1:624901/1‑150 (MQ=255)
                 aCCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCtactac                                                                                                                            <  2:624901/150‑1 (MQ=255)
                          cccAGCTTGCGCGCCCCCGAGCAGGCACTGTACGCGGGCCTTGTCGCCACCAGTACCAACGGTTGGCCGAACGTCAACAACGACGTCATCGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCTGGCTTCTACTACAAAACCTTc                                                                                                                   <  2:1079040/150‑1 (MQ=255)
                          cccAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTc                                                                                                                   <  2:1270150/150‑1 (MQ=255)
                          cccAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTc                                                                                                                   <  2:613205/150‑1 (MQ=255)
                          cccAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTc                                                                                                                   >  1:1244497/1‑150 (MQ=255)
                           ccAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCa                                                                                                                  >  1:1252877/1‑150 (MQ=255)
                            cAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCGGCCGGGCTTCTACTACAAACCCTTCAt                                                                                                                 <  2:575768/150‑1 (MQ=255)
                            cAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGc                                                                                                                                     >  2:838156/1‑130 (MQ=255)
                            cAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGc                                                                                                                                     <  1:838156/130‑1 (MQ=255)
                                gTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTAc                                                                                                             <  2:1244497/150‑1 (MQ=255)
                                    gcgcCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGTGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCa                                                                                                         >  1:797661/1‑150 (MQ=255)
                                      gcCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCaaa                                                                                                       >  2:58715/1‑150 (MQ=255)
                                         aCCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCg                                                                                                    >  1:1346838/1‑150 (MQ=255)
                                         aCCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCg                                                                                                    >  1:546314/1‑150 (MQ=255)
                                           ccAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACGATAGGTTGGGCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGtt                                                                                                  >  2:760810/1‑150 (MQ=255)
                                            cAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACGCCAAATCGGTc                                                                                                 >  2:396812/1‑150 (MQ=255)
                                             aaCAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGtt                                                                                                  <  2:364375/148‑1 (MQ=255)
                                             aaCAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGtt                                                                                                  >  1:364375/1‑148 (MQ=255)
                                                aGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGa                                                                                             >  2:1389393/1‑150 (MQ=255)
                                                  gCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGACGCCAGGCTTCCACTGCAAAACCTTCATGTACCCCAAATCGTTCTGGATg                                                                                           >  1:1026600/1‑150 (MQ=255)
                                                     cTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTc                                                                                                                   <  2:869196/123‑1 (MQ=255)
                                                     cTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTc                                                                                                                   >  1:869196/1‑123 (MQ=255)
                                                        tACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTAc                                                                                     <  2:1252877/150‑1 (MQ=255)
                                                          cgcgGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACga                                                                                   <  2:1346838/150‑1 (MQ=255)
                                                             gggCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGaa                                                                                <  2:797661/150‑1 (MQ=255)
                                                             gggCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGaa                                                                                <  2:203269/150‑1 (MQ=255)
                                                                ccTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGa                                                                                          >  1:832844/1‑137 (MQ=255)
                                                                ccTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGa                                                                                          <  2:832844/137‑1 (MQ=255)
                                                                ccTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTa                                                                             >  1:472819/1‑150 (MQ=255)
                                                                  ttGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACa                                                                           >  1:1138641/1‑150 (MQ=255)
                                                                        ccaccaGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTAccc                                                                                                           >  1:409593/1‑112 (MQ=255)
                                                                        ccaccaGCACCAACGGTTGGCCGAACGTCAAAAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTAccc                                                                                                           <  2:409593/112‑1 (MQ=255)
                                                                         caccaGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGc                                                                                                                                     >  2:192509/1‑85 (MQ=255)
                                                                         caccaGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGc                                                                                                                                     <  1:192509/85‑1 (MQ=255)
                                                                         caccaGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACga                                                                                   <  1:490362/135‑1 (MQ=255)
                                                                         caccaGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACga                                                                                   >  2:490362/1‑135 (MQ=255)
                                                                            caGCAACAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAgc                                                                 >  1:836841/1‑150 (MQ=255)
                                                                              gCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGt                                                                              <  1:209656/135‑1 (MQ=255)
                                                                              gCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGt                                                                              >  2:209656/1‑135 (MQ=255)
                                                                                  cAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCAt                                                                                                                 <  1:224676/96‑1 (MQ=255)
                                                                                  cAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCAt                                                                                                                 >  2:224676/1‑96 (MQ=255)
                                                                                  cAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCAACCTCGGCAAGGTAGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACTTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCgg                                                           >  2:520557/1‑150 (MQ=255)
                                                                                       gTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTgg                                                      >  1:875199/1‑150 (MQ=255)
                                                                                       gTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTgg                                                      >  2:755556/1‑150 (MQ=255)
                                                                                              gAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAAACCGCAACTACCCCAAATCGTTCTGGATGACTTAAGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTgc                                               >  1:386297/1‑150 (MQ=255)
                                                                                                   tCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTgcgcgcc                                          <  1:760810/150‑4 (MQ=255)
                                                                                                              gTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGAt                               >  2:324696/1‑150 (MQ=255)
                                                                                                               tCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATc                              >  2:721356/1‑150 (MQ=255)
                                                                                                               tCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATc                              <  1:58715/150‑1 (MQ=255)
                                                                                                                 aTGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCt                            >  1:1031577/1‑150 (MQ=255)
                                                                                                                   gggCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGa                          >  2:913472/1‑150 (MQ=255)
                                                                                                                    ggCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGCTCCTGAc                         <  2:1026600/150‑1 (MQ=255)
                                                                                                                     gCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACa                        <  1:721356/150‑1 (MQ=255)
                                                                                                                        tCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCt                     <  2:1031577/150‑1 (MQ=255)
                                                                                                                         ccTCGGCAAGGTTGGCGGCAGCATGATGCCGACGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTa                    >  2:874801/1‑150 (MQ=255)
                                                                                                                          cTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTAc                   >  1:142738/1‑150 (MQ=255)
                                                                                                                             ggCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGAc                <  1:910755/150‑1 (MQ=255)
                                                                                                                                 aGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTACa            <  1:755556/150‑1 (MQ=255)
                                                                                                                                  gcccggcggcATCATGCTGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTACAt           >  1:416659/5‑150 (MQ=255)
                                                                                                                                    ccggcggcATCATGCTGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTgg                                                      <  2:1381482/103‑1 (MQ=255)
                                                                                                                                    ccggcggcATCATGCTGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTgg                                                      >  1:1381482/3‑105 (MQ=255)
                                                                                                                                     gcccggcggCATCATGCTGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGaa                                   >  1:888403/4‑123 (MQ=255)
                                                                                                                                     gcccggcggCATCATGCTGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGaa                                   <  2:888403/120‑1 (MQ=255)
                                                                                                                                     cggcggcATCATGCTGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAgccgcc                                                             <  1:934791/96‑1 (MQ=255)
                                                                                                                                     cggcggcATCATGCTGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAgccgcc                                                             >  2:934791/2‑97 (MQ=255)
                                                                                                                                     cggcggcATCATGCTGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTACATGaa        >  2:965623/2‑150 (MQ=255)
                                                                                                                                     cggcggcATCATGCTGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTACATGaa        >  1:1269074/2‑150 (MQ=255)
                                                                                                                                       gcggcAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTACATGAAcc      <  1:1389393/150‑1 (MQ=255)
                                                                                                                                         ggcAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTaaa                                                                   >  1:971812/1‑87 (MQ=255)
                                                                                                                                         ggcAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTaaa                                                                   <  2:971812/87‑1 (MQ=255)
                                                                                                                                          gcagcaTGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGccgccg                                                            <  2:129883/93‑1 (MQ=255)
                                                                                                                                          gcagcaTGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGccgccg                                                            >  1:129883/1‑93 (MQ=255)
                                                                                                                                           catcaTGCTGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTACATGAACCGGCa  <  1:965623/147‑1 (MQ=255)
                                                                                                                                               atgatgCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTggg                                                     <  1:534200/95‑1 (MQ=255)
                                                                                                                                               atgatgCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTggg                                                     >  2:534200/1‑95 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
TCGGCTCCAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTACATGAACCGGCA  >  NC_002947/389950‑390238

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: