Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,232,481 C→T 6.8% intergenic (+27/+19) 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,4810CT 6.8% 316.5 / 8.4 102intergenic (+27/+19)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 C (50/45);  new base T (5/2);  total (55/47)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.83e-01

Only 100 of 104 total aligned reads displayed.

AGAGGGCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCA  >  NC_000913/1232334‑1232628
                                                                                                                                                   |                                                                                                                                                    
agagGGCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGATAAACGCGCTCCCGc                                                                                                                                                    >  2:101541/1‑150 (MQ=255)
 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:82755/1‑150 (MQ=255)
                                               aaTGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCt                                                                                                     <  1:157749/150‑1 (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)
                                                                            ccacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCag                                                                        <  1:1281220/150‑2 (MQ=255)
                                                                            ccacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAc                                                                        >  2:1575614/1‑150 (MQ=255)
                                                                            ccaaGTTATACAACGCGAAGTTCTGATTTACTGAGTGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCTGGAGCGCCGAATGGAGTAGCGACCGAACAGATCCCGTGTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAc                                                                        <  1:1047295/150‑1 (MQ=255)
                                                                               ggTTATACAACGCGAAGTTCTGATTTTCTGAGGTTCTAATAAATATTGCGGCAGGAAAATAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGATTTTTCGCTTTAAACGGCTGGGAAATCACCTCCAGCAC‑TGc                                                                    <  2:967099/149‑1 (MQ=255)
                                                                                gTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCCACCGAACAGATCACGTATTTTCGCTTTAAACGCCTCGGAAATCACCACCACCAC‑TGCg                                                                   >  2:424969/1‑150 (MQ=255)
                                                                                 ttATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGa                                                                  >  2:106371/1‑150 (MQ=255)
                                                                                  tataCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGAc                                                                 <  1:1313443/150‑1 (MQ=255)
                                                                                   ataCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAACGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGATGTAAACCGGGGGGAAATAACCACCAGCAA‑TGCGACa                                                                >  2:376817/1‑150 (MQ=255)
                                                                                    taCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAACAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACaa                                                               >  1:1108093/1‑150 (MQ=255)
                                                                                    taCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACaa                                                               >  1:1299425/1‑150 (MQ=255)
                                                                                     aCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAAt                                                              >  1:739011/1‑150 (MQ=255)
                                                                                      cAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCTCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAACCGGCTGGGAAATCACCACCAGCAC‑TGCTACAAta                                                             >  2:1374638/1‑149 (MQ=255)
                                                                                      cAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATc                                                             >  1:1733072/1‑150 (MQ=255)
                                                                                      cAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATc                                                             >  2:1364919/1‑150 (MQ=255)
                                                                                      cAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCAAGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATc                                                             >  1:202779/1‑150 (MQ=255)
                                                                                       aaCGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCa                                                            <  2:670607/150‑1 (MQ=255)
                                                                                       aaCGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCa                                                            >  2:1175176/1‑150 (MQ=255)
                                                                                         cgcgAAGTTCTGATTTACTGAGGGTCAAATTATTATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGTAATTCA‑CATCATCACGTGCGACAATAAgt                                                          >  2:231785/1‑149 (MQ=255)
                                                                                          gcgAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGt                                                         <  1:37900/150‑1 (MQ=255)
                                                                                           cgAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTa                                                        >  1:1595973/1‑150 (MQ=255)
                                                                                             aaGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAg                                                      >  2:1462281/1‑150 (MQ=255)
                                                                                                 tCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGAt                                                  >  1:650114/1‑150 (MQ=255)
                                                                                                    gATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGATaaa                                               >  2:277372/1‑150 (MQ=255)
                                                                                                          cTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGCGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGATAAAAATAcc                                         >  1:1384386/1‑150 (MQ=255)
                                                                                                           tGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGATAAAAATACCg                                        <  1:1582840/150‑1 (MQ=255)
                                                                                                           agAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGAGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACGACCAGCAC‑TGCGACAATCAGGTACGCGATAAAAATACCg                                        >  1:91227/2‑150 (MQ=255)
                                                                                                            gtgggTCAAATAAATATAGCGGCAGGAAAAAAGCGCTGCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGATAAAAATACCGa                                       <  1:718818/148‑1 (MQ=255)
                                                                                                             aGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGATAAAAATACCGag                                      <  2:1013153/150‑1 (MQ=255)
                                                                                                               ggTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGATAAAAATACCGagca                                    >  2:559440/1‑150 (MQ=255)
                                                                                                               ggTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGATAAAAATACCGagca                                    <  1:1982337/150‑1 (MQ=255)
                                                                                                                  cAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGcc                                 <  1:288116/150‑1 (MQ=255)
                                                                                                                  cAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGcc                                 <  1:2095567/150‑1 (MQ=255)
                                                                                                                    aataaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCAc                               >  2:89522/1‑150 (MQ=255)
                                                                                                                    aataaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCAc                               >  1:455154/1‑150 (MQ=255)
                                                                                                                    aataaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCAc                               >  2:1638143/1‑150 (MQ=255)
                                                                                                                          tataGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGc                         >  1:1680491/1‑150 (MQ=255)
                                                                                                                            taGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGCCCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCCGCAt                       >  1:1762936/1‑150 (MQ=255)
                                                                                                                               cGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTc                    <  1:1427018/150‑1 (MQ=255)
                                                                                                                                ggCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTcc                   <  2:1827074/150‑1 (MQ=255)
                                                                                                                                ggCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTcc                   <  1:1462281/150‑1 (MQ=255)
                                                                                                                                 gCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCa                  <  2:1204921/150‑1 (MQ=255)
                                                                                                                                  cAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAg                 >  2:1545224/1‑150 (MQ=255)
                                                                                                                                  cAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAg                 <  1:1175176/150‑1 (MQ=255)
                                                                                                                                  cAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAg                 >  1:1138017/1‑150 (MQ=255)
                                                                                                                                     gAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGAcc              <  1:1209445/150‑1 (MQ=255)
                                                                                                                                      cattcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCGTTTCCAGACCt             >  1:1502096/7‑150 (MQ=255)
                                                                                                                                      cattcggcgcTCCTGCGGGAGCGCCGAAGGGATTAGCGCCCGAACAGCTCCCGGTTTTTCGCGGTACACGGCTGGGAAATCACCCCCCGCAC‑TGCGACAAGCAGGGAAGCGATAAAAATACCGAGCAGCCGCTGCGGCATTTCCAGAccg             >  1:1510721/7‑149 (MQ=255)
                                                                                                                                       attcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGAGAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTa            >  2:1758406/6‑150 (MQ=255)
                                                                                                                                         tcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTaaa          >  1:626131/4‑150 (MQ=255)
                                                                                                                                          cggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGAGCACGTTTTTTCGCGTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCCCTGCGGCAGTGCAAGCCCTaaaa         >  2:1448665/3‑150 (MQ=255)
                                                                                                                                             cgcTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCCCGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAAt      <  1:106371/150‑1 (MQ=255)
                                                                                                                                             cgcTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAAt      >  2:152650/1‑150 (MQ=255)
                                                                                                                                              gcTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATc     <  1:960764/150‑1 (MQ=255)
                                                                                                                                                 cccGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGATATAAATCCCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAACTcaca  >  1:1088124/1‑150 (MQ=255)
                                                                                                                                                 cccGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATACCa  >  2:508520/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGAGGGCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAC‑TGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCA  >  NC_000913/1232334‑1232628

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