Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,984,873 A→G W379R (TGG→CGG)  aslA ← putative Ser‑type periplasmic non‑aryl sulfatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,984,8730AG100.0% 283.7 / NA 78W379R (TGG→CGG) aslAputative Ser‑type periplasmic non‑aryl sulfatase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (40/38);  total (40/38)

TTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGAAGAGATCTGCCAGATCGACAATACCGTCAGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCAGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCGGAGGTAAAGACGATCAGCGTGTTATCAAGCTGACCGTTTTTCTCCAGTGTTTTATACAGATTAGC  >  NC_000913/3984739‑3985008
                                                                                                                                      |                                                                                                                                       
tttCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGAAGAGATCTGCCAGATCGACAATACCGTCAGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTc                                                                                                                                     >  2:940701/1‑139 (MQ=255)
tttCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGAAGAGATCTGCCAGATCGACAATACCGTCAGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTc                                                                                                                                     >  2:1387856/1‑139 (MQ=255)
     cTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGAAGAGATCTGCCAGATCGACAATACCGTCAGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAAcc                                                                                                                                >  1:502710/1‑139 (MQ=255)
               cccTGCCAGATCCAGCGCGGTAGGGAAGAGATCTGCCAGATCGACAATACCGTCAGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGcacca                                                                                                                      <  1:1387856/139‑1 (MQ=255)
               cccTGCCAGATCCAGCGCGGTAGGGAAGAGATCTGCCAGATCGACAATACCGTCAGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGcacca                                                                                                                      >  2:1570449/1‑139 (MQ=255)
                ccTGCCAGATCCAGCGCGGTAGGGAAGAGATCTGCCAGATCGACAATACCGTCAGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCg                                                                                                                           <  1:1313447/133‑1 (MQ=255)
                ccTGCCAGATCCAGCGCGGTAGGGAAGAGATCTGCCAGATCGACAATACCGTCAGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCg                                                                                                                           >  2:1313447/1‑133 (MQ=255)
                    ccAGATCCAGCGCGGTAGGGAAGAGATCTGCCAGATCGACAATACCGTCAGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGaa                                                                                                                 <  1:1440718/139‑1 (MQ=255)
                      aGATCCAGCGCGGTAGGGAAGAGATCTGCCAGATCGACAATACCGTCAGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACg                                                                                                               <  1:468846/139‑1 (MQ=255)
                         tCCAGCGCGGTAGGGAAGAGATCTGCCAGATCGACAATACCGTCAGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACgggg                                                                                                            >  2:856607/1‑139 (MQ=255)
                              cgcgGTAGGGAAGAGATCTGCCAGATCGACAATACCGTCAGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGt                                                                                                       >  2:678562/1‑139 (MQ=255)
                              cgcgGTAGGGAAGAGATCTGCCAGATCGACAATACCGTCAGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGt                                                                                                       >  1:390397/1‑139 (MQ=255)
                                   tAGGGAAGAGATCTGCCAGATCGACAATACCGTCAGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTg                                                                                                  >  1:322625/1‑139 (MQ=255)
                                      ggAAGAGATCTGCCAGATCGACAATACCGTCAGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCAGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTgcgg                                                                                               <  1:1178550/139‑1 (MQ=255)
                                             aTCTGCCAGATCGACAATACCGTCAGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACt                                                                                        <  2:1064633/139‑1 (MQ=255)
                                             aTCTGCCAGATCGACAATACCGTCAGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACt                                                                                        >  1:55541/1‑139 (MQ=255)
                                                tGCCAGATCGACAATACCGTCAGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGaa                                                                                                                 <  1:938886/111‑1 (MQ=255)
                                                tGCCAGATCGACAATACCGTCAGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGaa                                                                                                                 >  2:938886/1‑111 (MQ=255)
                                                tGCCAGATCGACAATACCGTCAGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCg                                                                                     >  2:688478/1‑139 (MQ=255)
                                                     gATCGACAATACCGTCAGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTc                                                                                >  2:1206218/1‑139 (MQ=255)
                                                      aTCGACAATACCGTCAGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTcc                                                                               >  2:1581321/1‑139 (MQ=255)
                                                      aTCGACAATACCGTCAGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTcc                                                                               >  2:186228/1‑139 (MQ=255)
                                                          aCAATACCGTCAGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTc                                                                           >  1:1360284/1‑139 (MQ=255)
                                                                 cGTCAGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATc                                                                    <  1:856607/139‑1 (MQ=255)
                                                                  gTCAGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCg                                                                   <  2:1360284/139‑1 (MQ=255)
                                                                  gTCAGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCg                                                                   <  1:1596488/139‑1 (MQ=255)
                                                                  gTCAGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCg                                                                   <  1:678562/139‑1 (MQ=255)
                                                                   tCAGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCtt                                                                                 <  2:1480363/124‑1 (MQ=255)
                                                                   tCAGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCtt                                                                                 >  1:1480363/1‑124 (MQ=255)
                                                                   tCAGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCgg                                                                  >  2:1199150/1‑139 (MQ=255)
                                                                   tCAGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCgg                                                                  >  1:1303528/1‑139 (MQ=255)
                                                                     aGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCGGAg                                                                >  1:1397439/1‑139 (MQ=255)
                                                                     aGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCGGAg                                                                >  2:91638/1‑139 (MQ=255)
                                                                      gATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCGGAgg                                                               <  2:838795/139‑1 (MQ=255)
                                                                      gATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCGGAgg                                                               <  1:1570449/139‑1 (MQ=255)
                                                                                 ttGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCGGAGGTAAAGACGATc                                                    <  1:91638/139‑1 (MQ=255)
                                                                                 ttGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCGGAGGTAAAGACGATc                                                    <  2:613563/139‑1 (MQ=255)
                                                                                 ttGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCGGAGGTAAAGACGATc                                                    >  1:1253223/1‑139 (MQ=255)
                                                                                     atcatcCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCGGAGGTAAAGACGATCAGCg                                                <  2:264900/139‑1 (MQ=255)
                                                                                         tcCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCGGAGGTAAAGACGATCAGCGTGtt                                            <  2:390397/139‑1 (MQ=255)
                                                                                         tcCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCGGAGGTAAAGACGATCAGCGTGtt                                            >  1:1445180/1‑139 (MQ=255)
                                                                                           ccTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCg                                                                   <  2:1570787/114‑1 (MQ=255)
                                                                                           ccTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCg                                                                   >  1:1570787/1‑114 (MQ=255)
                                                                                              ttCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTc                                                                                      <  1:369833/92‑1 (MQ=255)
                                                                                              ttCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTc                                                                                      >  2:369833/1‑92 (MQ=255)
                                                                                               tCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGcaccac                                                                                                                     >  1:54434/1‑60 (MQ=255)
                                                                                               tCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGcaccac                                                                                                                     <  2:54434/60‑1 (MQ=255)
                                                                                               tCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCgg                                                                             >  1:839421/1‑100 (MQ=255)
                                                                                               tCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCgg                                                                             <  2:839421/100‑1 (MQ=255)
                                                                                                   gTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGaa                                                                                                                 <  2:1051987/60‑1 (MQ=255)
                                                                                                   gTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGaa                                                                                                                 >  1:1051987/1‑60 (MQ=255)
                                                                                                    tAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTc                                                                                >  2:191605/1‑92 (MQ=255)
                                                                                                    tAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTc                                                                                <  1:191605/92‑1 (MQ=255)
                                                                                                    tAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCGGAGGTAAAGACGATCAGCGTGTTATCAAGCTGAc                                 >  1:1111935/1‑139 (MQ=255)
                                                                                                    tAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCGGAGGTAAAGACGATCAGCGTGTTATCAAGCTGAc                                 <  1:186228/139‑1 (MQ=255)
                                                                                                    tAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCGGAGGTAAAGACGATCAGCGTGTTATCAAGCTGAc                                 >  1:1458593/1‑139 (MQ=255)
                                                                                                    tAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCGGAGGTAAAGACGATCAGCGTGTTATCAAGCTGAc                                 >  2:1214296/1‑139 (MQ=255)
                                                                                                           aaGTCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCGGAGGTAAAGACGATCAGCGTGTTATCAAGCTGACCGttttt                          >  2:1274593/1‑139 (MQ=255)
                                                                                                              tCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCg                                                                                                                           <  1:723752/39‑1 (MQ=255)
                                                                                                              tCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCg                                                                                                                           >  2:723752/1‑39 (MQ=255)
                                                                                                              tCGGTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCGGAGGTAAAGACGATCAGCGTGTTATCAAGCTGACCGTTTTtctc                       <  2:1253223/139‑1 (MQ=255)
                                                                                                                ggTACGCGAACGCGGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCGGAGGTAAAGACGATCAGCGTGTTATCAAGCTGACCGTTTTTCTCCa                     <  1:1305805/139‑1 (MQ=255)
                                                                                                                ggTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCGGAGGTAAAGACGATCAGCGTGTTATCAAGCTGACCGTTTTTCTCCa                     <  1:1199150/139‑1 (MQ=255)
                                                                                                                ggTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCGGAGGTAAAGACGATCAGCGTGTTATCAAGCTGACCGTTTTTCTCCa                     <  2:1445180/139‑1 (MQ=255)
                                                                                                                ggTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCGGAGGTAAAGACGATCAGCGTGTTATCAAGCTGACCGTTTTTCTCCa                     >  1:700353/1‑139 (MQ=255)
                                                                                                                 gTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTAt                                                                     >  2:1364805/1‑90 (MQ=255)
                                                                                                                 gTACGCGAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTAt                                                                     <  1:1364805/90‑1 (MQ=255)
                                                                                                                    cgcgAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCGGAGGTAAAGACGATCAGCGTGTTATCAAGCTGACCGTTTTTCTCCAgtgt                 <  2:1111935/139‑1 (MQ=255)
                                                                                                                     gcgAACGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTTTCAGAGGAAAAGACGATAAGTGTGTTATCAAGATGAACGTTTTTCTCCAGTGtt                <  1:1206218/139‑1 (MQ=255)
                                                                                                                        aaCGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACg                                                                                                               <  1:1468373/41‑1 (MQ=255)
                                                                                                                        aaCGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACg                                                                                                               >  2:1468373/1‑41 (MQ=255)
                                                                                                                        aaCGCCGCCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCGGAGGTAAAGACGATCAGCGTGTTATCAAGCTGACCGTTTTTCTCCAGTGTTTta             >  1:938118/1‑139 (MQ=255)
                                                                                                                          cgccgcCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCGGAGGTAAAGACGATCAGCGTGTTATCAAGCTGACCGTTTTTCTCCAGTGTTTtata           <  1:1214296/139‑1 (MQ=255)
                                                                                                                          cgccgcCTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCGGAGGTAAAGACGATCAGCGTGTTATCAAGCTGACCGTTTTTCTCCAGTGTTTtata           <  2:1430950/139‑1 (MQ=255)
                                                                                                                            ccgccTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCGGAGGTAAAGACGATCAGCGTGTTATCAAGCTGACCGTTTTTCTCCAGTGTTTTATACa         >  1:561380/1‑139 (MQ=255)
                                                                                                                               ccTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCGGAGGTAAAGACGATCAGCGTGTTATCAAGCTGACCGTTTTTCTCCAGTGTTTTATACAGAt      >  2:493718/1‑139 (MQ=255)
                                                                                                                                cTTCCCGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCGGAGGTAAAGACGATCAGCGTGTTATCAAGCTGACCGTTTTTCTCCAGTGTTTTATACAGAtt     <  2:1458593/139‑1 (MQ=255)
                                                                                                                                   cccGGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCGGAGGTAAAGACGATCAGCGTGTTATCAAGCTGACCGTTTTTCTCCAGTGTTTTATACAGATTAGc  <  2:502710/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                       
TTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGAAGAGATCTGCCAGATCGACAATACCGTCAGATTTACGCGGTTGGATCATCCCTTTCCAGTAAACGAAAGTCGGTACGCGAACGCCGCCTTCCCAGGTCGAACCTTTCGCACCACGGAACGGGGTGCGTCCGTGCGGCGGTACTTCGGCTTCCGGTCCGTTATCGGAGGTAAAGACGATCAGCGTGTTATCAAGCTGACCGTTTTTCTCCAGTGTTTTATACAGATTAGC  >  NC_000913/3984739‑3985008

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: