Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,681,133 C→G 5.0% V270L (GTG→CTG)  yhjJ ← peptidase M16 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,681,1330CG 5.0% 340.4 / 6.4 99V270L (GTG→CTG) yhjJpeptidase M16 family protein
Reads supporting (aligned to +/- strand):  ref base C (52/42);  new base G (3/2);  total (55/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

TTGCTGTTACTGGCACTTAACGCTTGCTGAACATGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATCCACGTTTCCCACCACC  >  NC_000913/3680985‑3681278
                                                                                                                                                    |                                                                                                                                                 
ttGCTGTTACTGGCACTTAACGCTTGCTGAACATGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCAcc                                                                                                                                                  >  2:416375/1‑150 (MQ=255)
   cTGTTACTGGCACTTAACGCTTGCTGAACATGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACcgcg                                                                                                                                               <  1:815540/150‑1 (MQ=255)
    tGTTACTGGCACTTAACGCTTGCTGAACATGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCgt                                                                                                                                              >  1:971121/1‑150 (MQ=255)
     gTTACTGGCACTTAACGCTTGCTGAACATGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCgtc                                                                                                                                             <  1:1389813/150‑1 (MQ=255)
       tACTGGCACTTAACGCTTGCTGAACATGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCGGTGCgtcag                                                                                                                                           <  2:1287652/150‑1 (MQ=255)
        aCTGGCACTTAACGCTTGCTGAACATGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCGGTGCgtcagg                                                                                                                                          <  2:613452/150‑2 (MQ=255)
        aCTGGCACTTAACGCTTGCTGAACATGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCgtcagg                                                                                                                                          >  2:406542/1‑149 (MQ=255)
             cACTTAACGCTTGCTGAACATGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCataa                                                                                                                                     >  1:672863/1‑150 (MQ=255)
             cACTTAACGCTTGCTGAACATGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCataa                                                                                                                                     >  1:48854/1‑150 (MQ=255)
             cACTTAACGCTTGCTGAACATGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCataa                                                                                                                                     >  1:42797/1‑150 (MQ=255)
              aCTTAACGCTTGCTGAACATGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCataat                                                                                                                                    <  2:780404/150‑1 (MQ=255)
                  aaCGCTTGCTGAACATGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTc                                                                                                                                >  1:1574716/1‑150 (MQ=255)
                     gCTTGCTGAACATGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCAcc                                                                                                                             <  1:1073164/150‑1 (MQ=255)
                        tGCTGAACATGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCt                                                                                                                          <  1:316723/150‑1 (MQ=255)
                              aCATGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCcgca                                                                                                                    <  1:854777/150‑1 (MQ=255)
                                 tGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCcgcacgc                                                                                                                 >  2:869241/1‑150 (MQ=255)
                                 tGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCcgcacgc                                                                                                                 >  1:354184/1‑150 (MQ=255)
                                    cAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGc                                                                                                              >  2:603369/1‑150 (MQ=255)
                                      gAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCgg                                                                                                            >  2:220334/1‑150 (MQ=255)
                                        aCAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGgag                                                                                                          <  2:971121/150‑1 (MQ=255)
                                         cAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGgaga                                                                                                         >  2:426858/1‑150 (MQ=255)
                                         cAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGgaga                                                                                                         >  1:1041156/1‑150 (MQ=255)
                                         cAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGgaga                                                                                                         >  2:593110/1‑150 (MQ=255)
                                          aGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGaa                                                                                                        >  2:1387229/1‑150 (MQ=255)
                                            cgcCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAg                                                                                                      >  1:248301/1‑150 (MQ=255)
                                               cTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTACGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGt                                                                                                   <  2:1366481/150‑1 (MQ=255)
                                                tCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCgtcagtca                                                                                                                                        <  1:543163/112‑1 (MQ=255)
                                                tCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCgtcagtca                                                                                                                                        >  2:543163/1‑112 (MQ=255)
                                                tCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTc                                                                                                  >  1:991817/1‑150 (MQ=255)
                                                 cgcgGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCg                                                                                                 <  2:354184/150‑1 (MQ=255)
                                                     ggCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCAc                                                                                             >  1:682098/1‑150 (MQ=255)
                                                      gCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCAcc                                                                                            >  1:1262993/1‑150 (MQ=255)
                                                       ccAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCg                                                                                           >  1:486577/1‑150 (MQ=255)
                                                        cAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCgg                                                                                          >  2:990463/1‑150 (MQ=255)
                                                        cAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCgg                                                                                          <  2:49414/150‑1 (MQ=255)
                                                             ccGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTg                                                                                     >  2:97087/1‑150 (MQ=255)
                                                              cGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTgg                                                                                    <  1:217203/150‑1 (MQ=255)
                                                              cGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTgg                                                                                    >  2:933377/1‑150 (MQ=255)
                                                              cGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTgg                                                                                    <  2:554693/150‑1 (MQ=255)
                                                               gCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGc                                                                                   <  2:244219/150‑1 (MQ=255)
                                                               gCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGc                                                                                   <  1:489621/150‑1 (MQ=255)
                                                                      aGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCAc                                                                            >  2:353792/1‑150 (MQ=255)
                                                                       gTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGGCGCACCGTGTTAGTTATAATAATGTCGTCTTCCGATAGCATCTGTTAAAGCGTCGGTATCGGATCTGGAGTTTCACg                                                                           <  2:1414981/150‑1 (MQ=255)
                                                                        tAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGt                                                                          <  1:1213080/150‑1 (MQ=255)
                                                                          gcgcAGCAGTGCGGCTGATTCACGAATTGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGttt                                                                        >  1:616394/1‑150 (MQ=255)
                                                                           cgcAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGGCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTg                                                                       <  2:1136101/150‑1 (MQ=255)
                                                                                  gTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCa                                                                <  2:1041156/150‑1 (MQ=255)
                                                                                    gCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGtt                                                                                <  2:797029/132‑1 (MQ=255)
                                                                                    gCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGtt                                                                                >  1:797029/1‑132 (MQ=255)
                                                                                      ggCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTc                                                            >  1:183170/1‑150 (MQ=255)
                                                                                      ggCTGAGTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTc                                                            <  1:616496/150‑1 (MQ=255)
                                                                                       gCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCg                                                           <  1:1495012/150‑1 (MQ=255)
                                                                                       gCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCg                                                           >  1:322200/1‑150 (MQ=255)
                                                                                         tGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGcc                                                         <  1:1436057/150‑1 (MQ=255)
                                                                                         tGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGcc                                                         <  1:636577/150‑1 (MQ=255)
                                                                                          gATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCa                                                        >  2:558870/1‑150 (MQ=255)
                                                                                          gATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCa                                                        <  2:248301/150‑1 (MQ=255)
                                                                                             tCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAAc                                                     <  1:974151/150‑1 (MQ=255)
                                                                                              cACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACg                                                    >  2:226900/1‑150 (MQ=255)
                                                                                              cACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACg                                                    >  1:1559926/1‑150 (MQ=255)
                                                                                              cACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACg                                                    >  1:1170458/1‑150 (MQ=255)
                                                                                               aCGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGt                                                   >  1:602257/1‑150 (MQ=255)
                                                                                               aCGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGt                                                   >  2:1268224/1‑150 (MQ=255)
                                                                                                 gAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGttt                                                 >  1:397958/1‑150 (MQ=255)
                                                                                                 gAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGttt                                                 <  2:84604/150‑1 (MQ=255)
                                                                                                  aaTCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGtttt                                                <  2:421275/150‑1 (MQ=255)
                                                                                                   aTCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTg                                               >  2:905164/1‑150 (MQ=255)
                                                                                                      ggCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGttt                                            >  2:129852/1‑150 (MQ=255)
                                                                                                          gCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTAttt                                        >  2:1136024/1‑150 (MQ=255)
                                                                                                          gCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTAttt                                        >  1:188986/1‑150 (MQ=255)
                                                                                                          gCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTAttt                                        >  1:604276/1‑150 (MQ=255)
                                                                                                          gCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTAttt                                        >  1:185332/1‑150 (MQ=255)
                                                                                                           ccACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTg                                       >  1:1094231/1‑150 (MQ=255)
                                                                                                             aCGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGt                                     <  2:1212315/150‑1 (MQ=255)
                                                                                                              cGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTc                                    >  1:915562/1‑150 (MQ=255)
                                                                                                                    aTCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAAc                              >  1:129083/1‑150 (MQ=255)
                                                                                                                      cccACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCg                            >  2:103286/1‑150 (MQ=255)
                                                                                                                       ccACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGa                           >  1:73427/1‑150 (MQ=255)
                                                                                                                         acaTGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAgc                         <  1:353792/150‑1 (MQ=255)
                                                                                                                            tGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGTGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAgcgcg                      >  2:28503/1‑150 (MQ=255)
                                                                                                                            tGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAgcgcg                      <  2:616394/150‑1 (MQ=255)
                                                                                                                                agataACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATc                  <  2:1627114/150‑1 (MQ=255)
                                                                                                                                 gataACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATcc                 >  1:755547/1‑150 (MQ=255)
                                                                                                                                 gataACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATcc                 >  1:1487770/1‑150 (MQ=255)
                                                                                                                                 cataACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATcc                 <  1:1101905/149‑1 (MQ=255)
                                                                                                                                  ataACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATCCa                <  2:115780/150‑1 (MQ=255)
                                                                                                                                      ttatgacTGACGCGGTGCGTCAGTCATAAGACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTGCCAGTTCGCCAAACGTTTTGTGTATTTGGTCGACACCCGGGCGCGCAGCCACGtt            >  1:1136823/7‑150 (MQ=255)
                                                                                                                                      ccGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCCCGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATCCACGtt            >  1:1599873/1‑150 (MQ=255)
                                                                                                                                        ggTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATCCACGTTTc          <  2:1447158/150‑1 (MQ=255)
                                                                                                                                        ggTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATCCACGTTTc          <  2:48854/150‑1 (MQ=255)
                                                                                                                                        ggTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATCCACGTTTc          <  2:42797/150‑1 (MQ=255)
                                                                                                                                        ggTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATCCACGTTTc          >  2:1179064/1‑150 (MQ=255)
                                                                                                                                         gTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATCCACGTTTcc         <  2:486577/150‑1 (MQ=255)
                                                                                                                                         gTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATCCACGTTTcc         <  2:129083/150‑1 (MQ=255)
                                                                                                                                          gacTGACGCGGTGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATCCACGTTTccc        >  2:518188/3‑150 (MQ=255)
                                                                                                                                          gacTGACGCGGTGCGTCAGTCATAAAACTAACCGCTTCCCCACGCAGCGGAGAAACCGACGGCACCCGAGCTGGCGTTTCACGTTTGCCATTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATCCACGTTTccc        >  2:330095/3‑150 (MQ=255)
                                                                                                                                           ccTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATCCACGTTTCcca       <  2:1263799/150‑1 (MQ=255)
                                                                                                                                               aCGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCACATCCACGTTTCccaccac   <  2:1574716/150‑1 (MQ=255)
                                                                                                                                                cGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATCCACGTTTCccaccacc  <  2:392919/150‑1 (MQ=255)
                                                                                                                                                cGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATCCACGTTTCccaccacc  >  2:198893/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TTGCTGTTACTGGCACTTAACGCTTGCTGAACATGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATCCACGTTTCCCACCACC  >  NC_000913/3680985‑3681278

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