Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,417,140 C→T P63P (CCC→CCT yjfP → carboxylesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,417,1400CT100.0% 230.0 / NA 66P63P (CCC→CCTyjfPcarboxylesterase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (35/31);  total (35/31)

AGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTCGCCACCCCACGGTGAGATGTACCGCCAGCATGATGGGATCGGGCTATTTTACA  >  NC_000913/4416895‑4417383
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   
agcagAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCa                                                                                                                                                                                                                                                <  1:49413/251‑1 (MQ=255)
       aCTTGTGCTAAGGCTTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTACCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGc                                                                                                                                                                                                                                         >  1:188840/1‑251 (MQ=255)
       aCTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGc                                                                                                                                                                                                                                         <  1:642356/251‑1 (MQ=255)
             gCTAAGGATTTCTTCATGTCATACCAGACAAAAATGATAGACAGATGATTGTAATAGAATCACGCGAGCTGGCAGATATGCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATATCTATTTTGCCGTTGCGCTGGCGCATGCTGTTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGtttt                                                                                                                                                                                                                                   <  1:83881/251‑1 (MQ=255)
               tAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGtttt                                                                                                                                                                                                                                   >  1:447255/1‑249 (MQ=255)
                aaGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTa                                                                                                                                                                                                                                  <  1:294668/249‑1 (MQ=255)
                   gATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAgtggt                                                                                                                                                                                                                             <  1:617113/251‑1 (MQ=255)
                        cttcATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGc                                                                                                                                                                                                                        >  1:23289/1‑251 (MQ=255)
                                        aCAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCGGGCAGATATTCCCGTTCTTCATGCGTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCGGGTGTATAGCTATTTTGCCGTTGCGCGGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCTCGGTAGCCGGGTTAGTGGTGACGCAGCGCGGCGGTTAAAt                                                                                                                                                                                                        <  1:30820/251‑1 (MQ=255)
                                          aaaaaGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCa                                                                                                                                                                                                      <  1:103150/251‑1 (MQ=255)
                                           aaaaGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCa                                                                                                                                                                                                      >  1:364774/1‑250 (MQ=255)
                                              aGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTc                                                                                                                                                                                                  <  2:123854/251‑1 (MQ=255)
                                                         atgatTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCtt                                                                                                                                                                                       >  1:207927/1‑251 (MQ=255)
                                                                aaaTAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCCAAACTTGCTACa                                                                                                                                                                                 >  1:66178/1‑250 (MQ=255)
                                                                 aaTAGAATCACGCGAGCTGTCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACaaa                                                                                                                                                                               <  1:246985/251‑1 (MQ=255)
                                                                 aaTAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACaaa                                                                                                                                                                               <  2:647271/251‑1 (MQ=255)
                                                                    aGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGt                                                                                                                                                                             >  1:407508/1‑250 (MQ=255)
                                                                       aTCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATgg                                                                                                                                                                         >  2:500557/1‑250 (MQ=255)
                                                                           cgcgAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAgg                                                                                                                                                                      >  1:129224/1‑250 (MQ=255)
                                                                             cgAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAAt                                                                                                                                                                   <  2:70965/251‑1 (MQ=255)
                                                                              gAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCCAATCTTGCTTCAAAGTATGCAGGAAtt                                                                                                                                                                  >  1:357640/1‑251 (MQ=255)
                                                                              gAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAAtt                                                                                                                                                                  <  1:361743/251‑1 (MQ=255)
                                                                              gAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAAtt                                                                                                                                                                  <  2:332186/251‑1 (MQ=255)
                                                                                   ggCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTACCTCGGCTTTACTTCATCCATTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGAGCATGCCGGATGCGCCTGAGCACGGTAGCCGTTTTAGTGGGGACGCAGCGCGGCGGTTAAATCAATGCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCacg                                                                                                                                                              <  2:549763/250‑2 (MQ=255)
                                                                                                            tATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCga                                                                                                                                       >  1:221821/1‑248 (MQ=255)
                                                                                                               ccTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGAGCCTGATCACGGTAGCCGTTTTAGTGGGGACGCAGCGCGGCGGGTAAATCAATTTTGGCAAATCTTGCTACAAAGGATGCAGGGATTTACTACTTTAAGTGCGGCAATAGCCGAAGaaa                                                                                                                                  >  1:555223/1‑250 (MQ=255)
                                                                                                                        cAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCa                                                                                                                                                                                             <  2:156640/182‑1 (MQ=255)
                                                                                                                        cAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCa                                                                                                                                                                                             >  1:156640/1‑182 (MQ=255)
                                                                                                                          aaaaGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGGATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTtg                                                                                                                      >  1:601714/1‑251 (MQ=255)
                                                                                                                               aTACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAg                                                                                                                                                                                                                      >  1:463593/1‑150 (MQ=255)
                                                                                                                               aTACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAg                                                                                                                                                                                                                      <  2:463593/150‑1 (MQ=255)
                                                                                                                                    ccGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAATATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCt                                                                                                            >  1:416629/1‑251 (MQ=255)
                                                                                                                                     cGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCCATTATGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTTCGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTg                                                                                                           >  2:189905/1‑251 (MQ=255)
                                                                                                                                      gTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAAATGGCTGGTTGATGACCGTCTgg                                                                                                          >  2:187720/1‑251 (MQ=255)
                                                                                                                                        tACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATACATTCTGGCAAATCTTGCTACAAAGAATGCAGGAATTTAATACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTngg                                                                                                         >  2:603187/1‑249 (MQ=255)
                                                                                                                                          ccGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGt                                                                                                      >  2:469201/1‑251 (MQ=255)
                                                                                                                                            gTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTCGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCCATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCg                                                                                                    >  2:479359/1‑251 (MQ=255)
                                                                                                                                                gTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTgg                                                                                                <  1:537865/251‑1 (MQ=255)
                                                                                                                                                  aaTTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGgc                                                                                               <  2:348072/250‑1 (MQ=255)
                                                                                                                                                          aTCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGAt                                                                                       <  2:364774/250‑1 (MQ=255)
                                                                                                                                                                       ttCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACgg                                                                         >  1:108349/1‑251 (MQ=255)
                                                                                                                                                                          aTCCAGTCGGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTTGTTTGCGGGTGATCATGCCGGATGCGCCTGATCGCGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTTGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGTCAc                                                                      <  2:188840/251‑1 (MQ=255)
                                                                                                                                                                           tCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTTCTACAAAATATTCAg                                                                                                                                                                       >  1:430069/1‑153 (MQ=255)
                                                                                                                                                                           tCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAg                                                                                                                                                                       <  2:430069/153‑1 (MQ=255)
                                                                                                                                                                           tCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGGATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAAATGGCTGCTTGATGACCGTCTTGCAGTCGGTGGCGCGTCGATGGGCGCGGTGACGGCAcc                                                                     >  2:251589/1‑250 (MQ=255)
                                                                                                                                                                            ccAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACt                                                                     <  2:221821/250‑1 (MQ=255)
                                                                                                                                                                              aGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAACTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTTGGTGGCGCGTCGATGGGCGCGATGACGGCACTggg                                                                  >  1:423257/1‑251 (MQ=255)
                                                                                                                                                                              aGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGGTAAATCAATTCTGGCAAATATTGCTACAAAGAATGCAGGAATTTAATACTTTACGTGGGGCAATAGGCGAAGAAAAATGGCTGATTGATGAACGTATGGCAGTCGGTGGGGCGTCGAGGGGCGCGCTGACGGCACTggg                                                                  >  1:582285/1‑251 (MQ=255)
                                                                                                                                                                                      gtATAGCTATTTTGCAGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGGATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGGAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTAccc                                                          >  1:483909/1‑250 (MQ=255)
                                                                                                                                                                                           gCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAAACAATTCTGGCAAATATTGCTAACAAGGATTAAGGAATTCACTACCTTACGTGAGGCAATATCCGAAGAAAACTGGCTGCTTGTTGAACGTCTTGCAGTCGGTGGCGCGTCGATGGGGGGGATGATGGcgcctggggttaccccccgc                                                     >  2:616839/1‑232 (MQ=255)
                                                                                                                                                                                                 ttGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATgg                                                                                     <  1:367203/213‑1 (MQ=255)
                                                                                                                                                                                                   gCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTCGCCACCCCAc                                             >  1:459468/1‑251 (MQ=255)
                                                                                                                                                                                                      gTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCCATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCCCTGGGGGTTACCGCTCGCCACCCCCCggg                                          >  1:100307/1‑250 (MQ=255)
                                                                                                                                                                                                          cgcTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGGCCTGGGGATTACCGCTCGCCACCCCACGGTgagg                                      >  2:183309/1‑250 (MQ=255)
                                                                                                                                                                                                           gcTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGGGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTCGCCACCCCACGGTg                                         >  1:475867/1‑247 (MQ=255)
                                                                                                                                                                                                             tGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGAGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTTCTACAAAGTATGCAGGAATTCACTACTTTACGGGCGGCAATAGCCTAAGAAAACTGGCTGCTTGATTACCGTCTGGTAGTCTGTGGCGCGGCGAGGGGCGCGATGACGGCACTGGGGATTAACGATCGCCAACCCACTGTGAGgtgt                                   >  1:202653/1‑251 (MQ=255)
                                                                                                                                                                                                                 gcAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTCGCCACCCCACGGTGAGATGTACCg                               <  1:124952/251‑1 (MQ=255)
                                                                                                                                                                                                                   aGGNTGGTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTCGCCACCCCACGGTGAGATGTACCGcc                             <  2:601714/251‑1 (MQ=255)
                                                                                                                                                                                                                        ggTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCCCTGGGGATTACCGCTCGCCACCCCACGGTGAGATGTACCGCCAGCat                        >  1:114690/1‑251 (MQ=255)
                                                                                                                                                                                                                        ggTTTGCGGGTGATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTCGCCACCCCACGGTGAGATGTACCGCCAGCat                        <  1:603187/251‑1 (MQ=255)
                                                                                                                                                                                                                           ttGCGGGTGGTCGTGCCGGATGCGCTTGATCCCGGTAGCCGTTTTAGGGGGGAAGCAGCGCGGCGGTTAAATAAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTATTGGCCGTCTGCCAGTCGGTGGGGGTTCGATGGGCGCGATGACGGCACTGGGGATTACCTCTCGCCTCCCCACGGTGAGATGTACCGCCAGCatga                      <  2:207927/250‑1 (MQ=255)
                                                                                                                                                                                                                                   gATCATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGAGGGGCGCGATGACGGCACTGGGGATCACCGCTCGCCCCCCCACGGGGAGATGAACCGCCAGCATGATGGGATCgg             <  1:187720/251‑1 (MQ=255)
                                                                                                                                                                                                                                      cATGCCGGATGCGCCTGATCCCGGTAGCCGTTTTAGTGGTGACGCAGCGCTGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTCGCCACCCCACGGTGAGATGTACCGCCAGCATGATGGGATCGGGCt          <  1:528760/251‑1 (MQ=255)
                                                                                                                                                                                                                                      cATGCCGGATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTGCCGCTCGCCACCCCACGGTGAGATGTACCGCCAGCATGATGGGATCGGGCt          <  2:407508/251‑1 (MQ=255)
                                                                                                                                                                                                                                             gATGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAAACAATTCTGGCAAATCTTTCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAACCGAAGAAAACTGGCTGCTTTATGACCGTCTGGCAGTCCGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTCGCCACCCCACGGGGAGATGTACCGCCAGCATGATGGGATCGGGCTATTTTa    >  1:262045/1‑250 (MQ=255)
                                                                                                                                                                                                                                              aTGCGCCTGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTCGCCACCCCACGGTGAGATGTACCGCCAGCATGATGGGATCGGGCTATTTTACa  <  1:506141/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   
AGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATACCAGACAAAAAGGAGAGACAGATGATTGAAATAGAATCACGCGAGCTGGCAGATATTCCCGTTCTTCATGCTTATCCTGTCGGGCAAAAAGATACCCCGTTACCGTGCGTAATTTTTTATCACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTCGCCACCCCACGGTGAGATGTACCGCCAGCATGATGGGATCGGGCTATTTTACA  >  NC_000913/4416895‑4417383

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