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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947390,0910GT 9.0% 245.9 / 10.5 78S121I (AGC→ATC) soxAsarcosine oxidase subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (39/32);  new base T (5/2);  total (44/34)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

CCTGCAGATCGGCTCCAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTACATGAACC  >  NC_002947/389942‑390234
                                                                                                                                                     |                                                                                                                                               
ccTGCAGATCGGCTCCAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGgcag                                                                                                                                                 <  1:1151546/150‑1 (MQ=255)
 cTGCAGATCGGCTCCAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGgcagc                                                                                                                                                >  1:575768/1‑150 (MQ=255)
        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                                                                                                                              <  2:1326031/150‑1 (MQ=255)
                   gAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTc                                                                                                                              <  1:1187345/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:613205/150‑1 (MQ=255)
                                  cccAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTc                                                                                                               <  2:1270150/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                                                                                                                                 <  1:838156/130‑1 (MQ=255)
                                    cAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGc                                                                                                                                 >  2:838156/1‑130 (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                                                                               >  2:490362/1‑135 (MQ=255)
                                                                                 caccaGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACga                                                                               <  1:490362/135‑1 (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                                                  >  2:755556/1‑150 (MQ=255)
                                                                                               gTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTgg                                                  >  1:875199/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                          <  1:58715/150‑1 (MQ=255)
                                                                                                                       tCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATc                          >  2:721356/1‑150 (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                                                               <  2:971812/87‑1 (MQ=255)
                                                                                                                                                 ggcAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTaaa                                                               >  1:971812/1‑87 (MQ=255)
                                                                                                                                                  gcagcaTGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGccgccg                                                        >  1:129883/1‑93 (MQ=255)
                                                                                                                                                  gcagcaTGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGccgccg                                                        <  2:129883/93‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
CCTGCAGATCGGCTCCAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCCAACAGGCACTGTACGCGGGCCTTGTCGCCACCAGCACCAACGGTTGGCCGAACGTCAACAACGACGTCATGGGCATCCTCGGCAAGGTTGGCGGCAGCATGATGCCGCCGGGCTTCTACTACAAAACCTTCATGTACCCCAAATCGTTCTGGATGACTTACGAGAAGTACATCCGTAAAGCCGCCGGCCTGGGCCGTGCGCCGCTGCAGAACGATCCTGACAGCTACGACTACATGAACC  >  NC_002947/389942‑390234

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: