Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,020,511 C→G 5.2% intergenic (+22/‑29) ssnA → / → ygfM putative chlorohydrolase/aminohydrolase/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,020,5110CG 5.2% 405.0 / 9.4 115intergenic (+22/‑29)ssnA/ygfMputative chlorohydrolase/aminohydrolase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (53/56);  new base G (3/3);  total (56/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 115 total aligned reads displayed.

GCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAACCC  >  NC_000913/3020365‑3020657
                                                                                                                                                   |                                                                                                                                                  
gCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTGACTTCGATTGCGACTCCAGTTATGCACAAGCCAGAAAAGCCGATGCCAGTATGTGGCGTCGGATGGATGCGGTGGCATAAATGACAGATGCCCTCTTCCC‑GCAg                                                                                                                                                  >  3:802613/1‑149 (MQ=255)
gCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCa                                                                                                                                                   >  1:492369/1‑148 (MQ=255)
gCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAg                                                                                                                                                  <  4:837804/149‑1 (MQ=255)
    gagggTCAATAGTGTGATGGTCTATGAAGAACGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACATGACAGAGAAGCCGCTGCCAGTATGTGGTGTCGCATCGATGCGCTGGCATACATGACAGATGCCCTCTTCCC‑GCAGGGaa                                                                                                                                              <  3:307756/146‑1 (MQ=255)
    gATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGa                                                                                                                                               <  1:644585/148‑1 (MQ=255)
    gATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGa                                                                                                                                               >  1:675008/1‑148 (MQ=255)
    gATGGTCAATGGTGTGATGGTCTATGAAGACCATCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCCGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGaa                                                                                                                                              <  4:777195/149‑1 (MQ=255)
       ggTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAagag                                                                                                                                           <  4:2506/149‑1 (MQ=255)
       ggTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAagag                                                                                                                                           <  3:416187/149‑1 (MQ=255)
       ggTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATACTCTCTTCCC‑TGCGGGAagag                                                                                                                                           <  4:602216/149‑1 (MQ=255)
         tCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAggg                                                                                                                                         <  4:616040/149‑1 (MQ=255)
          cAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCTGC‑GGGAagag                                                                                                                                           <  2:202732/146‑1 (MQ=255)
          cAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑TGCGGGAagag                                                                                                                                           >  1:202732/1‑146 (MQ=255)
          cAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAggg                                                                                                                                         >  2:713874/1‑148 (MQ=255)
          cAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAggg                                                                                                                                         <  1:689361/148‑1 (MQ=255)
           aaTGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGtt                                                                                                                                       >  4:622475/1‑149 (MQ=255)
            aTGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGtt                                                                                                                                       >  2:668946/1‑148 (MQ=255)
              ggTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTaaa                                                                                                                                    <  4:708997/149‑1 (MQ=255)
               gtgtGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGAAGCGCTGGCATAAATGACAGATGCCCTCTGCCC‑GCAGGGAAGAGGGTTAAAt                                                                                                                                   <  4:802613/149‑1 (MQ=255)
                 gtgATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATa                                                                                                                                  <  2:492369/148‑1 (MQ=255)
                 gtgATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATa                                                                                                                                  <  1:567132/148‑1 (MQ=255)
                     tGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGa                                                                                                                              <  1:668946/148‑1 (MQ=255)
                          tATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAaagaag                                                                                                                        >  4:410233/1‑149 (MQ=255)
                          tATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAaagaag                                                                                                                        >  3:102414/1‑149 (MQ=255)
                          tATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCGCTTCCC‑GCAGGGAAGAGGGTTAAATaa                                                                                                                                 <  3:639428/140‑1 (MQ=255)
                           aTGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAaagaag                                                                                                                        >  2:212435/1‑148 (MQ=255)
                           aTGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAaagaag                                                                                                                        >  1:619929/1‑148 (MQ=255)
                             gAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAgatg                                                                                                                     >  4:349511/1‑149 (MQ=255)
                              aaGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAgatga                                                                                                                    <  4:265127/149‑1 (MQ=255)
                                gACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAgatgat                                                                                                                   <  2:675008/148‑1 (MQ=255)
                                    gTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGaa                                                                                                               <  1:69197/148‑1 (MQ=255)
                                     tCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAAc                                                                                                              >  1:810193/1‑148 (MQ=255)
                                           tAACTTCGATTGCTATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAAttttt                                                                                                       <  4:696072/149‑1 (MQ=255)
                                            aaCTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAAttttt                                                                                                       <  2:826524/148‑1 (MQ=255)
                                             aCTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCa                                                                                                     <  3:410233/149‑1 (MQ=255)
                                             aCTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCa                                                                                                     >  3:907581/1‑149 (MQ=255)
                                               ttCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAgg                                                                                                   >  3:929438/1‑149 (MQ=255)
                                               ttCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAg                                                                                                    >  2:238030/1‑148 (MQ=255)
                                                    ttGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCg                                                                                               >  2:522873/1‑148 (MQ=255)
                                                     tGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGa                                                                                              <  1:841412/148‑1 (MQ=255)
                                                        gATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGActct                                                                                          >  3:773286/1‑149 (MQ=255)
                                                           tCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGt                                                                                        >  2:819435/1‑148 (MQ=255)
                                                             cATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGaa                                                                                                               <  3:715111/123‑1 (MQ=255)
                                                             cATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGaa                                                                                                               >  4:715111/1‑123 (MQ=255)
                                                              aTTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGaa                                                                                    >  4:257210/1‑149 (MQ=255)
                                                              aTTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGaa                                                                                    >  4:557799/1‑149 (MQ=255)
                                                              aTTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGaa                                                                                    >  3:951680/1‑149 (MQ=255)
                                                              aTTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGa                                                                                     >  1:312309/1‑148 (MQ=255)
                                                               tttATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAAc                                                                                   <  3:726474/149‑1 (MQ=255)
                                                                        aaGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTgg                                                                          >  3:278077/1‑149 (MQ=255)
                                                                        aaGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTg                                                                           <  1:862268/148‑1 (MQ=255)
                                                                         aGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTgg                                                                          >  2:876572/1‑148 (MQ=255)
                                                                          gCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGaa                                                                        >  4:680030/1‑149 (MQ=255)
                                                                             aGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAAc                                                                                   <  4:338350/135‑1 (MQ=255)
                                                                             aGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAAc                                                                                   >  3:338350/1‑135 (MQ=255)
                                                                               aaaaGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGaa                                                                   <  3:349511/149‑1 (MQ=255)
                                                                                     cGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGcgcc                                                              <  2:567514/148‑1 (MQ=255)
                                                                                           cAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCa                                                       <  3:745589/149‑1 (MQ=255)
                                                                                            aGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAg                                                      <  3:157767/149‑1 (MQ=255)
                                                                                            aGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAg                                                      <  3:902157/149‑1 (MQ=255)
                                                                                                  tgGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGa                                                 <  2:29978/148‑1 (MQ=255)
                                                                                                  tgGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGa                                                 <  2:738534/148‑1 (MQ=255)
                                                                                                   ggCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAg                                               >  4:601382/1‑149 (MQ=255)
                                                                                                     cGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGa                                                                                                                              <  2:854636/68‑1 (MQ=255)
                                                                                                     cGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGa                                                                                                                              >  1:854636/1‑68 (MQ=255)
                                                                                                       tCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCg                                            >  1:170679/1‑148 (MQ=255)
                                                                                                        cGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGt                                           <  1:674956/148‑1 (MQ=255)
                                                                                                         ggatggatGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTc                                          <  2:130676/148‑1 (MQ=255)
                                                                                                         ggatggatGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCt                                         >  4:713555/1‑149 (MQ=255)
                                                                                                         ggatggatGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGACGAAGCCGTc                                          <  1:87984/148‑1 (MQ=255)
                                                                                                          gatggatgCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCt                                         >  2:642908/1‑148 (MQ=255)
                                                                                                          gatggatgCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTg                                        <  4:48518/149‑1 (MQ=255)
                                                                                                             ggatgCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGtt                                     <  4:773286/149‑1 (MQ=255)
                                                                                                                   gcTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCgc                               <  4:655181/149‑2 (MQ=255)
                                                                                                                     tGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCgggg                             <  3:622475/149‑1 (MQ=255)
                                                                                                                        cATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGgc                           >  1:438528/1‑148 (MQ=255)
                                                                                                                         aTAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGTAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGgcag                         >  3:136526/1‑149 (MQ=255)
                                                                                                                          tAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGGGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGgcag                         >  2:866212/1‑148 (MQ=255)
                                                                                                                            aaTGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCaa                      <  3:325295/149‑1 (MQ=255)
                                                                                                                              tGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAAc                                                                                   <  3:25652/86‑1 (MQ=255)
                                                                                                                               gACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACt                   >  3:575028/1‑149 (MQ=255)
                                                                                                                                  aGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCaa                <  3:463627/149‑1 (MQ=255)
                                                                                                                                   gATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCaa                >  2:123647/1‑148 (MQ=255)
                                                                                                                                   gATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGCCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCaa                <  1:212435/148‑1 (MQ=255)
                                                                                                                                    aTGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAAt               <  1:908259/148‑1 (MQ=255)
                                                                                                                                     tGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATg              <  2:730812/148‑1 (MQ=255)
                                                                                                                                      gCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGc             >  1:174207/1‑148 (MQ=255)
                                                                                                                                        ccTCTTCCCTGC‑GGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAAttttt                                                                                                       >  4:892203/1‑56 (MQ=25)
                                                                                                                                        ccTCTTCCCTGC‑GGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTac          >  3:443710/1‑149 (MQ=255)
                                                                                                                                        ccTCTTCCC‑TGCGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAAttttt                                                                                                       <  3:892203/56‑1 (MQ=37)
                                                                                                                                        ccTCTTCCC‑TGCGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCt                                                                            <  2:363153/83‑1 (MQ=39)
                                                                                                                                        ccTCTTCCC‑TGCGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCt                                                                            >  1:363153/1‑83 (MQ=39)
                                                                                                                                        ccTCTTCCC‑TGCGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTac          >  4:755121/1‑149 (MQ=255)
                                                                                                                                         ctctTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTaca         >  4:521277/1‑149 (MQ=255)
                                                                                                                                         ctctTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTaca         >  4:255292/1‑149 (MQ=255)
                                                                                                                                         ctctTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTac          <  2:427184/148‑1 (MQ=255)
                                                                                                                                                c‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGtt                                     >  4:396170/1‑114 (MQ=255)
                                                                                                                                                c‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGtt                                     <  3:396170/114‑1 (MQ=255)
                                                                                                                                                c‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAAccc  >  4:654957/1‑149 (MQ=255)
                                                                                                                                                  gCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAAccc  >  2:119813/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAACCC  >  NC_000913/3020365‑3020657

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑