Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,238,465 G→T 6.5% intergenic (+168/‑115) ygjR → / → alx putative oxidoreductase YgjR/putative membrane‑bound redox modulator Alx

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,238,4650GT 6.5% 309.8 / 7.8 92intergenic (+168/‑115)ygjR/alxputative oxidoreductase YgjR/putative membrane‑bound redox modulator Alx
Reads supporting (aligned to +/- strand):  ref base G (39/47);  new base T (3/3);  total (42/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.04e-01

GGTATGTAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGG  >  NC_000913/3238322‑3238611
                                                                                                                                               |                                                                                                                                                  
ggTATGTAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTAt                                                                                                                                              >  2:15215/1‑150 (MQ=255)
    tGTAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAGGAATACGCAGACGTATTCCtt                                                                                                                                         >  2:388046/1‑150 (MQ=255)
     gTAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCtt                                                                                                                                         >  1:1155992/1‑150 (MQ=255)
      tAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCttt                                                                                                                                        >  1:982492/1‑150 (MQ=255)
       aaCAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGTCTGCGTATTCCtttc                                                                                                                                       <  2:355397/150‑2 (MQ=255)
       aaCAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCtttt                                                                                                                                       >  1:269174/1‑150 (MQ=255)
       aaCAGACCATTGACTGGCTAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCttttt                                                                                                                                      >  2:1246170/1‑150 (MQ=255)
        aCAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGTTGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGTCTGCGTATTCCtttcg                                                                                                                                      <  1:873968/150‑3 (MQ=255)
          agcccATTGACTGGCTGAAGGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTttg                                                                                                                                    <  1:349030/147‑1 (MQ=255)
          aGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTttg                                                                                                                                    <  1:245217/150‑1 (MQ=255)
          aGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTttg                                                                                                                                    >  1:1294575/1‑150 (MQ=255)
           gACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTttgt                                                                                                                                   <  2:188047/150‑1 (MQ=255)
           gACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAACACATCCGGTCACAGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTttgtt                                                                                                                                  >  2:1524104/1‑150 (MQ=255)
             ccATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGTCTGCGTATTCCtttcgggttg                                                                                                                                 <  2:434789/150‑8 (MQ=255)
                   aCTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTg                                                                                                                           <  2:1633763/150‑1 (MQ=255)
                   aCTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTg                                                                                                                           >  1:171398/1‑150 (MQ=255)
                       gCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTtgttgt                                                                                                                                <  2:1263269/141‑1 (MQ=255)
                       gCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTtgttgt                                                                                                                                >  1:1263269/1‑141 (MQ=255)
                        cTGAATGGTCTGTCATACTTTGTTCCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACCCATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCt                                                                                                                     <  1:3404/150‑1 (MQ=255)
                          gAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCtt                                                                                                                    >  1:149366/1‑150 (MQ=255)
                          gAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCtt                                                                                                                    >  2:1257206/1‑150 (MQ=255)
                           aaTGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTg                                                                                                                   <  2:1396753/150‑1 (MQ=255)
                             tGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGcc                                                                                                                 <  1:1444269/150‑1 (MQ=255)
                                  tGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAg                                                                                                            <  2:380809/150‑1 (MQ=255)
                                   gTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAgg                                                                                                           >  1:1032663/1‑150 (MQ=255)
                                     cATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCg                                                                                                         <  1:1446573/150‑1 (MQ=255)
                                       tACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAg                                                                                                       >  2:50998/1‑150 (MQ=255)
                                        aCTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAgg                                                                                                      <  1:15215/150‑1 (MQ=255)
                                        aCTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGCATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAgg                                                                                                      <  2:1593583/150‑1 (MQ=255)
                                           ttGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTg                                                                                                                           <  2:758112/126‑1 (MQ=255)
                                           ttGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTg                                                                                                                           >  1:758112/1‑126 (MQ=255)
                                               tACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGc                                                                                               <  2:1414414/150‑1 (MQ=255)
                                               tACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGc                                                                                               <  1:1285263/150‑1 (MQ=255)
                                                  cTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATa                                                                                            <  2:153780/150‑1 (MQ=255)
                                                  cTGCAAAGGGGAGTAACATCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATa                                                                                            >  2:631998/1‑150 (MQ=255)
                                                   tGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATaa                                                                                           >  1:858261/1‑150 (MQ=255)
                                                       aaGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAag                                                                                       >  2:249612/1‑150 (MQ=255)
                                                       aaGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAagc                                                                                      <  1:388046/150‑1 (MQ=255)
                                                           ggAGTAACTTCATTGCCGGTCGATTGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAagcagc                                                                                   <  1:1669518/150‑1 (MQ=255)
                                                           ggAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAagcagc                                                                                   <  1:1246170/150‑1 (MQ=255)
                                                           ggAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAagcagc                                                                                   <  2:581322/150‑1 (MQ=255)
                                                             aGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCgg                                                                                 <  1:1081506/150‑1 (MQ=255)
                                                              gTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGc                                                                                <  2:630035/150‑1 (MQ=255)
                                                                aaCTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTg                                                                              >  2:350917/1‑150 (MQ=255)
                                                                  cTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGAc                                                                            >  2:1478734/1‑150 (MQ=255)
                                                                        tGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTc                                                                      <  1:1282269/150‑1 (MQ=255)
                                                                         gCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTcc                                                                     <  2:1032663/150‑1 (MQ=255)
                                                                          ccGGTCGAACGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCg                                                                    >  2:1067985/1‑150 (MQ=255)
                                                                             gTCGATGGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACg                                                                 <  2:1348461/150‑1 (MQ=255)
                                                                                gATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTg                                                              >  2:957835/1‑150 (MQ=255)
                                                                                 aTCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTgg                                                             <  2:1269417/150‑1 (MQ=255)
                                                                                        ttACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGtttt                                                      <  1:508547/150‑1 (MQ=255)
                                                                                          aCGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGtttttt                                                    <  2:1532444/150‑1 (MQ=255)
                                                                                                tgtgAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTAtgtg                                              <  1:1025739/150‑1 (MQ=255)
                                                                                                tgtgAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTAtgtg                                              <  1:1019562/150‑1 (MQ=255)
                                                                                                 gtgAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATgtgt                                             <  2:1390977/150‑1 (MQ=255)
                                                                                                  tgAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGc                                                                                               <  2:1343388/99‑1 (MQ=255)
                                                                                                  tgAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGc                                                                                               >  1:1343388/1‑99 (MQ=255)
                                                                                                   gAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTaa                                           >  1:312059/1‑150 (MQ=255)
                                                                                                   gAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTaa                                           >  1:191805/1‑150 (MQ=255)
                                                                                                   gAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTaa                                           >  1:1597538/1‑150 (MQ=255)
                                                                                                      aaaaCACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGa                                        <  2:287328/150‑1 (MQ=255)
                                                                                                        aaCACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAAc                                      >  2:1490791/1‑150 (MQ=255)
                                                                                                        aaCACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAAc                                      >  2:1103136/1‑150 (MQ=255)
                                                                                                          cacaTCCGGTCACCGGGCAACCCGAAAGTAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAg                                          >  2:1482155/1‑144 (MQ=255)
                                                                                                          cacaTCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAg                                          <  1:1482155/144‑1 (MQ=255)
                                                                                                           acaTCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTc                                   <  2:149366/150‑1 (MQ=255)
                                                                                                            caTCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGtt                                                        <  2:649689/128‑1 (MQ=255)
                                                                                                            caTCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGtt                                                        >  1:649689/1‑128 (MQ=255)
                                                                                                            caTCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCt                                  <  2:171398/150‑1 (MQ=255)
                                                                                                             aTCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTa                                 <  1:576200/150‑1 (MQ=255)
                                                                                                             aTCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTa                                 <  2:1155992/150‑1 (MQ=255)
                                                                                                              tCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTAt                                >  2:1180053/1‑150 (MQ=255)
                                                                                                               ccGGTCACCGGGCAACCCGAAAGGACTACGCAGACGTCTTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATg                               <  2:643559/150‑1 (MQ=255)
                                                                                                               ccGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCTGTAATGCGAGGTCTATGCATAAAAAGCATCTGCTGACGTCTTCCTACGTTGTCCGTTTTTTTATTTGTAATGAAGTTCTATg                               <  2:1058201/150‑1 (MQ=255)
                                                                                                               ccGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATg                               >  2:1824/1‑150 (MQ=255)
                                                                                                               ccGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATg                               <  2:858261/150‑1 (MQ=255)
                                                                                                                 ggTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGaa                             >  1:503738/1‑150 (MQ=255)
                                                                                                                   tCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATa                           <  2:269174/150‑1 (MQ=255)
                                                                                                                    cACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATAc                          >  1:1462834/1‑150 (MQ=255)
                                                                                                                     aCCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACt                         >  1:325947/1‑150 (MQ=255)
                                                                                                                        gggCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTc                      <  2:1484715/150‑1 (MQ=255)
                                                                                                                             aCCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCAc                 >  2:564723/1‑150 (MQ=255)
                                                                                                                             aCCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCAc                 >  2:915402/1‑150 (MQ=255)
                                                                                                                               aaaaaaGGAATACGTCTGCGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGc               >  1:567072/4‑150 (MQ=255)
                                                                                                                                 aaaaGGAATACGTCTGCGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCg             >  1:1100294/2‑150 (MQ=255)
                                                                                                                                  aaaGGAATACGTCTGCGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGt            >  1:420708/1‑150 (MQ=255)
                                                                                                                                     ggAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGc         <  1:962117/150‑1 (MQ=255)
                                                                                                                                      gAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTc                      <  1:1157291/136‑1 (MQ=255)
                                                                                                                                      gAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTc                      >  2:1157291/1‑136 (MQ=255)
                                                                                                                                        aTACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTAt      >  2:1569213/1‑150 (MQ=255)
                                                                                                                                            gCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATgggg  <  1:957835/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GGTATGTAACAGACCATTGACTGGCTGAATGGTCTGTCATACTTTGTTACCTGCAAAGGGGAGTAACTTCATTGCCGGTCGATCGTCATTACGATGTGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGG  >  NC_000913/3238322‑3238611

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: