Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,033,624 A→G 5.2% intergenic (+11/+33) idi → / ← lysS isopentenyl diphosphate isomerase/lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,033,6240AG 5.2% 466.7 / 11.0 156intergenic (+11/+33)idi/lysSisopentenyl diphosphate isomerase/lysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base A (69/77);  new base G (6/2);  total (76/80)
Fisher's exact test for biased strand distribution p-value = 1.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-02

Only 100 of 157 total aligned reads displayed.

GGATTATCAATGGTGTGATTTAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGC  >  NC_000913/3033475‑3033772
                                                                                                                                                     |                                                                                                                                                    
ggATTATCAATGGTGTGATTTAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGa                                                                                                                                                      >  1:123717/1‑150 (MQ=255)
    tacccaTGGTGTGATTTAGCAGATGTATAACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACAtt                                                                                                                                                  <  2:1605353/145‑1 (MQ=255)
      tCAATGGTGTGATTGAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCATAAAACGATTATCTGCATTAACCCAGCTTAAATAAAAAAACCCCGACATTTg                                                                                                                                                <  1:112082/150‑1 (MQ=255)
        aaTTGTGTGATTTAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGcc                                                                                                                                              <  2:1955590/150‑1 (MQ=255)
        aaTGGTGTGATTTAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGcc                                                                                                                                              <  2:210490/150‑1 (MQ=255)
         aTGGTGTGATTTAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCTGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCg                                                                                                                                             <  1:1173290/150‑1 (MQ=255)
            gtgtGATTTAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCgggg                                                                                                                                          >  1:1116762/1‑150 (MQ=255)
            gtgtGATTTAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCgggg                                                                                                                                          >  1:887919/1‑150 (MQ=255)
            gggtgATTTAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCgggg                                                                                                                                          <  2:551318/148‑1 (MQ=255)
               tgATTTAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTtg                                                                                                                                       <  1:1945708/150‑1 (MQ=255)
                 aTTTAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTtgtg                                                                                                                                     <  1:413138/150‑1 (MQ=255)
                 aTTTAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTtgtg                                                                                                                                     <  2:634117/150‑1 (MQ=255)
                   ttAGCATATGTATTACACGGTATTGATGCCACGCAGTGGGCGTTCAGTCCGTGGATGTTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAg                                                                                                                                   <  1:54277/150‑1 (MQ=255)
                   ttAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGGCAAATGTCGGGGttttttt                                                                                                                                   <  2:137912/150‑6 (MQ=255)
                     aGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCa                                                                                                                                 <  1:1556499/150‑1 (MQ=255)
                        aGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATaa                                                                                                                              <  1:874435/150‑1 (MQ=255)
                           tGTATTACAAGTTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCTGCTTAAATAAAAAAACCCCGACATTTGCCGTGGTTGTGAGCATAACGt                                                                                                                           <  1:967600/150‑1 (MQ=255)
                             tATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTaa                                                                                                                         >  2:1419659/1‑150 (MQ=255)
                              aTTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAAt                                                                                                                        <  2:885014/150‑1 (MQ=255)
                               ttACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATg                                                                                                                       <  2:144936/150‑1 (MQ=255)
                                tccacGGTCTTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGc                                                                                                                      <  1:984723/148‑1 (MQ=255)
                                tACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGc                                                                                                                      <  1:1736301/150‑1 (MQ=255)
                                   acGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCAAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTa                                                                                                                   >  1:1921173/1‑150 (MQ=255)
                                    cGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGCCAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTAGCCGGGGTTGTGAGCATAACGTAATGCTTAt                                                                                                                  >  1:1578314/1‑150 (MQ=255)
                                    cGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTAt                                                                                                                  >  1:142290/1‑150 (MQ=255)
                                    cGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTAt                                                                                                                  >  2:663641/1‑150 (MQ=255)
                                     ggTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTAtt                                                                                                                 <  1:1857598/150‑1 (MQ=255)
                                      gTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTAttt                                                                                                                >  2:1530782/1‑150 (MQ=255)
                                      gTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCATAttt                                                                                                                >  2:1202318/1‑150 (MQ=255)
                                        tttGATGCCACTCCGTGGGCGTTCATTCCGTGGATGTTTATGAAGGCGACAAATCGCGAAGCCAGAAAACGATTATCCGAATTTACCCAGCTTAAATAAAAAAACCCCGAAATTTGGCGGGGTTGTGAGCATTACGTAATGCTTATTTTa                                                                                                              <  2:34242/149‑1 (MQ=255)
                                        aTTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTa                                                                                                              >  2:1538257/1‑150 (MQ=255)
                                          tGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTAcc                                                                                                            <  1:640534/150‑1 (MQ=255)
                                          tGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGGGATGCAGCCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTAcc                                                                                                            <  1:387934/150‑1 (MQ=255)
                                           gATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCg                                                                                                           >  2:260739/1‑150 (MQ=255)
                                            aTGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGCCATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCgg                                                                                                          >  1:1995403/1‑150 (MQ=255)
                                            aTGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCgg                                                                                                          >  2:1562222/1‑150 (MQ=255)
                                             tGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGa                                                                                                         >  2:436765/1‑150 (MQ=255)
                                                cACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGc                                                                                                      >  1:109032/1‑150 (MQ=255)
                                                  cGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCAt                                                                                                    <  1:740814/150‑1 (MQ=255)
                                                    ccatGGGCGTTCAGTCCGTGGATGGGGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTAACTGCCTTTACCAAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCAGATCGTAATGCGTATTTTACCGGCCGCATcc                                                                                                  >  2:346036/4‑149 (MQ=255)
                                                    ccGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACACAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCg                                                                                                  >  2:1807611/1‑150 (MQ=255)
                                                         ggCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATCACGTAATGCTTATTTTACCGGACGCATCGCCggg                                                                                             <  1:254744/150‑1 (MQ=255)
                                                         ggCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCggg                                                                                             <  2:123717/150‑1 (MQ=255)
                                                              tCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAg                                                                                        <  2:1937740/150‑1 (MQ=255)
                                                               cAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACAATTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGTAACAGa                                                                                       >  1:1326633/1‑150 (MQ=255)
                                                                aGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGaa                                                                                      <  2:707697/150‑1 (MQ=255)
                                                                   ccGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCAGTGAACAGAATaa                                                                                   >  2:355368/1‑150 (MQ=255)
                                                                         aTGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTcgc                                                                             <  2:777634/150‑1 (MQ=255)
                                                                          tGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTcgcg                                                                            <  2:812302/150‑1 (MQ=255)
                                                                          tGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTcgcg                                                                            <  2:1927967/150‑1 (MQ=255)
                                                                              gATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATg                                                                        >  1:1021331/1‑150 (MQ=255)
                                                                                tGCAGTCGACAAATCGCGAAGCCAGAAAACGATTAGCTGCATTTACCGAGCTTAAATAAAAATACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGTCGCATCGCCGGGAACAGAATAACGTCGCGGATGGt                                                                      <  2:1130131/150‑1 (MQ=255)
                                                                                tGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGt                                                                      <  2:1880161/150‑1 (MQ=255)
                                                                                 gCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTTATGCTTATTATACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTa                                                                     >  2:826276/1‑150 (MQ=255)
                                                                                    ggCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCAGGAACAGAATAACGTCGCGGATGGTATgt                                                                  >  2:530376/1‑149 (MQ=255)
                                                                                    ggCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAAAGCTTATTTTACCGGACGCAGCGCCGGGAACAGAATAACGTCGCGGATGGTATgg                                                                  >  2:880924/1‑150 (MQ=255)
                                                                                     gCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATACCGTCGCGGATGGTATGGc                                                                 <  2:97458/150‑1 (MQ=255)
                                                                                       gACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTg                                                               <  1:930462/150‑1 (MQ=255)
                                                                                         cAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGtt                                                             >  2:1951745/1‑150 (MQ=255)
                                                                                          aaaTCGAGAAGCCAGAATACGATTATCTGCTATTACACAGCTTAAATAAAAAACCCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGGTg                                                            <  1:282660/150‑1 (MQ=255)
                                                                                           actCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTCAATAAAAACACCCCGACATTTGCCGGGGTTGTGAGCATACCGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTgg                                                           <  2:1256136/148‑1 (MQ=255)
                                                                                           aaTCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTgg                                                           <  1:1669373/150‑1 (MQ=255)
                                                                                           aaTCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTgg                                                           <  1:485236/150‑1 (MQ=255)
                                                                                               gcgAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGaa                                                       <  2:1682354/150‑1 (MQ=255)
                                                                                                 gAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACa                                                     <  2:2007376/150‑1 (MQ=255)
                                                                                                      cAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTTAACAGCAtt                                                >  2:1633574/1‑150 (MQ=255)
                                                                                                      cAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCAtt                                                >  2:2001803/1‑150 (MQ=255)
                                                                                                      cAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCAtt                                                >  2:1821707/1‑150 (MQ=255)
                                                                                                      cAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACACCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCAtt                                                >  2:225719/1‑150 (MQ=255)
                                                                                                       aGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTa                                               <  1:412345/150‑1 (MQ=255)
                                                                                                       aGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTa                                               <  2:801212/150‑1 (MQ=255)
                                                                                                          aaaaGATTATCTGCATTTACCAAGCTTAAATAAAAAAACCCCGTCATTTGGCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGTTGGTATGGCTGTTGGTGAACAGCATTACCa                                            <  1:1269631/150‑1 (MQ=255)
                                                                                                             cGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATAc                                         >  1:502751/1‑150 (MQ=255)
                                                                                                                    cTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGAt                                  >  2:1772776/1‑150 (MQ=255)
                                                                                                                    cTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGAt                                  <  1:1433389/150‑1 (MQ=255)
                                                                                                                       cATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATAcc                               >  1:1011214/1‑150 (MQ=255)
                                                                                                                           tACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAAtt                           >  1:1247575/1‑150 (MQ=255)
                                                                                                                           tACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACAAAtt                           >  1:1086941/1‑150 (MQ=255)
                                                                                                                              ccAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTccc                        <  1:935106/150‑1 (MQ=255)
                                                                                                                                  cTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGAc                    <  2:1921173/150‑1 (MQ=255)
                                                                                                                                  cTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGAc                    >  1:1938752/1‑150 (MQ=255)
                                                                                                                                     aaataaaAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACctg                 <  1:1472728/150‑1 (MQ=255)
                                                                                                                                      aataaaAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACctgc                >  2:1523611/1‑150 (MQ=255)
                                                                                                                                       ataaaAAAACCCCGACATATGCCGGGGTTGTGAGAATAACGAAAAGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGATGTTGGTGAACACCATTAAGAAACGGTCGATCCCAATTAGCAGAActgat               >  1:951549/1‑148 (MQ=255)
                                                                                                                                        taaaAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACctgctg              >  2:447142/1‑150 (MQ=255)
                                                                                                                                         ctcacaaCCCCGGCAAATGTCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGt             >  2:2018030/6‑150 (MQ=255)
                                                                                                                                         aaaaaaaCCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGt             <  2:1116762/150‑1 (MQ=255)
                                                                                                                                          aaaaaaCCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATAGCCGGGAACAGAATAACGTCGAGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACGTGCTGTc            >  2:251460/1‑150 (MQ=255)
                                                                                                                                           cacaaCCCCGGCAAATGTCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGACCAGCATTACCATACGGGCGATCCCAATTCCCACACCTGCTGTcg           >  2:1206600/4‑150 (MQ=255)
                                                                                                                                             caaCCCCGGCAAATGTCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTACCAGACCTGCTGTcggc         >  1:888926/2‑150 (MQ=255)
                                                                                                                                              aaCCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTcggcg        <  1:161941/150‑1 (MQ=255)
                                                                                                                                              aaCCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTcggcg        <  1:822667/150‑1 (MQ=255)
                                                                                                                                               aCCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTcggcgg       <  1:1538257/150‑1 (MQ=255)
                                                                                                                                               aCCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTcggcgg       <  1:1480565/150‑1 (MQ=255)
                                                                                                                                               aCCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTcggcgg       <  1:1913561/150‑1 (MQ=255)
                                                                                                                                                ccccGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGt      <  1:436765/150‑1 (MQ=255)
                                                                                                                                                ccccGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTAGATACCAATTCCCAGACCTGCTGTCGGCGGt      >  1:568626/1‑150 (MQ=255)
                                                                                                                                                  ccGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTaa    <  1:1555919/150‑1 (MQ=255)
                                                                                                                                                  ccGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCCGTaa    >  1:729247/1‑150 (MQ=255)
                                                                                                                                                    gACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGc  >  1:374866/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGATTATCAATGGTGTGATTTAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGC  >  NC_000913/3033475‑3033772

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