Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,325,701 G→A A412V (GCT→GTT)  ftsH ← ATP‑dependent zinc metalloprotease FtsH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,325,7010GA100.0% 203.3 / NA 58A412V (GCT→GTT) ftsHATP‑dependent zinc metalloprotease FtsH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (34/24);  total (34/24)

AATATCGTTGGACGCACCGGTAGATACATGTTCCGGCCCGTAGATGATCTCTTCTGCCAGACGACCACCGTACAGCGTAGAAATCTGGCTTTCCAGTTTCTGACGGCTGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCACGAGCAGCGAACAGTGCCGCTTCGTTCACCAGGTTCGCCAGGTCAGCACCGGAGAAACCAGGAGTACCACGGGCAATGATTGCCGCGTCGATATCGGGTGCCAATGGTACGCGACGCATGTGAACTTTCAGG  >  NC_000913/3325457‑3325946
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
aaTATCGTTGGACGCACCGGTAGATACATGTTCCGGCCCGTAGATGATCTCTTCTGCCAGACGACCACCGTACAGCGTAGAAATCTGGCTTTCCAGTTTCTGACGGCTGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTc                                                                                                                                                                                                                                                  <  2:101900/250‑1 (MQ=255)
     cGTTGGACGCACCGGTAGATACATGTTCCGGCCCGTAGATGATCTCTTCTGCCAGACGACCACCGTACAGCGTAGAAATCTGGCTTTCCAGTTTCTGACGGCTGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGAttct                                                                                                                                                                                                                                            <  1:19822/251‑1 (MQ=255)
            cGCACCGGTAGATACATGTTCCGGCCCGTAGATGATCTCTTCTGCCAGACGACCACCGTACAGCGTAGAAATCTGGCTTTCCAGTTTCTGACGGCTGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCGCTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTgcg                                                                                                                                                                                                                                     <  1:267904/251‑1 (MQ=255)
            cGCACCGGTAGATACATGTTCCGGCCCGTAGATGATCTCTTCTGCCAGACGACCACCGTACAGCGTAGAAATCTGGCTTTCCAGTTTCTGACGGCTGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTgcg                                                                                                                                                                                                                                     >  1:49390/1‑251 (MQ=255)
                          cATGTTCCGGCCCGTAGATGATCTCTTCTGCCAGACGACCACCGTACAGCGTAGAAATCTGGCTTTCCAGTTTCTGACGGCTGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTTGATTCTTTCTGCGCTTCCGTCATCAcc                                                                                                                                                                                                                       >  1:258919/1‑251 (MQ=255)
                           aTGTTCCGGCCCGTAGATGATCTCTTCTGCCAGACGACCACCGTACAGCGTAGAAATCTGGCTTTCCAGTTTCTGACGGCTGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGCATCTTTCTGCGCTTCCCGCATCAcc                                                                                                                                                                                                                       >  2:228355/1‑250 (MQ=255)
                                ccGGCCCGTAGATGATCTCTTCTGCCAGACGACCACCGTACAGCGTAGAAATCTGGCTTTCCAGTTTCTGACGGCTGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGa                                                                                                                                                                                                                  <  1:150659/250‑1 (MQ=255)
                                 cGGCCCGTAGATGATCTCTTCTGCCAGACGACCACCGTACAGCGTAGAAATCTGGCTTTCCAGTTTCTGACGGCTGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGCAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACACGCTTCGTGGAAAACCCTCGTTTCTTTCTGCGCTTCAGTCATCACCTTGGaac                                                                                                                                                                                                                >  1:270168/1‑249 (MQ=255)
                                 cGGCCCGTAGATGATCTCTTCTGCCAGACGACCACCGTACAGCGTAGAAATCTGGCTTTCCAGTTTCTGACGGCTGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGc                                                                                                                                                                                                                >  1:34848/1‑251 (MQ=255)
                                  ggCCCGTAGATGATCTCTTCTGCCAGACGACCACCGTACAGCGTAGAAATCTGGCTTTCCAGTTTCTGACGGCTGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCg                                                                                                                                                                                                               >  1:249498/1‑251 (MQ=255)
                                  ggCCCGTAGATGATCTCTTCTGCCAGACGACCACCGTACAGCGTAAAAATCTGGCTTTCCAGTTTCTGACGGCTGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCCCTTCCGTCATCACCATGGTGCg                                                                                                                                                                                                               >  2:55713/1‑251 (MQ=255)
                                    cccGTAGATGATCTCTTCTGCCAGACGACCACCGTACAGCGTAGAAATCTGGCTTTCCAGTTTCTGACGGCTGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGAGCTTCCGTCATCACCATTGAGCGAc                                                                                                                                                                                                             >  1:9238/1‑251 (MQ=255)
                                       gTAGATGATCTCTTCTGCCAGACGACCACCGTACAGCGTAGAAATCTGGCTTTCCAGTTTCTGACGGCTGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGGCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGtt                                                                                                                                                                                                          <  1:13310/251‑1 (MQ=255)
                                               tctcTTCTGCCAGAAGACCACCGTACAGCGAGGAAATCTGGCTTTCCAGTTTCTGACGGCTGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCAccc                                                                                                                                                                                                  >  1:247969/1‑251 (MQ=255)
                                                        ccAGACGACCACCGTACAGCGTAGAAATCTGGCTTTCCAGTTTCTGACGGCTGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGAtt                                                                                                                                                                                         <  2:208454/251‑1 (MQ=255)
                                                              gacCACCGTACAGCGTAGAAATCTGGCTTTCCAGTTTCTGACGGCTGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCCTCATGATTTTGTCt                                                                                                                                                                                   >  1:260441/1‑251 (MQ=255)
                                                                    cGTACAGCGTAGAAATCTGGTTTTCCAGTTTCTGACGGCTGGCGCTGATTGCGTCGCCCGCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCt                                                                                                                                                                             <  2:187693/251‑1 (MQ=255)
                                                                     gTACAGCGTAGAAATCTGGCTTTCCAGTTTCTGACGGCTGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCtt                                                                                                                                                                            >  2:178633/1‑251 (MQ=255)
                                                                     gTACAGCGTAGAAATCTGGCTTTCCAGTTTCTGACGGCTGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCtt                                                                                                                                                                            >  1:124904/1‑251 (MQ=255)
                                                                      tACAGCGTAGAAATCTGGCTTTCCAGTTTCTGACGGCTGGCGCTGATTGCGTAGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCTCTTTCTGCACCGGAGCGTGTTCCGGCACCAGTCGTCCGCTAACCGCATTACCCCCTTCGTGGTAAACCGTCAATTCGTTCTGAGCTTCCGTCCTCACCATGTCGCGACGTTCCGCCCCCAGCATGCTTTTGTCTTTCCCttt                                                                                                                                                                           >  2:103644/1‑251 (MQ=255)
                                                                       aCAGCGTAGAAATCTGGCTTTCCAGTTTCTGACGGCTGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCttt                                                                                                                                                                           <  2:198975/250‑1 (MQ=255)
                                                                           cGTAGAAATCTGGCTTTCCAGTTTCTGACGGCTGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGa                                                                                                                                                                      <  1:178633/251‑1 (MQ=255)
                                                                               gAAATCTGGCTTTCCAGTTTCTGACGGCTGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTCCTCGAACTc                                                                                                                                                                  <  1:63252/251‑1 (MQ=255)
                                                                                 aaTCTGGCTTTCCAGTTTCTGACGGCTGGCGCTGATTGCGTCGCCCTCAGTCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGGCTTTACCTTTCTCGACCTCa                                                                                                                                                                 >  2:77750/1‑250 (MQ=255)
                                                                                  agcTGGCTTTCCAGTTTCTTGCGGCTGGCGCTGATTGCGGCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCCCGACCGCGTGGGATAATCGTCACTTTGTGCACCTGGTCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACTCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCTCCTTGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAAc                                                                                                                                                               <  1:28409/249‑1 (MQ=255)
                                                                                     tGGCTTTCCAGTTTCTGACGGCTGGCGCCGATTGCGTCGCACTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGATTTCTCGAACTCACCCAt                                                                                                                                                            >  2:186624/1‑251 (MQ=255)
                                                                                     tGGCTTTCCAGTTTCTGACGGCTGGAGCAAATTGCGACGCGCTCAGGCAAGAAGAAAGGCACACCCTGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGGTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTAAATCTGCGCTTCCGTCATCACCATGGAGAGGCGTTCCGCACCCAACATGATTTTGTCTTTCGCTTTCTCGAACTCCACCAt                                                                                                                                                            >  2:206755/1‑251 (MQ=255)
                                                                                                     gACGGCTGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCTAACTCAACCATCGACACAACGCGt                                                                                                                                               <  2:30513/248‑1 (MQ=255)
                                                                                                            ggtgcTGATTGCGTCGTCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACGTTGTTCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCAGCGACACAACGCGTTTGTTGCCa                                                                                                                                      <  1:159356/247‑1 (MQ=255)
                                                                                                                     tGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCCCCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCACGAGCAGCg                                                                                                                             <  2:219820/250‑1 (MQ=255)
                                                                                                                           gCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCAGGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCACGAGCAGCGAACAGTg                                                                                                                      <  2:249498/251‑1 (MQ=255)
                                                                                                                             ccTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCACGAGCAGCGAACAGTGcc                                                                                                                    >  2:56678/1‑251 (MQ=255)
                                                                                                                              cTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCACGAGCAGCGAACAGTGCCg                                                                                                                   >  1:21499/1‑251 (MQ=255)
                                                                                                                                    cAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGCCCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCACGAGCAGCGAACAGAGCCGCTTCGt                                                                                                             >  2:274551/1‑251 (MQ=255)
                                                                                                                                                  cacCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCACGAGCAGCGAACAGTGCCGCTTCGTTCACCAGGTTCGcc                                                                                               >  2:115093/1‑251 (MQ=255)
                                                                                                                                                       aGCGCACGCCCGCGTGGGATAATCGTCTCTCTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCCCCCACCATGATTTTGTCTTTCGCTTTCTCGAACGCAACCATCTACACAACGCGTTTGTTGCCACGAGCATCGAACAGTGCCGCTTCGTTCACCAGGTTCGCCAGGTc                                                                                          <  1:108752/251‑1 (MQ=255)
                                                                                                                                                              gACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCACGAGCAGCGAACAGTGCCGCTTCGTTCACCAGGTTCGCCAGGTCAGCACCg                                                                                   >  1:24665/1‑251 (MQ=255)
                                                                                                                                                                  gngTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCACGAGCAGCGAACAGTGCCGCTTCGTTCACCAGGTTCGCCAGGTCAGCACCGga                                                                                 <  2:262749/249‑1 (MQ=255)
                                                                                                                                                                  gcgTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCACGAGCAGCGAACAGTGCCGCTTCGTTCACCAGGTTCGCCAGGTCAGCACCGga                                                                                 >  1:262749/1‑249 (MQ=255)
                                                                                                                                                                  gcgTGGGATAATCGTCACTGTGTGCTCCGTATTGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCCTCACCATTGTGCGACGTTCCGCACCCATCATGATTTTGTCCTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCACGAGCAGCGATCAGTGCCGCTTCGTTCACCAGGTTCGCCAGGTCAGCACCGgaga                                                                               <  2:271663/251‑1 (MQ=255)
                                                                                                                                                                   cgTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCACGAGCAGCGAACAGTGCCGCTTCGTTCACCAGGTTCGCCAGGTCAGCACCGGAGaa                                                                              <  2:56171/251‑1 (MQ=255)
                                                                                                                                                                     tGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATACTCGCATGACCCGCTTCGTGGTAAACCGTAGATTCTTTCTGCGCTTCCGTCATCCCCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCACGAGCAGCGAACAGTGCCGCTTCGTTCACCCGGTTCGCCAGGTCAGCACCGGAGAAAc                                                                            >  1:255723/1‑251 (MQ=255)
                                                                                                                                                                         aTAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCACGAGCAGCGAACAGTGCCGCTTCGTTCACCAGGTTCGCCAGGTCAGCACCGGAGAAACCAgg                                                                        >  1:177740/1‑251 (MQ=255)
                                                                                                                                                                                  aCTTTGTGCACCGCATCGCGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGATTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACAATGGAGCGACCTTCCGCACCCATCATGATTATGTCTTTCGCTTTCTCGCACTAAACCATCGACACAACGCGTTTGTTGCCACGAGCAGCGAACAGTGCCGCTTCGTTCACCAGGTTCGCCAGGTCAGCACAGGCGAAACCAGGAGTACNACg                                                               >  2:260602/1‑251 (MQ=255)
                                                                                                                                                                                          cACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCca                                                                                                                                                                                                 >  1:257875/1‑113 (MQ=255)
                                                                                                                                                                                          cACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCGGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCca                                                                                                                                                                                                 <  2:257875/113‑1 (MQ=255)
                                                                                                                                                                                            ccGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCACGAGCAGCGAACAGTGCCGCTTCGTTCACCAGGTTCGCCAGGTCAGCACCGGAGAAACCAGGAGTACCACGGGCAATGAtt                                                     <  1:104042/251‑1 (MQ=255)
                                                                                                                                                                                                     gtTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCACGAGCAGCGAACAGTGCCGCTTCGTTCACCAGGTTCGCCAGGTCAGCACCGGAGAAACCAGGAGTACCACGGGCAATGATTGCCGCGTCg                                            >  1:27411/1‑251 (MQ=255)
                                                                                                                                                                                                               cAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCACGAGCAGCGAACAGTGCCGCTTCGTTCACCAGGTTCGCCAGGTCAGCACCGGAGAAACCAGGAGTACCACGGGCAATGATTGCCGCGTCGATATCGGGTg                                  >  2:110517/1‑251 (MQ=255)
                                                                                                                                                                                                                 ggCGACCGATAATCGCATGACCCGCTTCGTTGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCCCCATGGAGTGTCGTTCCGCTCCCATCAGGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCCCGAGCAGCGAACAGTGCCGCTTCGTTCACCAGGTTCGCCAGGTCAGCACCGGAGAAACCAGGATTACCACGGTCAATGATTGCCGCGTCGATATCGGGTGcc                                <  1:206755/251‑1 (MQ=255)
                                                                                                                                                                                                                       cGATAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCACGAGCAGCGAACAGTGCCGCTTCGTTCACCAGGTTCGCCAGGTCAGCACCGGAGAAACCAGGAGTACCACGGGCAATGATTGCCGCGTCGATATCGGGTGCCAATGGt                          <  2:34848/251‑1 (MQ=255)
                                                                                                                                                                                                                         aTAATCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCACGAGCAGCGAACAGTGCCGCTTCGTTCACCAGGTTCGCCAGGTCAGCACCGGAGAAACCAGGAGTACCACGGGCAATGATTGCCGCGTCGATATCGGGTGCCAATgnaac                        >  2:257910/1‑247 (MQ=255)
                                                                                                                                                                                                                           aaTCGCATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCACGAGCAGCGAACAGTGCCGCTTCGTTCACCAGGTTCGCCAGGTCAGCACCGGAGAAACCAGGAGTACCACGGGCAATGATTGCCGCGTCGATATCGGGTGCCAATGGTAcgc                      <  1:193181/251‑1 (MQ=255)
                                                                                                                                                                                                                                cATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCACGAGCAGCGAACAGTGCCGCTTCGTTCACCAGGTTCGCCAGGTCAGCACCGGAGAAACCAGGAGTACCACGGGCAATGATTGCCGCGTCGATATCGGGTGCCAATGGTACGCGACGc                 >  2:185120/1‑251 (MQ=255)
                                                                                                                                                                                                                                cATGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCACGAGCAGCGAACAGTGCCGCTTCGTTCACCAGGTTCGCCAGGTCAGCACCGGAGAAACCAGGAGTACCACGGGCAATGATTGCCGCGTCGATATCGGGTGCCAATGGTACGCGACGc                 >  2:162042/1‑251 (MQ=255)
                                                                                                                                                                                                                                 aTGACCCGCTTCGTGGTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGCCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCACGAGCAGCGAACAGTGCCGCTTCGTTCACCAGGTTCGCCAGGTCAGCACCGGAGAAACCAGGAGTACCACGGGCAATGATTGCCGCGTCGATATCGGGGGCCAATGGTACGCGACGc                 >  1:51021/1‑250 (MQ=255)
                                                                                                                                                                                                                                       cGCTTCGTGGTAAAACGTCGATTCTTTATGCGCTTACGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCACGAGCAGCGAACAGTGCCGCTTCGTTCACCAGGTTCGCCAGGTCAGCACCGGAGAAACCAGGAGTACCACGGGCAATGATTGCCGCGTCGATATCGGGTGCCAATGGTACGCGACGCATGNGa           >  2:254119/1‑250 (MQ=255)
                                                                                                                                                                                                                                                gTAAACCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCACGAGCAGCGAACAGTGCCGCTTCGTTCACCAGGTTCGCCAGGTCAGCACCGGAGAAACCAGGAGTACCACGGGCAATGATTGCCGCGTCGATATCGGGTGCCAATGGTACGCGACGCATGTGAACTTTCAgg  >  1:187083/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
AATATCGTTGGACGCACCGGTAGATACATGTTCCGGCCCGTAGATGATCTCTTCTGCCAGACGACCACCGTACAGCGTAGAAATCTGGCTTTCCAGTTTCTGACGGCTGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCACGAGCAGCGAACAGTGCCGCTTCGTTCACCAGGTTCGCCAGGTCAGCACCGGAGAAACCAGGAGTACCACGGGCAATGATTGCCGCGTCGATATCGGGTGCCAATGGTACGCGACGCATGTGAACTTTCAGG  >  NC_000913/3325457‑3325946

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