Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,977,962 G→T 80.8% T79K (ACG→AAG)  flhD ← flagellar class II regulon transcriptional activator, with FlhC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,977,9620GT80.8% 287.9 / 78.7 125T79K (ACG→AAG) flhDflagellar class II regulon transcriptional activator, with FlhC
Reads supporting (aligned to +/- strand):  ref base G (12/12);  new base T (48/53);  total (60/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

Only 100 of 125 total aligned reads displayed.

TTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCG  >  NC_000913/1977818‑1978105
                                                                                                                                                |                                                                                                                                               
ttCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCaa                                                                                                                                             <  1:71153/149‑1 (MQ=255)
    tGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGa                                                                                                                                         >  1:221090/1‑149 (MQ=255)
    tGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGa                                                                                                                                         >  2:340182/1‑149 (MQ=255)
    tGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGa                                                                                                                                         >  3:675613/1‑149 (MQ=255)
    tGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGa                                                                                                                                         >  1:122390/1‑149 (MQ=255)
      aaCAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGt                                                                                                                                       >  1:417570/1‑149 (MQ=255)
        cAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTaa                                                                                                                                     <  2:187836/149‑1 (MQ=255)
           tGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCg                                                                                                                                  >  4:483180/1‑149 (MQ=255)
           tGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCg                                                                                                                                  >  2:23685/1‑149 (MQ=255)
           tGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCg                                                                                                                                  >  4:526494/1‑149 (MQ=255)
            gCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGt                                                                                                                                 <  1:466622/149‑1 (MQ=255)
                  ttCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCtggtg                                                                                                                           <  3:326485/149‑1 (MQ=255)
                    cACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGc                                                                                                                         >  4:359323/1‑149 (MQ=255)
                     aCTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCt                                                                                                                        <  4:675613/149‑1 (MQ=255)
                       tCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATGATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGt                                                                                                                      <  2:770913/149‑1 (MQ=255)
                       tCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGt                                                                                                                      >  2:617909/1‑149 (MQ=255)
                       tCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGt                                                                                                                      >  4:741264/1‑149 (MQ=255)
                       tCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCCGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGt                                                                                                                      >  4:797191/1‑149 (MQ=255)
                           gATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCaaa                                                                                                                  >  2:495230/1‑149 (MQ=255)
                               aGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATc                                                                                                                                   <  3:775893/128‑1 (MQ=255)
                               aGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATc                                                                                                                                   >  4:775893/1‑128 (MQ=255)
                               aGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACgg                                                                                                              <  2:346455/149‑1 (MQ=255)
                               aGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGc                                                                                                                         <  4:582678/138‑1 (MQ=255)
                               aGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGc                                                                                                                         >  3:582678/1‑138 (MQ=255)
                                 gCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGaa                                                                                                            <  3:37690/149‑1 (MQ=255)
                                   ccTTTTCTTGTGCAGCGCTTCTTCAGTCTGATTATCTTCATTCAGCAAGCGTCTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGACTGTCAAAACGGagtt                                                                                                          <  4:895758/149‑4 (MQ=255)
                                   ccTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAGCATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGCAGt                                                                                                          <  2:141710/149‑1 (MQ=255)
                                   ccTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGt                                                                                                          >  2:516099/1‑149 (MQ=255)
                                     ttttCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGTGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGa                                                                                                        >  4:375399/1‑149 (MQ=255)
                                      tttCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGAc                                                                                                       <  1:617909/149‑1 (MQ=255)
                                              gcAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGt                                                                                               >  4:375483/1‑149 (MQ=255)
                                               cAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGtt                                                                                              <  2:569704/149‑1 (MQ=255)
                                                 gcgcTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGa                                                                                            <  3:67642/149‑1 (MQ=255)
                                                  cgcTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAAGCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGAt                                                                                           <  3:603539/149‑1 (MQ=255)
                                                   gcTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGAtt                                                                                          <  1:23685/149‑1 (MQ=255)
                                                   gcTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGAtt                                                                                          <  3:798034/149‑1 (MQ=255)
                                                   gcTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGAtt                                                                                          <  3:526494/149‑1 (MQ=255)
                                                          cAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCtg                                                                                                                              <  2:291679/106‑1 (MQ=255)
                                                          cAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCtg                                                                                                                              >  1:291679/1‑106 (MQ=255)
                                                           aGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTg                                                                                  <  3:640477/149‑1 (MQ=255)
                                                           aGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTg                                                                                  >  2:469516/1‑149 (MQ=255)
                                                           aGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTg                                                                                  >  4:375997/1‑149 (MQ=255)
                                                             gCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGcc                                                                                <  2:417570/149‑1 (MQ=255)
                                                             gCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGcc                                                                                >  1:702548/1‑149 (MQ=255)
                                                              cTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCa                                                                               <  4:891799/149‑1 (MQ=255)
                                                                 aTTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCt                                                                            >  4:736371/1‑149 (MQ=255)
                                                                  ttAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCtg                                                                                                                              <  1:93547/98‑1 (MQ=255)
                                                                  ttAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCtg                                                                                                                              >  2:93547/1‑98 (MQ=255)
                                                                    aaCATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTaa                                                                         <  3:229390/149‑1 (MQ=255)
                                                                     aCATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAAc                                                                        <  2:879474/149‑1 (MQ=255)
                                                                      catcatTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTTCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAAcc                                                                       >  3:600152/1‑149 (MQ=255)
                                                                      catcatTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAAcc                                                                       >  1:775512/1‑149 (MQ=255)
                                                                       atcatTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGtt                                                                                      >  4:887716/1‑133 (MQ=255)
                                                                       atcatTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGtt                                                                                      <  3:887716/133‑1 (MQ=255)
                                                                       atcatTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCa                                                                      >  2:633053/1‑149 (MQ=255)
                                                                          tttCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCAtt                                                                    >  3:902871/2‑148 (MQ=255)
                                                                          atTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCAttt                                                                   <  3:359323/149‑1 (MQ=255)
                                                                          atTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCAtt                                                                    <  4:902871/148‑1 (MQ=255)
                                                                           ttCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGt                                                                                               <  1:322970/120‑1 (MQ=255)
                                                                           ttCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTg                                                                  >  4:623499/1‑149 (MQ=255)
                                                                           ttCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTg                                                                  <  3:623499/149‑1 (MQ=255)
                                                                            tCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGc                                                                 <  4:123248/149‑1 (MQ=255)
                                                                            tCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGc                                                                 <  1:203302/149‑1 (MQ=255)
                                                                             cAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCg                                                                >  3:422297/1‑149 (MQ=255)
                                                                                cAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGaa                                                             <  3:375997/149‑1 (MQ=255)
                                                                                   gCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAagag                                                          >  2:794830/1‑149 (MQ=255)
                                                                                        ttGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGt                                                     <  2:122390/149‑1 (MQ=255)
                                                                                          gagaGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGc                                                   <  4:196412/149‑1 (MQ=255)
                                                                                              gCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCt                                               <  2:244342/149‑1 (MQ=255)
                                                                                                  gatgCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACg                                           >  1:782152/1‑149 (MQ=255)
                                                                                                      ccGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGt                                       <  2:839660/149‑1 (MQ=255)
                                                                                                       cGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTc                                      <  1:495230/149‑1 (MQ=255)
                                                                                                       cGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTc                                      <  2:221090/149‑1 (MQ=255)
                                                                                                         gTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGc                                    >  1:186176/1‑149 (MQ=255)
                                                                                                          tATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGcc                                   >  1:435982/1‑149 (MQ=255)
                                                                                                          tATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGcc                                   >  2:701917/1‑149 (MQ=255)
                                                                                                               aTTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATttc                              <  2:672461/149‑1 (MQ=255)
                                                                                                                tttGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATttct                             >  4:139658/1‑149 (MQ=255)
                                                                                                                 ttGCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATttctt                            <  1:35534/149‑1 (MQ=255)
                                                                                                                   gCTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCa                          <  2:702548/149‑1 (MQ=255)
                                                                                                                    cTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGTTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCAt                         <  3:120248/149‑1 (MQ=255)
                                                                                                                    cTGGAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCAt                         <  2:820088/149‑1 (MQ=255)
                                                                                                                     tGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCAtt                        >  4:160116/1‑149 (MQ=255)
                                                                                                                      ggAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCAttt                       >  1:419551/1‑149 (MQ=255)
                                                                                                                      ggAGATCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCAttt                       >  3:558173/1‑149 (MQ=255)
                                                                                                                        agaTCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTAt                     >  1:359924/1‑149 (MQ=255)
                                                                                                                        agaTCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTAt                     <  1:633053/149‑1 (MQ=255)
                                                                                                                        agaTCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTAt                     <  2:667930/149‑1 (MQ=255)
                                                                                                                          aTCGTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGc                   >  3:212345/1‑149 (MQ=255)
                                                                                                                             gTCAACGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCga                >  4:401296/1‑149 (MQ=255)
                                                                                                                              tCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCgag               >  1:517463/1‑149 (MQ=255)
                                                                                                                                aaCGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGTCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGAc             >  1:153247/1‑149 (MQ=255)
                                                                                                                                aaCGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATttc                              <  2:354769/132‑1 (MQ=255)
                                                                                                                                aaCGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATttc                              >  1:354769/1‑132 (MQ=255)
                                                                                                                                 aCGCGGGAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGTAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGTGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACg            >  4:621218/1‑149 (MQ=255)
                                                                                                                                 aCGCGGGAATCTAGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACg            >  3:68085/1‑149 (MQ=255)
                                                                                                                                     gggAATCTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCAttt                                                                   <  4:41733/90‑1 (MQ=255)
                                                                                                                                       gAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACAt      >  1:571036/1‑149 (MQ=255)
                                                                                                                                          tCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGc   <  2:517463/149‑1 (MQ=255)
                                                                                                                                           cTTGCTTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCg  >  3:703868/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                               
TTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCG  >  NC_000913/1977818‑1978105

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: