Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,548,400 G→A 5.4% C39Y (TGT→TAT) ‡ fimG → minor component of type 1 fimbriae

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,548,4000GA 5.4% 476.8 / 9.4 147C39Y (TGT→TAT) ‡fimGminor component of type 1 fimbriae
Reads supporting (aligned to +/- strand):  ref base G (66/73);  new base A (5/3);  total (71/76)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01

Only 100 of 155 total aligned reads displayed.

CTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTGGGAACGTCGAGGGTCACTGCCAGCTTCA  >  NC_000913/4548251‑4548549
                                                                                                                                                     |                                                                                                                                                     
cTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCgtg                                                                                                                                                       <  2:1361542/150‑1 (MQ=255)
cTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCgtg                                                                                                                                                       <  2:58445/150‑1 (MQ=255)
 ttCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCgtgt                                                                                                                                                      <  2:470597/150‑1 (MQ=255)
    aCTCTTGAATATCAGTCACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACg                                                                                                                                                   <  1:196666/150‑1 (MQ=255)
       ctTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGtt                                                                                                                                                >  2:653995/1‑150 (MQ=255)
        ttGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGttt                                                                                                                                               >  1:1152654/1‑150 (MQ=255)
        ttGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGttt                                                                                                                                               >  2:517904/1‑150 (MQ=255)
         tGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTc                                                                                                                                              <  2:598377/150‑1 (MQ=255)
          gAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTACACGGtttgg                                                                                                                                             <  2:415883/150‑3 (MQ=255)
          gAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGAGCTCGCAAGTGCGACGATACAGGCATCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTcc                                                                                                                                             >  1:614013/1‑150 (MQ=255)
           aaTATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTcca                                                                                                                                            <  2:922051/150‑1 (MQ=255)
             tatCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTccacc                                                                                                                                          <  1:1589887/150‑1 (MQ=255)
             tatCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTccacc                                                                                                                                          <  2:1329407/150‑1 (MQ=255)
             tatCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTccacc                                                                                                                                          <  1:1122162/150‑1 (MQ=255)
             tatCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTACACGGtttggcga                                                                                                                                          <  1:780787/150‑6 (MQ=255)
              atCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTccacaa                                                                                                                                         >  2:1165205/1‑148 (MQ=255)
                    aaCTGGAGATGCTCATGAAATGGGGCAAACGGGGGTATGTATTGGCCGCAATATTGGCGCTCGGAAGTGCGACGATACAGGCAGCCGATGGCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATg                                                                                                                                   <  2:426605/150‑1 (MQ=255)
                    aaCTGCAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGAGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTATCCACCACCAATg                                                                                                                                   >  1:1030239/1‑150 (MQ=255)
                     aCTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGc                                                                                                                                  >  2:994099/1‑150 (MQ=255)
                      cTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATCTTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGcc                                                                                                                                 <  1:686139/150‑1 (MQ=255)
                       tGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCCCCACCAATGCCa                                                                                                                                >  2:1003561/1‑150 (MQ=255)
                        ggAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCAc                                                                                                                               <  2:1688319/150‑1 (MQ=255)
                             tGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTg                                                                                                                          >  1:557832/1‑150 (MQ=255)
                                tCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACTGTTTCCACCACCAATGCCACGGTTGAta                                                                                                                       >  2:1452792/1‑149 (MQ=255)
                                 cATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGAtct                                                                                                                      >  1:1285963/1‑150 (MQ=255)
                                   tGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACAACCAATGCCACGGTTGATCTCg                                                                                                                    >  1:1571569/1‑150 (MQ=255)
                                    gAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCgg                                                                                                                   <  2:727787/150‑1 (MQ=255)
                                     aaaTGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGc                                                                                                                  <  2:886692/150‑1 (MQ=255)
                                      aaTGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCg                                                                                                                 <  2:945321/150‑1 (MQ=255)
                                         ggTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATc                                                                                                              >  2:398551/1‑150 (MQ=255)
                                          gTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCCCCAATGCCACGGTTGATCTCGGCGATCt                                                                                                             >  2:1602379/1‑150 (MQ=255)
                                           tGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCtt                                                                                                            >  2:1044261/1‑150 (MQ=255)
                                                aCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTc                                                                                                       <  1:849597/150‑1 (MQ=255)
                                                   tGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCCCCAATGCCACGGTTGATCTCGGCGATCTTTATTCttt                                                                                                    >  1:308841/1‑150 (MQ=255)
                                                    gggTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTc                                                                                                   <  1:887286/150‑1 (MQ=255)
                                                     ggTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCCCGGTTGATCTCGGCGATCTTTATTCTTTcc                                                                                                  >  2:980495/1‑149 (MQ=255)
                                                     ggTATGTAATGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGAACTCGGCGATCTTTATTCTTTCa                                                                                                  >  1:855323/1‑150 (MQ=255)
                                                      gtatgtatTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAg                                                                                                 <  1:1104661/150‑1 (MQ=255)
                                                       tatgtatTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGt                                                                                                >  1:1284762/1‑150 (MQ=255)
                                                         tgtatTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCt                                                                                              <  1:1612225/150‑1 (MQ=255)
                                                            atTGGCGGCATTATTGGCGCTCGTAAGTGCTACGATACAGGCAGCCGATGTCACCATCACGGTGATCGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTAt                                                                                           <  2:241144/150‑1 (MQ=255)
                                                             ttGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATg                                                                                          <  2:1634754/150‑1 (MQ=255)
                                                             ttGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATg                                                                                          <  1:994099/150‑1 (MQ=255)
                                                             ttGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGACGATGTCACCCTCACGGTGAACGGTAAGGTCGGCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATg                                                                                          <  2:1550664/150‑1 (MQ=255)
                                                              tGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGt                                                                                         >  2:50173/1‑150 (MQ=255)
                                                               ggcggcAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTc                                                                                        >  2:1389302/1‑150 (MQ=255)
                                                                   gcAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGcc                                                                                    >  1:1078519/1‑150 (MQ=255)
                                                                   gcAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGcc                                                                                    >  1:731802/1‑150 (MQ=255)
                                                                   gcAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGcc                                                                                    <  2:1078519/150‑1 (MQ=255)
                                                                    cAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCg                                                                                   <  2:1285963/150‑1 (MQ=255)
                                                                       tatTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCgggg                                                                                <  1:936503/150‑1 (MQ=255)
                                                                          tGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGggcgg                                                                             <  1:43103/150‑1 (MQ=255)
                                                                           ggCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGggcggc                                                                            >  1:788064/1‑150 (MQ=255)
                                                                              gcTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATc                                                                         <  1:1051179/150‑1 (MQ=255)
                                                                               cTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCg                                                                        >  2:54660/1‑150 (MQ=255)
                                                                               cTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCg                                                                        >  1:359616/1‑150 (MQ=255)
                                                                               cTCGCAAGTGCGACGAAACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCg                                                                        <  2:1074762/150‑1 (MQ=255)
                                                                                 cGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGc                                                                      >  1:163171/1‑150 (MQ=255)
                                                                                  gCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGcc                                                                     <  1:833451/150‑1 (MQ=255)
                                                                                   cAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGTATCGGCCt                                                                    <  2:582816/150‑1 (MQ=255)
                                                                                    aaGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTg                                                                   <  1:683307/150‑1 (MQ=255)
                                                                                     aGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTgg                                                                  <  1:808632/150‑1 (MQ=255)
                                                                                     aGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTgg                                                                  <  1:165895/150‑1 (MQ=255)
                                                                                         cgacgaTACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCatg                                                              >  2:808474/1‑150 (MQ=255)
                                                                                             gaTACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGt                                                          >  2:610169/1‑150 (MQ=255)
                                                                                              aTACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCCCCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCCGTCTTATGTCTGCCGGGGCGGCAGCGGCCTGGCATGATGtt                                                         >  2:828565/1‑150 (MQ=255)
                                                                                               tACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTg                                                        <  1:513084/150‑1 (MQ=255)
                                                                                               tACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTg                                                        <  1:913385/150‑1 (MQ=255)
                                                                                                aCAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTgc                                                       >  2:1167801/1‑150 (MQ=255)
                                                                                                 cAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTgcg                                                      <  1:1044261/150‑1 (MQ=255)
                                                                                                  aGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGAGCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTgcgc                                                     <  2:555035/150‑1 (MQ=255)
                                                                                                   ggCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCt                                                    <  1:27595/150‑1 (MQ=255)
                                                                                                         cGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGtt                                              >  2:1013041/1‑150 (MQ=255)
                                                                                                            tGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGAc                                           <  1:1464750/150‑1 (MQ=255)
                                                                                                               cACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTCGAGTTGACt                                          <  1:1354808/148‑1 (MQ=255)
                                                                                                                aCCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAAt                                       >  2:517097/1‑150 (MQ=255)
                                                                                                                 ccATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAAtt                                      <  1:42068/150‑1 (MQ=255)
                                                                                                                  cATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTg                                     <  2:359616/150‑1 (MQ=255)
                                                                                                                   aTCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACTGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCTTCGGACTGGCATGATGTTGCGCTTGAGTTTACTAATTGt                                    >  2:1587773/1‑150 (MQ=255)
                                                                                                                    tCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTc                                   <  1:1366825/150‑1 (MQ=255)
                                                                                                                      aCGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCg                                 <  2:1467857/150‑1 (MQ=255)
                                                                                                                       cGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTgcg                                                      >  1:1141439/1‑128 (MQ=255)
                                                                                                                       cGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTgcg                                                      <  2:1141439/128‑1 (MQ=255)
                                                                                                                             aCGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACAAATTGTCCGGTGGGaa                          <  2:972819/150‑1 (MQ=255)
                                                                                                                               ggTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTGGGAACg                        >  2:652066/1‑150 (MQ=255)
                                                                                                                                    ggTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTATTTGTCCGGTGGGAACGTCGAg                   >  1:836291/1‑150 (MQ=255)
                                                                                                                                      tcgtcgCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAg                                                >  1:369461/1‑119 (MQ=255)
                                                                                                                                      tcgtcgCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAg                                                <  2:369461/119‑1 (MQ=255)
                                                                                                                                      gcgtcgCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTGGGAACGTCGAggg                 <  2:1369306/149‑1 (MQ=255)
                                                                                                                                         ggtggaaaCCGTACACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTGGGAACGTCGAGGGTCa              >  2:1253175/6‑150 (MQ=255)
                                                                                                                                          cgCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGGCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTGGGAACGTCGAGGGTCAc             <  2:90135/150‑1 (MQ=255)
                                                                                                                                           tggaaaCCGTACACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTgg                             <  1:722674/130‑1 (MQ=255)
                                                                                                                                            ggaaaCCGTACACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTGGGAACGTCGAGGGTCACTg           >  1:131315/3‑150 (MQ=255)
                                                                                                                                            ggaaaCCGTACACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTGGGAACGTAGAGGGTCACTg           >  2:347993/3‑150 (MQ=255)
                                                                                                                                            ccAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGGTTGAGTTGACTAATTGTCCGGTGGGAACGTCGAGGGTCACTg           >  1:303415/1‑150 (MQ=255)
                                                                                                                                            ccAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTGGGAACGTCGAGGGTCACTg           >  1:579430/1‑150 (MQ=255)
                                                                                                                                            caaaaCCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCgg                                >  2:138809/3‑129 (MQ=255)
                                                                                                                                            caaaaCCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCgg                                <  1:138809/127‑1 (MQ=255)
                                                                                                                                                    tgtACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTGGGAACGTCGAGGGTCACTGCCAGCTTc   <  1:652066/150‑1 (MQ=255)
                                                                                                                                                    tgtACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTGGGAACGTCGAGGGTCACTGCCAGCTTc   >  2:1112578/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTGGGAACGTCGAGGGTCACTGCCAGCTTCA  >  NC_000913/4548251‑4548549

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