Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,232,481 C→T 5.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 5.8% 232.2 / 3.2 85intergenic (+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 (48/32);  new base T (3/2);  total (51/34)
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.50e-01

GGCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCA  >  NC_000913/1232338‑1232628
                                                                                                                                               |                                                                                                                                                   
ggCAGAGGAGCACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCGGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAgg                                                                                                                                                >  2:402173/1‑149 (MQ=255)
   agagGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAgc                                                                                                                                             <  2:342527/149‑1 (MQ=255)
     agGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAgcgc                                                                                                                                           <  2:319616/149‑1 (MQ=255)
      gtaaCACGCTATTTTTCCGGGCTGGGGATCCAGCAAGAATTTAAGATGAAGCGAGCCATGCTTTCACCACGTTATACATCGCGTCGTTCTGATTTACTGAGGGTCAAATAAATATAGCTGCAGGAAAACAGCGCTCCCGCAGGAGCGcc                                                                                                                                          <  2:679317/147‑1 (MQ=255)
        aaCACTCTCTATTGCCGGCCAGGGGATCCAGCAACAATGGCAGATGAACCGAGCCATGCTATCACCACGTTATACAACGCGAAGTTCTGAATTACTGAGGGGCAAATAAATCAAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgcttt                                                                                                                                        <  1:766879/149‑4 (MQ=255)
         acacTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGaa                                                                                                                                       <  1:788075/149‑1 (MQ=255)
          cacTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAAt                                                                                                                                      <  1:526566/149‑1 (MQ=255)
           acTCTATTTTGCCGGTCAGTGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgctttttt                                                                                                                                     <  2:1479651/149‑7 (MQ=255)
              ctATTTTGCCGGGCAGGGGATCCAGCAACAATGGAAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGAt                                                                                                                                  <  2:217576/149‑1 (MQ=255)
               tATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGAtt                                                                                                                                 >  1:1515701/1‑149 (MQ=255)
                aTTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTa                                                                                                                                >  2:140538/1‑149 (MQ=255)
                  tttGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGc                                                                                                                              <  1:1007162/149‑1 (MQ=255)
                      ccGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGTGCCATTCTTTCACCACGTTATACAACGCGATGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCTGGAAAAAAGCGCTACCGCAGGAGCGCCGAATGTTTTAGCGAcc                                                                                                                          <  2:7750/149‑1 (MQ=255)
                         ggCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATGTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGaa                                                                                                                       <  1:1077230/149‑1 (MQ=255)
                         ggCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGATTTCACCACGTTATACAACGCGAAGTTATGATTTACTGAGTGTCAAATAAATATAGCGGCTGGAAAAAATCGCTCCCGCAGGAGAGCCGAATGGAGTAGAGACCGaa                                                                                                                       >  2:918256/1‑149 (MQ=255)
                          gCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAAc                                                                                                                      >  2:345416/1‑149 (MQ=255)
                           cAGGGGATCCAGCAACAATGGCAGATGCAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACa                                                                                                                     <  2:1366574/149‑1 (MQ=255)
                           cAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGCCCGAACAg                                                                                                                    >  2:1416215/1‑149 (MQ=255)
                            aGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAACTATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAg                                                                                                                    >  1:1271606/1‑149 (MQ=255)
                             ggggATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCTGCAGGAAAACAGCGCTCCCTCAGGAGCGCCGAATGGATTAGCGACCGAACAGa                                                                                                                   >  1:1059066/1‑149 (MQ=255)
                             ggggATCCAGCAACAATGGCAGATGAAGCGAGACATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCCCCGAATGGAATAGCGACCGAACAGa                                                                                                                   >  1:564900/1‑149 (MQ=255)
                               ggATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATc                                                                                                                 <  1:1091028/149‑1 (MQ=255)
                                    cAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGtt                                                                                                            >  1:1330887/1‑149 (MQ=255)
                                     aGCAACAATGGCAGATGAAGCGAGCCAAGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGttt                                                                                                           <  2:584015/149‑1 (MQ=255)
                                       caacaaTGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGttttt                                                                                                         <  2:232172/149‑1 (MQ=255)
                                        aacaaTGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGtttttt                                                                                                        >  2:1088952/1‑149 (MQ=255)
                                        aacaaTGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGtttttt                                                                                                        <  1:347403/149‑1 (MQ=255)
                                        aacaaTGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGtttttt                                                                                                        >  1:249818/1‑149 (MQ=255)
                                        aacaaTGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGACGAATGGATTAGCGACCGAACAGATCACGtttttt                                                                                                        >  2:451292/1‑149 (MQ=255)
                                            aTGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACTACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCt                                                                                                    <  1:108035/149‑1 (MQ=255)
                                             tGGCAGATGAAGCTAGCAATGCTTTCACCACGTTATAAAACGCGAAGTTCTGATTTACTGTGGGTCAAATAAATATAGCGGCAGGAAAAACGCGCTCCCGCAGGAGCGCCGAATGTATTAGCGACCGAACAGATCACGTTTTTTCGCtt                                                                                                   <  1:592611/149‑1 (MQ=255)
                                                   aTGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACg                                                                                             >  2:1111916/1‑149 (MQ=255)
                                                     gAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTTATTTACTGAGGGTCAACTAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGc                                                                                           <  2:280858/149‑1 (MQ=255)
                                                        gCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTAGCTTTAAACGGATgg                                                                                        >  1:1267216/1‑149 (MQ=255)
                                                             ccATGCTTTCACCACGTTATACATGGCGAAGTTCTGATTTACGGAGGGTCAAATAAATTTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATAACGTTTTTTCGCTTTAAACGGCTGGGAAAt                                                                                   >  1:771515/1‑149 (MQ=255)
                                                             ccATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAAt                                                                                   >  1:715163/1‑149 (MQ=255)
                                                              cATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATc                                                                                  <  1:615446/149‑1 (MQ=255)
                                                               aTGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAATTATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTGTTCGCTTTAAACGGCTGGGAAATca                                                                                 >  2:1489720/1‑149 (MQ=255)
                                                               aTGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGCCCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATcc                                                                                 >  2:1357893/1‑148 (MQ=255)
                                                               aTGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATca                                                                                 >  2:1216964/1‑149 (MQ=255)
                                                                   tttCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAATTAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGATCAGCTCACGTTTTTTCGCTTTAAACGGCTGGGAAATcaccac                                                                             >  1:1230790/1‑149 (MQ=255)
                                                                     tCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGCGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCAccacca                                                                           >  2:1347998/1‑149 (MQ=255)
                                                                      caccacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAg                                                                          <  1:622784/149‑1 (MQ=255)
                                                                        ccacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCa                                                                        >  2:295902/1‑149 (MQ=255)
                                                                         cacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAc                                                                       >  2:596865/1‑149 (MQ=255)
                                                                                 aCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATAAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACaa                                                               >  1:314428/1‑149 (MQ=255)
                                                                                   aaCGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATc                                                             >  2:405378/1‑149 (MQ=255)
                                                                                   aaCGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATc                                                             >  1:335424/1‑149 (MQ=255)
                                                                                        gAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATAACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCATCACTGCGACAATCAGGTa                                                        >  2:867137/1‑149 (MQ=255)
                                                                                          aGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGGATGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACCGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAg                                                      <  1:735529/149‑1 (MQ=255)
                                                                                          aGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAg                                                      <  2:715163/149‑1 (MQ=255)
                                                                                           tttCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGc                                                     <  1:517124/148‑1 (MQ=255)
                                                                                              cTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGAt                                                  >  1:785016/1‑149 (MQ=255)
                                                                                               tGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGCACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATa                                                 >  2:983755/1‑149 (MQ=255)
                                                                                                gATTTACTGAGGGGCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATaa                                                <  2:473165/149‑1 (MQ=255)
                                                                                                 aTTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATaaa                                               >  1:672142/1‑149 (MQ=255)
                                                                                                      cTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCAGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATaa                                          >  2:507102/1‑148 (MQ=255)
                                                                                                       tGAGTGTCAAATAAATATAGCGGCAGGAAAAATGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACATATCACGTTTTTTCGCTTTAAACGGCTGGGATATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATAcc                                         <  1:453826/149‑1 (MQ=255)
                                                                                                            gTCAAATAAATATAGCGGCACGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTACGCTTTAAACCGCTGGGAAATCACCACCAGCACAGCGACAATCAGGTAAGCGATAAAAATACCGaaga                                    >  2:1429386/1‑146 (MQ=255)
                                                                                                             tCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGagcag                                   >  1:1400437/1‑149 (MQ=255)
                                                                                                                 ataaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCTATAAAAATACCGAGCAGCCAc                               >  2:392507/1‑149 (MQ=255)
                                                                                                                 ataaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCAc                               >  1:93723/1‑149 (MQ=255)
                                                                                                                 ataaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCAc                               >  2:430725/1‑149 (MQ=255)
                                                                                                                 ataaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGAGTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAAGCAGGTAAGCGATAAAAATCCAGCGCAGCCAc                               >  2:868254/1‑149 (MQ=255)
                                                                                                                     atatAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCg                           <  1:557516/149‑1 (MQ=255)
                                                                                                                      tctaGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTGCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCgg                          <  1:693968/147‑1 (MQ=255)
                                                                                                                        taGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGTAATCGCCACCAGCACTGCGACAAGCAGGTAAGCGATATAAATACCGAGCAGCCACTGCGGct                        >  2:1365754/1‑148 (MQ=255)
                                                                                                                             gCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGTGAAATCACCACCAGCACTGAGACACTCAGGTAAGCGATAAAAAAACCGAGCAGCCACTGCGGCATTTca                   >  2:899122/1‑148 (MQ=255)
                                                                                                                             gCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTcc                   >  2:57371/1‑149 (MQ=255)
                                                                                                                               aGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGTCCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCAGCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAg                 <  1:985667/149‑1 (MQ=255)
                                                                                                                               aGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTc                    >  1:1531187/1‑146 (MQ=255)
                                                                                                                               aGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTc                    <  2:1531187/146‑1 (MQ=255)
                                                                                                                               aGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGGTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAg                 <  2:1271606/149‑1 (MQ=255)
                                                                                                                                ggAAAACCGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGCCCGAACAGATCACGTTTATTCGCTATAAACGGCTGGGAACTCACCACCAGCACTGCGACAATCAGGTAAGGGATAATAATACGGAGCAGCCACTGCGGCATTTCCAGa                >  1:99599/1‑149 (MQ=255)
                                                                                                                                 gAAAACAGCGCTCCCGCAGGACAGCCGAATGGATTAGCGCCCGATCTGATCACGTTTTTTCGCTTTAAACGGCTGGGAATTCACAACCAGCACTGCGTCAAAGAGGTAAGCGATAAAAATACCGAGCAGCCAATGCGGCATTTCCAGAc               >  2:923569/1‑149 (MQ=255)
                                                                                                                                 gAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGCTCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAGTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGAc               >  1:352648/1‑149 (MQ=255)
                                                                                                                                 gAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGAc               >  2:1374951/1‑149 (MQ=255)
                                                                                                                                  aaaaaaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGCCCGAACAGATCACGTTTTTCCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCCCTGCGGCATTTCCAGAcc              >  2:1467527/1‑149 (MQ=255)
                                                                                                                                  aaaaaaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGAcc              >  1:1460482/1‑149 (MQ=255)
                                                                                                                                  aaaaaaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGAcc              <  2:33553/149‑1 (MQ=255)
                                                                                                                                    ttcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAATTCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAGATACCGAGCAGCCACTGCGGCATTTCCAGACCTa            >  1:112312/5‑149 (MQ=255)
                                                                                                                                       ggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTaaaa         >  2:777323/2‑149 (MQ=255)
                                                                                                                                       ggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTaaaa         >  2:433173/2‑149 (MQ=255)
                                                                                                                                         cgcTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTaaaaaa       <  1:140538/149‑1 (MQ=255)
                                                                                                                                            tCCCGCAGGAGCGCCGAATGGATTAGCGACCTAACTGATCAGGTTTTTTCGCTTTAAACGGCAGGGAAATCACCAACAGCACTGCGACAATCATGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATcc    <  1:150983/149‑1 (MQ=255)
                                                                                                                                              ccGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCa  >  1:323134/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GGCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCA  >  NC_000913/1232338‑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: