Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 309,680 G→T 21.8% A116D (GCC→GAC)  ecpB ← ECP production pilus chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913309,6800GT21.8% 329.0 / 76.4 124A116D (GCC→GAC) ecpBECP production pilus chaperone
Reads supporting (aligned to +/- strand):  ref base G (40/57);  new base T (10/17);  total (50/74)
Fisher's exact test for biased strand distribution p-value = 8.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

Only 100 of 125 total aligned reads displayed.

GCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAAC  >  NC_000913/309533‑309828
                                                                                                                                                   |                                                                                                                                                    
gCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGc                                                                                                                                                     >  2:197782/1‑149 (MQ=255)
   tttCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTg                                                                                                                                                  >  2:207059/1‑149 (MQ=255)
    ttCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGTCCTGa                                                                                                                                                 <  2:672951/149‑1 (MQ=255)
     tCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGTCCTGAt                                                                                                                                                <  2:181701/149‑1 (MQ=255)
     tCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGAt                                                                                                                                                <  1:663864/149‑1 (MQ=255)
     tCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGAt                                                                                                                                                <  2:660733/149‑1 (MQ=255)
       cTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCa                                                                                                                                              <  3:431725/149‑1 (MQ=255)
        tGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGGGGTCCTGATCaa                                                                                                                                             <  1:532348/149‑1 (MQ=255)
          tATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGTCCTGATCAAAc                                                                                                                                           >  4:746018/1‑149 (MQ=255)
            tttGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCa                                                                                                                                         <  1:394855/149‑1 (MQ=255)
                tCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACa                                                                                                                                     <  3:793661/149‑1 (MQ=255)
                 cAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACaa                                                                                                                                    <  4:374872/149‑1 (MQ=255)
                  aGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAGTGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAAt                                                                                                                                   <  3:602769/149‑1 (MQ=255)
                    ggAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGc                                                                                                                                 >  2:845411/1‑149 (MQ=255)
                     gAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGTCATGATCAAACCAGACAATGCg                                                                                                                                <  2:847213/149‑1 (MQ=255)
                        ccGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGgta                                                                                                                             <  1:645837/149‑1 (MQ=255)
                         cGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGgtag                                                                                                                            <  4:128910/149‑1 (MQ=255)
                         cGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGgtag                                                                                                                            >  1:542912/1‑149 (MQ=255)
                          gTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGgtagt                                                                                                                           >  1:26616/1‑149 (MQ=255)
                          gTTGGCGGACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGTCCTGATCAAACCAGACAATGCGgtagt                                                                                                                           >  2:504133/1‑149 (MQ=255)
                           ttGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGTCCTGATCAAACCAGACAATGCGgtagta                                                                                                                          <  3:267607/149‑1 (MQ=255)
                                    tGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCttt                                                                                                                 <  3:286996/149‑1 (MQ=255)
                                    tGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCttt                                                                                                                 >  4:905976/1‑149 (MQ=255)
                                    tGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCttt                                                                                                                 >  1:357694/1‑149 (MQ=255)
                                     gAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGTCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCtttt                                                                                                                <  4:20963/149‑1 (MQ=255)
                                     gAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGTCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCtttt                                                                                                                <  4:505758/149‑1 (MQ=255)
                                      aaaGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCttttt                                                                                                               <  4:417297/149‑1 (MQ=255)
                                          tGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATc                                                                                                           >  1:430751/1‑149 (MQ=255)
                                            gTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGTCCTGATCAAACCAGACAAt                                                                                                                                   <  4:872345/123‑1 (MQ=255)
                                            gTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGTCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTg                                                                                                         >  3:926043/1‑149 (MQ=255)
                                               gTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCgg                                                                                                      <  4:344951/149‑1 (MQ=255)
                                                   gCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGTCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCtttt                                                                                                                >  3:170372/1‑135 (MQ=255)
                                                   gCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGTCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCtttt                                                                                                                <  4:170372/135‑1 (MQ=255)
                                                     cTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCt                                                                                                <  1:845411/149‑1 (MQ=255)
                                                         cGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGTGCATCACTGAGGTCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTtata                                                                                            <  2:133200/149‑1 (MQ=255)
                                                         cGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGTCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTtata                                                                                            <  1:785815/149‑1 (MQ=255)
                                                         cGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTtata                                                                                            <  1:207059/149‑1 (MQ=255)
                                                         cGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTtata                                                                                            <  1:324079/149‑1 (MQ=255)
                                                             gggCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATagaag                                                                                        >  4:365984/1‑149 (MQ=255)
                                                               gCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCa                                                                                                                                         <  1:410844/98‑1 (MQ=255)
                                                               gCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCa                                                                                                                                         >  2:410844/1‑98 (MQ=255)
                                                                 gacgacCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGTCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGc                                                                                    <  4:816430/149‑1 (MQ=255)
                                                                    gacCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGa                                                                                 <  4:64659/149‑1 (MQ=255)
                                                                    gacCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGa                                                                                 <  3:732683/149‑1 (MQ=255)
                                                                    gacCAGAACGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGa                                                                                 >  3:853662/1‑149 (MQ=255)
                                                                     acCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGAt                                                                                >  2:897693/1‑149 (MQ=255)
                                                                        aGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCAc                                                                             <  3:199654/149‑1 (MQ=255)
                                                                         gAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGTCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACt                                                                            >  3:612468/1‑149 (MQ=255)
                                                                         gAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACt                                                                            >  3:648351/1‑149 (MQ=255)
                                                                          aaTGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACtt                                                                           <  2:131536/149‑1 (MQ=255)
                                                                           aTGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGTCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTc                                                                          <  1:538591/149‑1 (MQ=255)
                                                                           aTGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTc                                                                          <  4:897644/149‑1 (MQ=255)
                                                                              gTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGTCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCt                                                                       >  1:261237/1‑149 (MQ=255)
                                                                                 ccGATGCGGGAGGAAGAAGTGGCCACAGCGCTGCGGTTGGCAGTATCGCGCTGCACATCACTGAGGGCCTGATCAAACCTGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGTTCCCTTATAGAAGAAGCTGATCACTTCGCTGGc                                                                    <  4:343873/149‑1 (MQ=255)
                                                                                    aTGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTg                                                                 >  1:487067/1‑149 (MQ=255)
                                                                                     tGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTgg                                                                <  4:809531/149‑1 (MQ=255)
                                                                                        gggcggAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCg                                                             >  4:226054/1‑149 (MQ=255)
                                                                                         ggcggAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCgg                                                            <  2:695427/149‑1 (MQ=255)
                                                                                          gcggAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCt                                                                   <  3:457067/141‑1 (MQ=255)
                                                                                          gcggAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCggg                                                           <  1:131998/149‑1 (MQ=255)
                                                                                            ggAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTa                                                         <  2:326366/149‑1 (MQ=255)
                                                                                            ggAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTa                                                         >  2:639272/1‑149 (MQ=255)
                                                                                                   gTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCaa                                                  <  3:159294/149‑1 (MQ=255)
                                                                                                    tGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGc                                                                                                        <  2:31202/94‑1 (MQ=255)
                                                                                                    tGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGc                                                                                                        >  1:31202/1‑94 (MQ=255)
                                                                                                    tGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAg                                                 >  3:82915/1‑149 (MQ=255)
                                                                                                       ccACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTg                                              <  4:82915/149‑1 (MQ=255)
                                                                                                       ccACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTg                                              <  4:314531/149‑1 (MQ=255)
                                                                                                           aGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCgg                                          >  3:579370/1‑149 (MQ=255)
                                                                                                              gcTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCt                                                                   >  2:167035/1‑121 (MQ=255)
                                                                                                              gcTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCt                                                                   <  1:167035/121‑1 (MQ=255)
                                                                                                              gcTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGgag                                       <  3:905976/149‑1 (MQ=255)
                                                                                                                tGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGgagtg                                     <  2:525781/149‑1 (MQ=255)
                                                                                                                tGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGAGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGgagtg                                     >  1:153716/1‑149 (MQ=255)
                                                                                                                 gCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGgagtga                                    <  3:305541/149‑1 (MQ=255)
                                                                                                                  cGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGgagtgag                                   <  2:89807/149‑1 (MQ=255)
                                                                                                                    gTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTa                                 >  2:270007/1‑149 (MQ=255)
                                                                                                                     ttGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAg                                <  2:487067/149‑1 (MQ=255)
                                                                                                                        gCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCaa                             <  1:197782/149‑1 (MQ=255)
                                                                                                                        gCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCaa                             <  1:639272/149‑1 (MQ=255)
                                                                                                                         cATTATCGCGCTGTGCATCACTGAGGTCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAAc                            >  3:734221/1‑149 (MQ=255)
                                                                                                                         cATTATCGCGCTGCGCATCACTGAGGTCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAAc                            >  2:124549/1‑149 (MQ=255)
                                                                                                                          actATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCATTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTag                                                        >  4:664499/3‑120 (MQ=255)
                                                                                                                          aTTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGAt                                                                                >  3:192369/1‑96 (MQ=255)
                                                                                                                          aTTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTag                                                        <  3:664499/120‑1 (MQ=255)
                                                                                                                          aTTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACt                           >  1:297270/1‑149 (MQ=255)
                                                                                                                             aTCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCAt                        >  1:244672/1‑149 (MQ=255)
                                                                                                                              tCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATc                       <  2:542912/149‑1 (MQ=255)
                                                                                                                               cgcgCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATcc                      >  4:842932/1‑149 (MQ=255)
                                                                                                                                  gcTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGc                   <  3:365984/149‑1 (MQ=255)
                                                                                                                                        cATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCg                                           <  4:872182/119‑1 (MQ=255)
                                                                                                                                        cATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCg                                           >  3:872182/1‑119 (MQ=255)
                                                                                                                                         aTCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCa            >  1:546125/1‑149 (MQ=255)
                                                                                                                                         aTCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCa            >  1:668979/1‑149 (MQ=255)
                                                                                                                                         aTCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGACGGGAGTGAGTAGCAACTCATCCGGCTTGTCCa            <  3:585389/149‑1 (MQ=255)
                                                                                                                                          tCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCAt           >  1:71417/1‑149 (MQ=255)
                                                                                                                                           cACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCAtt          <  2:134043/149‑1 (MQ=255)
                                                                                                                                            aCTGAGGGCCTGATCAAACGAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTg         >  3:543582/1‑149 (MQ=255)
                                                                                                                                            aCTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTg         <  3:334522/149‑1 (MQ=255)
                                                                                                                                              tGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTgag       >  1:727380/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATGCGGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAAC  >  NC_000913/309533‑309828

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: