Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,548,401 T→C 5.4% C39C (TGT→TGC) ‡ 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,4010TC 5.4% 441.5 / 10.1 147C39C (TGT→TGC) ‡fimGminor component of type 1 fimbriae
Reads supporting (aligned to +/- strand):  ref base T (66/73);  new base C (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 = 1.00e+00

Only 100 of 154 total aligned reads displayed.

TTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTGGGAACGTCGAGGGTCACTGCCAGCTTCAG  >  NC_000913/4548252‑4548550
                                                                                                                                                      |                                                                                                                                                     
ttCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCgtgt                                                                                                                                                       <  2:470597/150‑1 (MQ=255)
   aCTCTTGAATATCAGTCACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACg                                                                                                                                                    <  1:196666/150‑1 (MQ=255)
       ttGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGttt                                                                                                                                                >  1:1152654/1‑150 (MQ=255)
        tGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTc                                                                                                                                               <  2:598377/150‑1 (MQ=255)
         gAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTACACGGtttgg                                                                                                                                              <  2:415883/150‑3 (MQ=255)
         gAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGAGCTCGCAAGTGCGACGATACAGGCATCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTcc                                                                                                                                              >  1:614013/1‑150 (MQ=255)
          aaTATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTcca                                                                                                                                             <  2:922051/150‑1 (MQ=255)
            tatCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTccacc                                                                                                                                           <  2:1329407/150‑1 (MQ=255)
            tatCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTccacc                                                                                                                                           <  1:1122162/150‑1 (MQ=255)
            tatCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTACACGGtttggcga                                                                                                                                           <  1:780787/150‑6 (MQ=255)
             atCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTccacaa                                                                                                                                          >  2:1165205/1‑148 (MQ=255)
                 gTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCaa                                                                                                                                      >  2:513084/1‑150 (MQ=255)
                   aaCTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATg                                                                                                                                    >  2:1372764/1‑150 (MQ=255)
                   aaCTGCAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGAGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTATCCACCACCAATg                                                                                                                                    >  1:1030239/1‑150 (MQ=255)
                    aCTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGc                                                                                                                                   >  2:994099/1‑150 (MQ=255)
                     cTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATCTTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGcc                                                                                                                                  <  1:686139/150‑1 (MQ=255)
                      tGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCCCCACCAATGCCa                                                                                                                                 >  2:1003561/1‑150 (MQ=255)
                             gCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGa                                                                                                                          <  2:481353/150‑1 (MQ=255)
                              cTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGAt                                                                                                                         >  2:357025/1‑150 (MQ=255)
                                cATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGAtc                                                                                                                        >  2:1705918/1‑150 (MQ=255)
                                  tGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCg                                                                                                                     >  1:1611372/1‑150 (MQ=255)
                                  tGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACAACCAATGCCACGGTTGATCTCg                                                                                                                     >  1:1571569/1‑150 (MQ=255)
                                   gAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCgg                                                                                                                    <  2:727787/150‑1 (MQ=255)
                                    aaaTGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGc                                                                                                                   <  2:886692/150‑1 (MQ=255)
                                     aaTGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCg                                                                                                                  <  2:945321/150‑1 (MQ=255)
                                        ggTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATc                                                                                                               <  1:1180752/150‑1 (MQ=255)
                                          tGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCtt                                                                                                             >  2:1044261/1‑150 (MQ=255)
                                             aaaCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTAt                                                                                                          <  2:1538057/150‑1 (MQ=255)
                                               aCGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTc                                                                                                        <  1:849597/150‑1 (MQ=255)
                                                  tGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCCCCAATGCCACGGTTGATCTCGGCGATCTTTATTCttt                                                                                                     >  1:308841/1‑150 (MQ=255)
                                                   gggTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTc                                                                                                    <  1:887286/150‑1 (MQ=255)
                                                    ggTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCCCGGTTGATCTCGGCGATCTTTATTCTTTcc                                                                                                   >  2:980495/1‑149 (MQ=255)
                                                    ggTATGTAATGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGAACTCGGCGATCTTTATTCTTTCa                                                                                                   >  1:855323/1‑150 (MQ=255)
                                                     gtatgtatTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAg                                                                                                  <  1:1104661/150‑1 (MQ=255)
                                                      tatgtatTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGt                                                                                                 >  1:1467857/1‑150 (MQ=255)
                                                      tatgtatTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGt                                                                                                 >  1:1284762/1‑150 (MQ=255)
                                                        tgtatTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCt                                                                                               <  1:1612225/150‑1 (MQ=255)
                                                           atTGGCGGCATTATTGGCGCTCGTAAGTGCTACGATACAGGCAGCCGATGTCACCATCACGGTGATCGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTAt                                                                                            <  2:241144/150‑1 (MQ=255)
                                                            ttGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATg                                                                                           <  2:1634754/150‑1 (MQ=255)
                                                            ttGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATg                                                                                           <  1:994099/150‑1 (MQ=255)
                                                              ggcggcAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTc                                                                                         >  2:1389302/1‑150 (MQ=255)
                                                                  gcAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGcc                                                                                     <  2:1078519/150‑1 (MQ=255)
                                                                  gcAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGcc                                                                                     >  1:731802/1‑150 (MQ=255)
                                                                  gcAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGcc                                                                                     >  1:1078519/1‑150 (MQ=255)
                                                                   cAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCg                                                                                    <  2:1285963/150‑1 (MQ=255)
                                                                      tatTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCgggg                                                                                 <  1:936503/150‑1 (MQ=255)
                                                                             gcTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATc                                                                          <  1:1051179/150‑1 (MQ=255)
                                                                              cTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCg                                                                         >  1:359616/1‑150 (MQ=255)
                                                                              cTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCg                                                                         >  2:54660/1‑150 (MQ=255)
                                                                              cTCGCAAGTGCGACGAAACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCg                                                                         <  2:1074762/150‑1 (MQ=255)
                                                                               tCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCgg                                                                        >  2:42068/1‑150 (MQ=255)
                                                                               tCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTATTTCAGTCTTATGTCTGACGGGGCGGCATCgg                                                                        >  2:575555/1‑150 (MQ=255)
                                                                                cGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGc                                                                       <  1:398551/150‑1 (MQ=255)
                                                                                 gCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGcc                                                                      <  1:833451/150‑1 (MQ=255)
                                                                                  cAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGTATCGGCCt                                                                     <  2:582816/150‑1 (MQ=255)
                                                                                   aaGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTg                                                                    <  1:683307/150‑1 (MQ=255)
                                                                                    aGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTgg                                                                   <  1:165895/150‑1 (MQ=255)
                                                                                    aGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTgg                                                                   <  1:1372764/150‑1 (MQ=255)
                                                                                    aGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTgg                                                                   <  1:808632/150‑1 (MQ=255)
                                                                                     gTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGc                                                                  <  2:522746/150‑1 (MQ=255)
                                                                                      tGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCGACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCa                                                                 <  2:1021265/150‑1 (MQ=255)
                                                                                        cgacgaTACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCatg                                                               >  2:808474/1‑150 (MQ=255)
                                                                                            gaTACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGt                                                           >  2:610169/1‑150 (MQ=255)
                                                                                             aTACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCCCCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCCGTCTTATGTCTGCCGGGGCGGCAGCGGCCTGGCATGATGtt                                                          >  2:828565/1‑150 (MQ=255)
                                                                                              tACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTg                                                         <  1:513084/150‑1 (MQ=255)
                                                                                              tACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTg                                                         <  1:913385/150‑1 (MQ=255)
                                                                                               aCAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTgc                                                        >  2:1167801/1‑150 (MQ=255)
                                                                                                cAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTgcg                                                       <  1:1044261/150‑1 (MQ=255)
                                                                                                  ggCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCt                                                     <  1:27595/150‑1 (MQ=255)
                                                                                                      gCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAg                                                 <  2:1599623/150‑1 (MQ=255)
                                                                                                         gATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCCGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTg                                              >  2:1406753/1‑150 (MQ=255)
                                                                                                              cACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACt                                           >  2:1354808/1‑148 (MQ=255)
                                                                                                              cACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTCGAGTTGACt                                           <  1:1354808/148‑1 (MQ=255)
                                                                                                               aCCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAAt                                        >  2:517097/1‑150 (MQ=255)
                                                                                                                ccATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAAtt                                       <  1:42068/150‑1 (MQ=255)
                                                                                                                 cATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTg                                      <  2:359616/150‑1 (MQ=255)
                                                                                                                  aTCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACTGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCTTCGGACTGGCATGATGTTGCGCTTGAGTTTACTAATTGt                                     >  2:1587773/1‑150 (MQ=255)
                                                                                                                    cACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTcc                                   <  1:575555/150‑1 (MQ=255)
                                                                                                                     aCGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCg                                  <  2:1467857/150‑1 (MQ=255)
                                                                                                                      cGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTgcg                                                       >  1:1141439/1‑128 (MQ=255)
                                                                                                                      cGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTgcg                                                       <  2:1141439/128‑1 (MQ=255)
                                                                                                                          gAACGGTAAGGTCATCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTgcgc                                                      <  2:1653963/125‑1 (MQ=255)
                                                                                                                          gAACGGTAAGGTCATCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTgcgc                                                      >  1:1653963/1‑125 (MQ=255)
                                                                                                                            aCGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACAAATTGTCCGGTGGGaa                           <  2:972819/150‑1 (MQ=255)
                                                                                                                              ggTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTGGGAACg                         >  2:652066/1‑150 (MQ=255)
                                                                                                                                   ggTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTATTTGTCCGGTGGGAACGTCGAg                    >  1:836291/1‑150 (MQ=255)
                                                                                                                                         ggtggaaaCCGTACACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTGGGAACGTCGAGGGTCa               >  2:1253175/6‑150 (MQ=255)
                                                                                                                                         cgCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGGCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTGGGAACGTCGAGGGTCAc              <  2:90135/150‑1 (MQ=255)
                                                                                                                                          gtggaaaCCGTACACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTGGGAACGTCGAGGGTCAc              >  2:374170/5‑150 (MQ=255)
                                                                                                                                           tggaaaCCGTACACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTgg                              <  1:722674/130‑1 (MQ=255)
                                                                                                                                           tggaaaCCGTACACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTgg                              >  2:722674/4‑133 (MQ=255)
                                                                                                                                           cc‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGGTTGAGTTGACTAATTGTCCGGTGGGAACGTCGAGGGTCACTg            >  1:303415/1‑150 (MQ=255)
                                                                                                                                           cc‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTGGGAACGTCGAGGGTCACTg            >  1:579430/1‑150 (MQ=255)
                                                                                                                                            ggaaaCCGTACACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTGGGAACGTCGAGGGTCACTg            >  1:131315/3‑150 (MQ=255)
                                                                                                                                            caaaaCCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCgg                                 <  1:138809/127‑1 (MQ=255)
                                                                                                                                               aaCCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGAGCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTGGGAACGTCGAGGGTCACTGCCa         <  1:50173/150‑1 (MQ=255)
                                                                                                                                                    tgtACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTGGGAACGTCGAGGGTCACTGCCAGCTTc    <  1:652066/150‑1 (MQ=255)
                                                                                                                                                    tgtACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTGGGAACGTCGAGGGTCACTGCCAGCTTc    <  1:808474/150‑1 (MQ=255)
                                                                                                                                                    tgtACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTGGGAACGTCGAGGGTCACTGCCAGCTTc    >  2:1112578/1‑150 (MQ=255)
                                                                                                                                                     gtACGGTTTCAACCACCAATGCCACGGTTGATCTCGGCGATATTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTGGGAACGTCGAGGGTCACTGCCAGCTTCa   >  1:1295253/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TTCACTCTTGAATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCC‑AAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTGGGAACGTCGAGGGTCACTGCCAGCTTCAG  >  NC_000913/4548252‑4548550

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: ‑