Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,041,940 C→A 5.5% intergenic (+25/+90) appA → / ← etk periplasmic phosphoanhydride phosphatase/multiple inositol‑polyphosphate phosphatase/protein‑tyrosine kinase Etk

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,041,9400CA 5.5% 309.9 / 5.6 91intergenic (+25/+90)appA/etkperiplasmic phosphoanhydride phosphatase/multiple inositol‑polyphosphate phosphatase/protein‑tyrosine kinase Etk
Reads supporting (aligned to +/- strand):  ref base C (49/37);  new base A (3/2);  total (52/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

CGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCT  >  NC_000913/1041791‑1042089
                                                                                                                                                     |                                                                                                                                                     
cgccgcCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTAc                                                                                                                                                       >  2:693493/1‑150 (MQ=255)
cgccgcCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTAc                                                                                                                                                       >  2:60565/1‑150 (MQ=255)
 gccgccCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTAcc                                                                                                                                                      <  2:200080/150‑1 (MQ=255)
    gccCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGa                                                                                                                                                   >  2:130457/1‑150 (MQ=255)
    gccCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGa                                                                                                                                                   <  2:786112/150‑1 (MQ=255)
        gtagagGTGAAACTGACCGTGGCAGGATGTGAAGAGGGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGGAGTTTGTAATGCATAAAAAAGAGCATTGAGTTACCTGAATGc                                                                                                                                               <  2:709893/148‑1 (MQ=255)
        ggAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGc                                                                                                                                               <  1:681392/150‑1 (MQ=255)
           gagGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGGTAACTGAATGctct                                                                                                                                            <  1:1523137/150‑1 (MQ=255)
             ggTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGa                                                                                                                                          >  1:1989501/1‑150 (MQ=255)
              gTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAg                                                                                                                                         >  1:1586576/1‑150 (MQ=255)
               tGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAgg                                                                                                                                        >  1:1892355/1‑150 (MQ=255)
                gAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGc                                                                                                                                       >  1:1361615/1‑150 (MQ=255)
                gAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGc                                                                                                                                       <  2:691376/150‑1 (MQ=255)
                 aaaCTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCt                                                                                                                                      >  2:1098051/1‑150 (MQ=255)
                 aaaCTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGGTAACTGAATGctctttttta                                                                                                                                      <  1:1948608/150‑7 (MQ=255)
                  aaCTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCtg                                                                                                                                     >  1:1227390/1‑150 (MQ=255)
                     tGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCtgatg                                                                                                                                  <  1:591024/150‑1 (MQ=255)
                      gACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGGGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCtgatga                                                                                                                                 <  1:80894/150‑1 (MQ=255)
                        cccTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACa                                                                                                                               <  1:1970077/150‑1 (MQ=255)
                           tGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAAc                                                                                                                            >  1:703167/1‑150 (MQ=255)
                                ggATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACgaaga                                                                                                                       <  2:1643084/150‑1 (MQ=255)
                                ggATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACgaaga                                                                                                                       >  1:1496974/1‑150 (MQ=255)
                                   tgtgAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACt                                                                                                                    >  2:558053/1‑150 (MQ=255)
                                     tgAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGAATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGt                                                                                                                  >  2:1874764/1‑150 (MQ=255)
                                      gAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTc                                                                                                                 <  2:25270/150‑1 (MQ=255)
                                      gAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTc                                                                                                                 >  1:157935/1‑150 (MQ=255)
                                      gAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTc                                                                                                                 >  1:1009212/1‑150 (MQ=255)
                                         gagCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTaa                                                                                                              <  2:1227390/150‑1 (MQ=255)
                                          agCGTAATTCGCAGGGGATGTGTTCGTTGGCATGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAAt                                                                                                             <  1:1955992/150‑1 (MQ=255)
                                           gCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATg                                                                                                            >  2:314931/1‑150 (MQ=255)
                                           gCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACACATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATg                                                                                                            <  2:790124/150‑1 (MQ=255)
                                             gAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCg                                                                                                          >  2:1722283/1‑150 (MQ=255)
                                             gAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCg                                                                                                          >  1:1708248/1‑150 (MQ=255)
                                              aaaTGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGt                                                                                                         <  1:314931/150‑1 (MQ=255)
                                                 tGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGa                                                                                                      >  2:885266/1‑150 (MQ=255)
                                                 tGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGa                                                                                                      <  1:1450132/150‑1 (MQ=255)
                                                  gcgcAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGAc                                                                                                     >  1:779996/1‑150 (MQ=255)
                                                          cATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAg                                                                                             >  1:27633/1‑150 (MQ=255)
                                                          cATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAg                                                                                             >  1:237813/1‑150 (MQ=255)
                                                           aTGTGTTCGTTGGCAGGTTTTACTCAAATCGTGAATTAAGCAGGCATACCGGCGTGCAGTTTGTAATGAATAAAAAATACCATTGAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAgg                                                                                            <  1:562738/150‑1 (MQ=255)
                                                            tgtgTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGc                                                                                           >  1:1243710/1‑150 (MQ=255)
                                                             gtgtTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCg                                                                                          >  1:1395021/1‑150 (MQ=255)
                                                              tgtTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGt                                                                                         >  1:349763/1‑150 (MQ=255)
                                                              tgtTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGt                                                                                         >  1:1711161/1‑150 (MQ=255)
                                                               gtTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGtt                                                                                        <  1:130457/150‑1 (MQ=255)
                                                                  cGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCAc                                                                                     >  2:1579421/1‑150 (MQ=255)
                                                                   gTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGGCCGGAAAAGGCGTTCAcg                                                                                    >  1:1169479/1‑150 (MQ=255)
                                                                   gTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCAcg                                                                                    <  2:1361615/150‑1 (MQ=255)
                                                                      ggCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCAcgccg                                                                                 <  2:1601260/150‑1 (MQ=255)
                                                                          ggTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATc                                                                             >  1:1350072/1‑150 (MQ=255)
                                                                                aCGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCa                                                                       <  1:210649/150‑1 (MQ=255)
                                                                                   cAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACtt                                                                    >  2:798291/1‑150 (MQ=255)
                                                                                     aaTCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTc                                                                  <  1:558053/150‑1 (MQ=255)
                                                                                      aTCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCa                                                                 >  2:1289011/1‑150 (MQ=255)
                                                                                          tgaatgaaGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGcc                                                                        <  1:680238/139‑1 (MQ=255)
                                                                                          tgaatgaaGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGcc                                                                        >  2:680238/1‑139 (MQ=255)
                                                                                          tgaatgaaGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGttt                                                             <  2:703167/150‑1 (MQ=255)
                                                                                           gaatgaaGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGtttt                                                            <  1:1722283/150‑1 (MQ=255)
                                                                                           gaatgaaGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTCTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGtttt                                                            <  2:1989501/150‑1 (MQ=255)
                                                                                             atgaaGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTAc                                                          <  1:1032875/150‑1 (MQ=255)
                                                                                             atgaaGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTAc                                                          >  2:1235863/1‑150 (MQ=255)
                                                                                                  gCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCttt                                                     >  1:1446884/1‑150 (MQ=255)
                                                                                                     cGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAAATGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTtctc                                                  <  1:1098051/150‑1 (MQ=255)
                                                                                                       cATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGTAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCgg                                                <  2:1708248/150‑1 (MQ=255)
                                                                                                       cATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCgg                                                >  2:1216963/1‑150 (MQ=255)
                                                                                                         tACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAg                                              <  1:1285100/150‑1 (MQ=255)
                                                                                                           ccGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTa                                            >  2:1851988/1‑150 (MQ=255)
                                                                                                              gCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACt                                         <  2:1446884/150‑1 (MQ=255)
                                                                                                                     gTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCg                                  <  2:1243710/150‑1 (MQ=255)
                                                                                                                     ctttGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCg                                  >  1:1234808/2‑150 (MQ=255)
                                                                                                                      tttGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGt                                 >  2:151293/1‑150 (MQ=255)
                                                                                                                      tttGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGt                                 <  2:1586576/150‑1 (MQ=255)
                                                                                                                      tttGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGt                                 <  1:688961/150‑1 (MQ=255)
                                                                                                                       ttGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTa                                <  1:1289011/150‑1 (MQ=255)
                                                                                                                          tAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATa                             >  1:498116/1‑150 (MQ=255)
                                                                                                                           aaTGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAg                            <  2:1350072/150‑1 (MQ=255)
                                                                                                                             tGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGtt                          <  2:370884/150‑1 (MQ=255)
                                                                                                                               cATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTtat                        >  1:456731/1‑150 (MQ=255)
                                                                                                                                aTAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTtata                       >  1:407915/1‑150 (MQ=255)
                                                                                                                                aTAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTtata                       >  1:1851684/1‑150 (MQ=255)
                                                                                                                                aTAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTtata                       >  1:1352816/1‑150 (MQ=255)
                                                                                                                                  gcctcagagCATTCAGGTAACTAAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGc                     >  1:515782/6‑150 (MQ=255)
                                                                                                                                  aaaaaaGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGc                     <  2:27633/150‑1 (MQ=255)
                                                                                                                                   aaaaaGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGcc                    >  1:878675/1‑150 (MQ=255)
                                                                                                                                    ctcagagCATTCAGGTAACTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCg                   >  1:199972/4‑150 (MQ=255)
                                                                                                                                    aaaaGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCg                   >  1:1791681/1‑150 (MQ=255)
                                                                                                                                    aaaaGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCg                   >  1:16788/1‑150 (MQ=255)
                                                                                                                                     tcagagCATTCAGGTAACTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGt                  >  2:1201603/3‑150 (MQ=255)
                                                                                                                                     aaaGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGt                  >  2:710950/1‑150 (MQ=255)
                                                                                                                                            aTTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTa           <  2:1234808/150‑1 (MQ=255)
                                                                                                                                               cAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGc        >  2:1807521/1‑150 (MQ=255)
                                                                                                                                               cAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGc        >  1:1603915/1‑150 (MQ=255)
                                                                                                                                                 gTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGt     >  2:1170551/1‑150 (MQ=255)
                                                                                                                                                    aCCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGc   >  2:1839767/1‑150 (MQ=255)
                                                                                                                                                     ccTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTGCTCTGTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTACCTGTAAGCGGTGCt  >  2:825981/1‑150 (MQ=255)
                                                                                                                                                     ccTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCt  <  1:885266/150‑1 (MQ=255)
                                                                                                                                                     ccTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCt  >  1:1038943/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCT  >  NC_000913/1041791‑1042089

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: