Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,232,484 A→G 6.9% intergenic (+30/+16) umuC → / ← dsbB translesion error‑prone DNA polymerase V subunit; DNA polymerase activity/oxidoreductase that catalyzes reoxidation of DsbA protein disulfide isomerase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,232,4840AG 6.9% 300.8 / 10.0 102intergenic (+30/+16)umuC/dsbBtranslesion error‑prone DNA polymerase V subunit; DNA polymerase activity/oxidoreductase that catalyzes reoxidation of DsbA protein disulfide isomerase I
Reads supporting (aligned to +/- strand):  ref base A (51/43);  new base G (5/2);  total (56/46)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.89e-01

Only 100 of 103 total aligned reads displayed.

GAGGGCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCA  >  NC_000913/1232335‑1232628
                                                                                                                                                     |                                                                                                                                                 
gagGGCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCa                                                                                                                                                   <  2:1866063/150‑1 (MQ=255)
  gggCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAgg                                                                                                                                                 >  2:758887/1‑150 (MQ=255)
    gCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAg                                                                                                                                               >  1:1204921/1‑150 (MQ=255)
     cAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAAGGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAACCGCGCTCCTGCGGGAgc                                                                                                                                              <  1:318999/150‑1 (MQ=255)
        agGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAATCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGcc                                                                                                                                           <  1:334158/150‑1 (MQ=255)
         ggAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCg                                                                                                                                          <  2:969598/150‑1 (MQ=255)
          gAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgcttt                                                                                                                                         <  1:990377/150‑4 (MQ=255)
            acacTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAAt                                                                                                                                       <  2:513668/150‑1 (MQ=255)
            acacTCTATTTTGCCGGGCAGGGGATCCAGAAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAAt                                                                                                                                       <  1:1932452/150‑1 (MQ=255)
               ctctATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCCGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGa                                                                                                                                    <  1:417244/150‑1 (MQ=255)
                tctATTTTGCCTGGCAGGGTATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTAAACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGAt                                                                                                                                   <  1:188407/150‑1 (MQ=255)
                     tttGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCg                                                                                                                              <  2:534697/150‑1 (MQ=255)
                       tGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGAc                                                                                                                            >  2:413883/1‑150 (MQ=255)
                        gCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGAcc                                                                                                                           >  1:670607/1‑150 (MQ=255)
                          cGGGCAGGGGATCCAGCAACATTGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGa                                                                                                                         <  1:439675/150‑1 (MQ=255)
                          cGGGCAGGGGATCCAGCAACAATGGCCGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGa                                                                                                                         <  1:272059/150‑1 (MQ=255)
                              cAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAg                                                                                                                     >  1:1552619/1‑150 (MQ=255)
                                ggggATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGAt                                                                                                                   >  2:2070892/1‑150 (MQ=255)
                                  ggATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGCCCGAACAGATcc                                                                                                                 >  1:1011981/1‑149 (MQ=255)
                                  ggATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCa                                                                                                                 >  1:248984/1‑150 (MQ=255)
                                     tCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGt                                                                                                              >  1:1827074/1‑150 (MQ=255)
                                      ccAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAACGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGtt                                                                                                             >  1:91980/1‑150 (MQ=255)
                                        aGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGtttt                                                                                                           >  1:1014717/1‑150 (MQ=255)
                                         gCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGttttt                                                                                                          <  2:126590/150‑1 (MQ=255)
                                          caacaaTGGCAGATGAAGCGAGCCATGCTATCACCACCTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAACTAAATATAGCGGAAGGAAAAAAGCGCTACCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGtttttt                                                                                                         <  1:137621/150‑1 (MQ=255)
                                           aacaaTGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTc                                                                                                        <  1:1521417/150‑1 (MQ=255)
                                           aacaaTGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTc                                                                                                        >  2:858431/1‑150 (MQ=255)
                                           aacaaTCGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGACGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTc                                                                                                        <  2:1737758/150‑1 (MQ=255)
                                             caaTGGTAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGGTTTTTCGc                                                                                                      >  2:600602/1‑150 (MQ=255)
                                             caaTGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGc                                                                                                      >  1:49314/1‑150 (MQ=255)
                                              aaTGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCt                                                                                                     <  1:157749/150‑1 (MQ=255)
                                              aaTGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCt                                                                                                     >  1:82755/1‑150 (MQ=255)
                                               aTGGCAGATGAATCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCtt                                                                                                    <  2:1196758/150‑1 (MQ=255)
                                               aTGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCtt                                                                                                    <  1:1298100/150‑1 (MQ=255)
                                               aTGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCtt                                                                                                    >  2:1672259/1‑150 (MQ=255)
                                                    aGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAAc                                                                                               <  1:1760612/150‑1 (MQ=255)
                                                     gATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACg                                                                                              <  1:2016812/150‑1 (MQ=255)
                                                        gAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATACCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTGGGTTTAACCGGCt                                                                                           >  2:207519/1‑150 (MQ=255)
                                                         aaGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTg                                                                                          >  1:280011/1‑150 (MQ=255)
                                                           gCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGCTCACGTTTTTTCGCTTTAAACGGCTggg                                                                                        >  1:771138/1‑150 (MQ=255)
                                                                ccATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATc                                                                                   <  2:948957/150‑1 (MQ=255)
                                                                 cATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATca                                                                                  <  2:1552619/150‑1 (MQ=255)
                                                                      tttCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCaccacc                                                                             <  1:1883765/150‑1 (MQ=255)
                                                                         caccacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGc                                                                          <  2:2048560/150‑1 (MQ=255)
                                                                         caccacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGc                                                                          >  1:1678214/1‑150 (MQ=255)
                                                                           ccacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCag                                                                        <  1:1281220/150‑2 (MQ=255)
                                                                           ccacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAc                                                                        >  2:1575614/1‑150 (MQ=255)
                                                                           ccaaGTTATACAACGCGAAGTTCTGATTTACTGAGTGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCTGGAGCGCCGAATGGAGTAGCGACCGAACAGATCCCGTGTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAc                                                                        <  1:1047295/150‑1 (MQ=255)
                                                                              ggTTATACAACGCGAAGTTCTGATTTTCTGAGGTTCTAATAAATATTGCGGCAGGAAAATAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGATTTTTCGCTTTAAACGGCTGGGAAATCACCTCCAGCACTGc                                                                    <  2:967099/149‑1 (MQ=255)
                                                                               gTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCCACCGAACAGATCACGTATTTTCGCTTTAAACGCCTCGGAAATCACCACCACCACTGCg                                                                   >  2:424969/1‑150 (MQ=255)
                                                                                ttATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGa                                                                  >  2:106371/1‑150 (MQ=255)
                                                                                 tataCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGAc                                                                 <  1:1313443/150‑1 (MQ=255)
                                                                                  ataCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAACGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGATGTAAACCGGGGGGAAATAACCACCAGCAATGCGACa                                                                >  2:376817/1‑150 (MQ=255)
                                                                                   taCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAACAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACaa                                                               >  1:1108093/1‑150 (MQ=255)
                                                                                   taCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACaa                                                               >  1:1299425/1‑150 (MQ=255)
                                                                                    aCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAAt                                                              >  1:739011/1‑150 (MQ=255)
                                                                                     cAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATc                                                             >  2:1364919/1‑150 (MQ=255)
                                                                                     cAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATc                                                             >  1:1733072/1‑150 (MQ=255)
                                                                                     cAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCAAGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATc                                                             >  1:202779/1‑150 (MQ=255)
                                                                                      aaCGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCa                                                            >  2:1175176/1‑150 (MQ=255)
                                                                                      aaCGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCa                                                            <  2:670607/150‑1 (MQ=255)
                                                                                       aCGCGAAGTTCTGATTTACTGAGTGTCAAATAATTATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAg                                                           <  1:665521/150‑1 (MQ=255)
                                                                                        cgcgAAGTTCTGATTTACTGAGGGTCAAATTATTATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGTAATTCACATCATCACGTGCGACAATAAgt                                                          >  2:231785/1‑149 (MQ=255)
                                                                                         gcgAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGt                                                         <  1:37900/150‑1 (MQ=255)
                                                                                          cgAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTa                                                        >  1:1595973/1‑150 (MQ=255)
                                                                                            aaGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAg                                                      >  2:1462281/1‑150 (MQ=255)
                                                                                                tCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGAt                                                  >  1:650114/1‑150 (MQ=255)
                                                                                                   gATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATaaa                                               >  2:277372/1‑150 (MQ=255)
                                                                                                      ttACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAAt                                            <  2:1643498/150‑1 (MQ=255)
                                                                                                         cTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGCGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATAcc                                         >  1:1384386/1‑150 (MQ=255)
                                                                                                          tGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCg                                        <  1:1582840/150‑1 (MQ=255)
                                                                                                          agAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGAGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACGACCAGCACTGCGACAATCAGGTACGCGATAAAAATACCg                                        >  1:91227/2‑150 (MQ=255)
                                                                                                           gtgggTCAAATAAATATAGCGGCAGGAAAAAAGCGCTGCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGa                                       <  1:718818/148‑1 (MQ=255)
                                                                                                            aGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGag                                      <  2:1013153/150‑1 (MQ=255)
                                                                                                              ggTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGagca                                    >  2:559440/1‑150 (MQ=255)
                                                                                                                 cAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGcc                                 <  1:288116/150‑1 (MQ=255)
                                                                                                                 cAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGcc                                 <  1:2095567/150‑1 (MQ=255)
                                                                                                                   aataaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCAc                               >  2:89522/1‑150 (MQ=255)
                                                                                                                   aataaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCAc                               >  1:455154/1‑150 (MQ=255)
                                                                                                                   aataaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCAc                               >  2:1638143/1‑150 (MQ=255)
                                                                                                                         tataGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGc                         >  1:1680491/1‑150 (MQ=255)
                                                                                                                           taGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGCCCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCCGCAt                       >  1:1762936/1‑150 (MQ=255)
                                                                                                                              cGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTc                    <  1:1427018/150‑1 (MQ=255)
                                                                                                                               ggCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTcc                   <  1:1462281/150‑1 (MQ=255)
                                                                                                                               ggCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTcc                   <  2:1827074/150‑1 (MQ=255)
                                                                                                                                gCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCa                  <  2:1204921/150‑1 (MQ=255)
                                                                                                                                 cAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAg                 >  1:1138017/1‑150 (MQ=255)
                                                                                                                                 cAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAg                 >  2:1545224/1‑150 (MQ=255)
                                                                                                                                 cAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAg                 <  1:1175176/150‑1 (MQ=255)
                                                                                                                                    gAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGAcc              <  1:1209445/150‑1 (MQ=255)
                                                                                                                                     cattcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCGTTTCCAGACCt             >  1:1502096/7‑150 (MQ=255)
                                                                                                                                     cattcggcgcTCCTGCGGGAGCGCCGAAGGGATTAGCGCCCGAACAGCTCCCGGTTTTTCGCGGTACACGGCTGGGAAATCACCCCCCGCACTGCGACAAGCAGGGAAGCGATAAAAATACCGAGCAGCCGCTGCGGCATTTCCAGAccg             >  1:1510721/7‑149 (MQ=255)
                                                                                                                                      attcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGAGAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTa            >  2:1758406/6‑150 (MQ=255)
                                                                                                                                        tcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTaaa          >  1:626131/4‑150 (MQ=255)
                                                                                                                                         cggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGAGCACGTTTTTTCGCGTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCCCTGCGGCAGTGCAAGCCCTaaaa         >  2:1448665/3‑150 (MQ=255)
                                                                                                                                            cgcTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCCCGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAAt      <  1:106371/150‑1 (MQ=255)
                                                                                                                                            cgcTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAAt      >  2:152650/1‑150 (MQ=255)
                                                                                                                                             gcTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATc     <  1:960764/150‑1 (MQ=255)
                                                                                                                                                cccGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATATAAATCCCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAACTcaca  >  1:1088124/1‑150 (MQ=255)
                                                                                                                                                cccGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATACCa  >  2:508520/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GAGGGCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCA  >  NC_000913/1232335‑1232628

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: