Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 388,725 T→A 7.3% L198Q (CTG→CAG)  soxB → sarcosine oxidase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947388,7250TA 7.3% 205.6 / 7.5 70L198Q (CTG→CAG) soxBsarcosine oxidase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (28/36);  new base A (2/3);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.34e-01

GAACACTGCCCAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCAAGCT  >  NC_002947/388576‑388872
                                                                                                                                                     |                                                                                                                                                   
gAACACTGCCCAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCt                                                                                                                                                     >  1:784560/1‑150 (MQ=255)
 aaCACTGCCCAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTg                                                                                                                                                    >  2:905526/1‑150 (MQ=255)
      tGCCCAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTGCGCCCGCGCTGCCGACGCCCTGGGCGt                                                                                                                                               <  2:855057/150‑1 (MQ=255)
        cccAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCg                                                                                                                                             >  1:704000/1‑150 (MQ=255)
           aGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGTCGACGCCCAGGGCGTCggca                                                                                                                                          <  2:687919/150‑4 (MQ=255)
           aGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGAGGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCAGGGCGTCggca                                                                                                                                          <  2:270653/150‑4 (MQ=255)
             gTTGCGGCCGAAATCTCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCAGGGCGTCggcagc                                                                                                                                        <  2:1353611/150‑6 (MQ=255)
                   gCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGAc                                                                                                                                           >  2:1576963/1‑141 (MQ=255)
                   gCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGAc                                                                                                                                           <  1:1576963/141‑1 (MQ=255)
                     cGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCa                                                                                                                                >  1:1621871/1‑150 (MQ=255)
                      gAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCaa                                                                                                                               >  2:836769/1‑150 (MQ=255)
                      gAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCaa                                                                                                                               >  1:324579/1‑150 (MQ=255)
                       aaaTCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCaaa                                                                                                                              >  2:1818200/1‑150 (MQ=255)
                        aaTCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCCGGGGCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAAc                                                                                                                             <  2:1720289/150‑1 (MQ=255)
                            ccGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGaa                                                                                                                         <  2:1611348/150‑1 (MQ=255)
                            ccGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGaa                                                                                                                         <  1:349005/150‑1 (MQ=255)
                                  cTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATc                                                                                                                   <  2:1091630/150‑1 (MQ=255)
                                  cTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATc                                                                                                                   <  2:331195/150‑1 (MQ=255)
                                   tGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCg                                                                                                                  >  1:649922/1‑150 (MQ=255)
                                   tGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCg                                                                                                                  >  2:1633885/1‑150 (MQ=255)
                                           tCGAAGAAAACCCGCTAACCGAACCTCGGCGCCACCGTACGGAGCAGTGGTAGCGTTGCCCGTCACGATGACGTGGCCTGGGGCTTCGCCCACGCTGCCGACGCCCCGGGGGGCGACCTGAACCAGAAAACCGAAGTTATCGGCTTACGc                                                                                                          >  2:624628/1‑150 (MQ=255)
                                             gaagaaCACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCaa                                                                                                        >  1:1743781/1‑150 (MQ=255)
                                                   cacCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGaaaa                                                                                                  >  1:693404/1‑150 (MQ=255)
                                                         cTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGgcgc                                                                                            <  2:404276/150‑1 (MQ=255)
                                                         cTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGgcgc                                                                                            >  1:1609479/1‑150 (MQ=255)
                                                                   cTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGgt                                                                                  <  1:1876055/150‑1 (MQ=255)
                                                                    tCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGaa                                                                                                                         <  1:1681833/110‑1 (MQ=255)
                                                                    tCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGaa                                                                                                                         >  2:1681833/1‑110 (MQ=255)
                                                                    tCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGTAGTGATCGGCTTCCGCAAGGAAAACTGCGCGGTCATCGgtg                                                                                 <  2:1621871/150‑1 (MQ=255)
                                                                    tCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGgtg                                                                                 >  1:915491/1‑150 (MQ=255)
                                                                     cGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGgtgt                                                                                <  2:704000/150‑1 (MQ=255)
                                                                     cGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGgtgt                                                                                <  1:958273/150‑1 (MQ=255)
                                                                     cGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGgtgt                                                                                <  1:905526/150‑1 (MQ=255)
                                                                     cGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGgtgt                                                                                <  2:649922/150‑1 (MQ=255)
                                                                     cGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGgtgt                                                                                <  1:1389002/150‑1 (MQ=255)
                                                                         gcCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGaa                                                                            >  1:1422036/1‑150 (MQ=255)
                                                                           cACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGCTCCCGTGGTCTCGGGGTCCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAAc                                                                          >  1:1008261/1‑150 (MQ=255)
                                                                            aCCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAg                                                                                                                        >  1:984401/1‑103 (MQ=255)
                                                                            aCCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAg                                                                                                                        <  2:984401/103‑1 (MQ=255)
                                                                                tACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACgg                                                                                               >  2:987696/1‑124 (MQ=255)
                                                                                tACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACgg                                                                                               <  1:987696/124‑1 (MQ=255)
                                                                                 aCAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTCGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACcaacaa                                                                    <  2:1324217/150‑1 (MQ=255)
                                                                                   aGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAg                                                                  <  1:1498098/150‑1 (MQ=255)
                                                                                   aGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAg                                                                  <  2:205187/150‑1 (MQ=255)
                                                                                    gcgcCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAgg                                                                 <  1:1725780/150‑1 (MQ=255)
                                                                                       ccGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCtt                                                              <  2:1008261/150‑1 (MQ=255)
                                                                                        cGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTc                                                             <  1:1004620/150‑1 (MQ=255)
                                                                                        cGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTc                                                             >  2:748359/1‑150 (MQ=255)
                                                                                              ggCGTTGCCCGTCACGATGCCGTGGCCTGCCGACACGCCGGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGgc                                                       <  2:1609479/150‑1 (MQ=255)
                                                                                                   tGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTACGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCaa                                                  <  2:784560/150‑1 (MQ=255)
                                                                                                       cGTCACGATGCCATGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATACAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAgcgc                                              >  1:1494601/1‑150 (MQ=255)
                                                                                                             gATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGc                                        >  2:294964/1‑150 (MQ=255)
                                                                                                              aTGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCg                                       <  1:1633885/150‑1 (MQ=255)
                                                                                                                gCCGTGGCCTGTGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTg                                     <  1:1062393/150‑1 (MQ=255)
                                                                                                                gCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTg                                     <  2:915491/150‑1 (MQ=255)
                                                                                                                   gTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGc                                                                >  1:1298675/1‑120 (MQ=255)
                                                                                                                   gTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGc                                                                <  2:1298675/120‑1 (MQ=255)
                                                                                                                      gCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGACAAGCGCGTCGGCGTGGTCAcc                               <  2:1345139/150‑1 (MQ=255)
                                                                                                                          ggggCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCAccgccg                           <  1:836769/150‑1 (MQ=255)
                                                                                                                           gggCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCTAGCAAACCGAAGTGATCGGCTTCCGCaa                                                                                                        >  1:1290319/1‑72 (MQ=255)
                                                                                                                           gggCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCTAGCAAACCGAAGTGATCGGCTTCCGCaa                                                                                                        <  2:1290319/72‑1 (MQ=255)
                                                                                                                           gggCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCgg                          <  2:1422036/150‑1 (MQ=255)
                                                                                                                           gggCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCgg                          <  1:294964/150‑1 (MQ=255)
                                                                                                                                tCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACt                     <  2:693404/150‑1 (MQ=255)
                                                                                                                                    ccGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCgg                 >  1:996825/1‑150 (MQ=255)
                                                                                                                                     cgcgCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCggg                >  1:1801106/1‑150 (MQ=255)
                                                                                                                                         aggtcGACGCCCAGGGCGTCGGCCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGcaca            >  2:1404052/5‑150 (MQ=255)
                                                                                                                                          ggtcGACGCCCAGGGCGTCGGCCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACAt           >  1:421007/4‑150 (MQ=255)
                                                                                                                                              gACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGcc       >  1:1314317/1‑150 (MQ=255)
                                                                                                                                                 gCCCTGGGCGTCGACCTGATCCAGCAAGCCCAGGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCAAg    >  1:1807914/1‑150 (MQ=255)
                                                                                                                                                   ccTGGGCGTCGACCTGATCCAGCAAACGGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGACACCGCCGGTAACTCCGGGCACATAGCCAAGCt  >  2:399489/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GAACACTGCCCAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCAAGCT  >  NC_002947/388576‑388872

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: