Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 727,117 C→T G539D (GGC→GAC)  kdpA ← K(+) transporting P‑type ATPase subunit KdpA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913727,1170CT98.3% 194.9 / NA 59G539D (GGC→GAC) kdpAK(+) transporting P‑type ATPase subunit KdpA
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (24/33);  total (25/34)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GCTTGCCATCGCGATGCTAATACAGGTGGTCAGCAGACTGCCGATCCAGACGATAAACATCACCGGATTGCGCCATTGCGCCTGCGGGTTTAATTTTTTCACCGCTTCTTTCAGCGCCTGAACGACAAGTGTTGGTTCGAATAGCGCCAGTTGTTTACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGCGCCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTACTCACCAGCGAACCGGCAATTGCCATCACCGGGATAATCACCCCGAAGCGACCGACAAACATGCAGAACGCCAGTAAACAGTTCCAGAACGGAGAGTTGGCGCTTAATCCGGCAAAGGCGCTGCCGTTGTTGTTAGCGGCGGATGA  >  NC_000913/726872‑727352
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                           
gCTTGCCATCGCGATGCTAATACAGGTGGTCAGCAGACTGCCGATCCAGACGATAAACATCACCGGATTGCGCCATTGCGCCTGCGGGTTTAATTTTTTCACCGCTTCTTTCAGCGCCTGAACGACAAGTGTTGGTTCGAATAGCGCCAGTTGTTTACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAAcc                                                                                                                                                                                                                                        <  2:184994/251‑1 (MQ=255)
gCTTGCCATCGCGATGCTAATACAGGTGGTCAGCAGACTGCCGATCCAGACGATAAACATCACCGGATTGCGCCATTGCGCCTGCGGGTTTAATTTTTTCACCGCTTCTTTCAGCGCCTGAACGACAAGTGTTGGTTCGAATAGCGCCAGTTGTTTACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAAcc                                                                                                                                                                                                                                        >  2:644973/1‑251 (MQ=255)
   tGCCATCGCGATGCTAATACAGGTGGTCAGCAGACTGCCGATCCAGACGATAAACATCACCGGATTGCGCCATTGCGCCTGCGGGTTTAATTTTTTCACCGCTTCTTTCAGCGCCTGAACGACAAGTGTTGGTTCGAATAGCGCCAGTTGTTTACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCa                                                                                                                                                                                                                                       <  2:223532/249‑1 (MQ=255)
      cATCGCGATGCTAATACAGGTGGTCAGCAGACTGCCGATCCAGACGATAAACATCACCGGATTGCGCCATTGCGCCTGCGGGTTTAATTTTTTCACCGCTTCTTTCAGCGCCTGAACGACAAGTGTTGGTTCGAATAGCGCCAGTTGTTTACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAac                                                                                                                                                                                                                                  <  2:12753/251‑1 (MQ=255)
       aTCGCGATGCTAATACAGGTGGTCAGCAGACTGCCGATCCAGACGATAAACATCACCGGATTGCGCCATTGCGCCTGCGGGTTTAATTTTTTCACCGCTTCTTTCAGCGCCTGAACGACAAGTGTTGGTTCGAATAGCGCCAGTTGTTTACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAac                                                                                                                                                                                                                                  >  2:103871/1‑250 (MQ=255)
            gATGCTAATACCGGTGGTCAGCAGACTGCCGATCCAGACGATAAACATCACCGGATTGCGCCATTGCGCCTGCGGGTTTAATTTTTTCACCGCTTCTTTCAGCGCCTGAACGACAAGTGTTGGTTCGAATAGCGCCAGTTGTTTACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTg                                                                                                                                                                                                                            >  2:276869/1‑251 (MQ=255)
                 tAATACAGGTGGTCAGCAGACTGCCGATCCAGCCGATAAACATCACCGGATTGCGCCAATGCGCCTGCGGGATTAATTTTTTCACCGCTTATTTCAGCGCCTGAACGACAAGGGCTGGTACGAATAGCGCCAGTTGTTTACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGATAAGAGGGTCAACCAGCAACACggggc                                                                                                                                                                                                                           >  2:604345/1‑244 (MQ=255)
                  aaTACAGTTGGTCAGCAGACTGCCGATCCAGACGATAAACATCACCGGATTGCGCCATTGCGCCTGCGGGTTTAATTTTTTCACCGCTTCTTTCAGCTCCTGAACGACAAGTGTTGGTTCGAATAGCGCCAGTTGTTTACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACTGTGCCGAtt                                                                                                                                                                                                                      <  2:353121/251‑1 (MQ=255)
                               agcggACTGCCGATCCAGACTCTATACATCCCCGCATTTCGCCATTGCGCCTGCGGCTTTAATTTTTTCACCGCTTCTTTCAGCGCCTGAACGACAAGTGTTGGTTCGAATAGCGCCAGTTGTTTACGACTCATCTTCAGTGCTCACTCCATATCATCAGGCGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCCGTGTGTCACCCAGCAACACGGTGCCGTTTAACAGGCcagcaa                                                                                                                                                                                                         <  2:19180/247‑1 (MQ=255)
                                  agACTGCCGATCCAGACGATAAACATCACCGGATTGCGCCATTGCGCCTGCGGGTTTAATTTTTTCACCGCTTCTTTCAGCGCCTGAACGACAAGTGTTGGTTCGAATAACGCCAGTTGTGTACGACTCATATTCAGTGCTCACTCACTATCATCAGGAGAGGTATTCCGCCACCGGACCAAGCACCAGGGCAGGGATAAAGGGCAGGGGGTCAACCAGCAACACGGGGCCGATTAACAGGCCaacaaac                                                                                                                                                                                                       >  1:475529/1‑250 (MQ=255)
                                       gCCGATCCAGACGATAAACATCACCGGATTGCGCCATTGCGCCTGCGGGTTTAATTTTTTCACCGCTTCTTTCAGCGCCTGAACGACAAGTGTTGGTGCGAATAGCGCCAGTTGTTTACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCg                                                                                                                                                                                                   <  1:542937/249‑1 (MQ=255)
                                       gCCGATCCAGACGATAAACATCACCGGATTGCGCCATTGCGCCTGCGGGTTTAATTTTTGCACCGCTTCTTTCAGCGCCTGAACGACAAGTGTTGGTTCGAATAGCGCCAGTTGTTTACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCg                                                                                                                                                                                                   <  1:270799/249‑1 (MQ=255)
                                                gACGATAAACATCACCGGATTGCGCCATTGCGCCTGCGGGTTTAATTTTTTCACCGCTTCTTTCAGCGCCTGAACGACAAGTGTTGGTTCGAATAGCGCCAGTTGTTTACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGtt                                                                                                                                                                                        <  1:153322/251‑1 (MQ=255)
                                                  cGATAAACATCACCGGATTGCGCCATTGCGCCTGCGGGTTTAATTTTTTCACCGCTTCTTTCAGCGCCTGAACGACAAGTGTTGGTTCGAATAGCGCCAGTTGTTTACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTgg                                                                                                                                                                                      >  1:374645/1‑251 (MQ=255)
                                                   gATAAAAATCGCCGGATTTCGCCATTGCGTCTGCGGGTTTAAGTTTTTCACCGCTTCTTTCAGCGCCTGAACGGCAAGTGTTGGTTCGAATAGCGCCAGTTGTTTACGACTCATAGTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACTGTGCCGATTAACAGGCCATCAAACAGCGGGCCGTGCGTTGGc                                                                                                                                                                                     <  2:421315/251‑1 (MQ=255)
                                                          aTCACCGGATTGCGCCATTGCGCCTGCGGGTTTAATTTTTTCACCGCTTCTTTCAGCGCCTGAACGACAAGTGTTGGTTCGAATAGCGCCAGTTGTTTACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGc                                                                                                                                                                              <  1:47691/251‑1 (MQ=255)
                                                               cGGATTGCGCCATTGCGCCTGCGGGTTTAATTTTTTCACCGCTTCTTTCAGCGCCTGAACGACAAGTGTTGGTTCGAATAGCGCCAGTTGTTTACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCgg                                                                                                                                                                           >  2:475634/1‑249 (MQ=255)
                                                                                          tANTTTTTTCACCGCTTCTTTCAGCGCCTGAACGACAAGTGTTGGTTCGAATAGCGCCAGTTGTTTACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATGCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAGCCAGGAACACGTTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCGGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTACTCAc                                                                                                                                               <  2:603369/250‑1 (MQ=255)
                                                                                             ttttttCACCGCTTCTTTCAGCGCCTGAACGACAAGTGTTGGTTCGAATAGCGCCAGTTGTTTACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTACTCACCAGc                                                                                                                                           >  1:237813/1‑251 (MQ=255)
                                                                                              tttttCACCGCTTCTTTCAGCGCCTGAACGACAAGTGTTGGTTCGAATAGCGCCAGTTGTTTACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTACTCACCAgt                                                                                                                                           <  2:225223/250‑2 (MQ=255)
                                                                                                 ttCACCGCTTCTTTCAGCGCCTGAACGACAAGTGTTGGTTCGAATAGCGCCAGTTGTTTACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATGCCGCCACCGGACCAAGCGCCAGGGTAGGGATAAAGGTCAGTGGGCCAACCAGCAACACGGTGCCGATTAACAGGCCCACAAACAGCGGGCCGTGCGTTGGCAGCGTGCAGGAGCTGGCGGCTTGGCTCTTTTTACTCACCAGCGaa                                                                                                                                        >  2:250138/1‑250 (MQ=255)
                                                                                                   cACCGCTTCTTTCAGCGCCTGAACGACAAGTGTTGGTTCGAATAGCGCCAGTTGTTTACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCt                                                                                                                                                                <  1:472330/224‑1 (MQ=255)
                                                                                                   cACCGCTTCTTTCAGCGCCTGAACGACAAGTGTTGGTTCGAATAGCGCCAGTTGTTTACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCt                                                                                                                                                                >  2:472330/1‑224 (MQ=255)
                                                                                                       gCTTCTTTCAGCGCCTGAACGACAAGTGTTGGTTCGAATAGCGCCAGTTGTTTACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTACTCACCAGCGAACCGGCaa                                                                                                                                 <  1:444849/251‑1 (MQ=255)
                                                                                                           cttTCAGCGACTGAACGACAAGTGTTGGTTCGAATAGAGCCAGTTGTTTACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTACTCACCAGCGAACCGGCAATTGc                                                                                                                             >  2:653748/1‑251 (MQ=255)
                                                                                                               cAGCGCCTGAACGACAAGTGTTGGTTCGAATAGCGCCAGTTGTTTACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTACTCACCAGCGAACCGGCAATTGCCAt                                                                                                                          <  1:432628/250‑1 (MQ=255)
                                                                                                                       gAACGACAAGTGTTGGTTCGAATAGCGCCAGTTGTTTACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTACTCACCAGCGAACCGGCAATTGCCATCACCGGGAt                                                                                                                 <  1:151706/251‑1 (MQ=255)
                                                                                                                           gACAAGTGTTGGTTCGAATAGCGCCAGTTGTTTACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTACTCACCAGCGAACCGGCAATTGCCATCACCGGGATAATc                                                                                                             >  1:419105/1‑251 (MQ=255)
                                                                                                                             cAAGTGTTGGTTCGAATAGCGCCAGTTGTTTACGACTAATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGACAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTACTCACCAGCGAACCGGCAATTGCCATCACCGGGATAATCAc                                                                                                           >  2:558249/1‑251 (MQ=255)
                                                                                                                              aaGTGTTGATTCGAATAGCGCCAGTTGTTTACGACCCATATTCAGTGCTCACTCACTATCATCAGGAGAGATCTTCCGCCACCGGACCAAGCGCCAGGGCAGGGTTCAAGGTCAGTGTGTCAACCCGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGGGCGGGCGCCCTGGTGCTTTTTACTCGCCAGCGAACCGGCAATTGCCATTACCGGCATACTCAcc                                                                                                          >  2:643830/1‑251 (MQ=255)
                                                                                                                              aaGTATTGGTTCGAATAGCGCCAGTTGTTTACGACTCATATTCAGTGCTCGCTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGAGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCAGTGCGTTGGCAACGTGCCGGAGCAGGCGGCTTGGCTCTTTTTACTCACCAGCGAACCGGCAATTGCCATCACCCGGCTAATCAcc                                                                                                          >  1:77546/1‑251 (MQ=255)
                                                                                                                                   ttGGTTCGAATAGCGCCAGTTGTTTACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTACTCACCAGCGAACCGGCAATTGCCATCACCGGGATAATCACCCCGa                                                                                                      >  1:50457/1‑250 (MQ=255)
                                                                                                                                      gTTCGAATAGCGCCAGTTATTTACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTACTCACCAGCGAACCGGCAATTGCCATCACCGGGATACTCACCCCGAAGcg                                                                                                  <  2:643833/251‑1 (MQ=255)
                                                                                                                                       ttCGAATAGCGCCAGTTGTTTACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTACTCACCAGCGAACCGGCAATTGCCATCACCGGGATAATCACCCCGAAGcg                                                                                                  >  2:489215/1‑250 (MQ=255)
                                                                                                                                               gcgcCAGTTGTTTACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATAGTCCGCCCCCGGCCCAAGCGCCAGGGGAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTCCTCACCAGCGAACCGGCAATTGCCATCACCGGGATAATCACCCCGAAGCGACCGACaaa                                                                                         <  1:296901/251‑1 (MQ=255)
                                                                                                                                                cgcCAGTTGTTTACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTACTCACCAGCGAACCGGCAATTGCCATCACCGGGATAATCACCCCGAAGCGACCGACAAAc                                                                                        <  1:103871/251‑1 (MQ=255)
                                                                                                                                                    aGTTGTTTACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTACTCACCAGCGAACCGGCAATTGCCATCACCGGGATAATCACCCCGAAGCGACCGACAAACATGc                                                                                    <  1:644973/251‑1 (MQ=255)
                                                                                                                                                          ttACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTACTCACCAGCGAACCGGCAATTGCCATCACCGGGATAATCACCCCGAAGCGACCGACAAACATGCAGAAc                                                                               <  1:598734/250‑1 (MQ=255)
                                                                                                                                                          ttACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTACTCACCAGCGAACCGGCAATTGCCATCACCGGGATAATCACCCCGAAGCGACCGACAAACATGCAGAAc                                                                               <  2:374645/250‑1 (MQ=255)
                                                                                                                                                                 tCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTACTCACCAGCGAACCGGCAATTGCCATCACCGGGATAATCACCCCGAAGCGACCGACAAACATGCAGAACGCCAGTa                                                                        >  1:494648/1‑250 (MQ=255)
                                                                                                                                                                     atTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCACCCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTACTCACCAGCGAACCGGCAATTGCCATCACCGGGATAATCACCCCGAAGCGACCGACAAACATGCAGAACGCCAGTAAAc                                                                     >  2:97994/1‑249 (MQ=255)
                                                                                                                                                                        cAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTACTCACCAGCGAACCGGCAATTGCCATCACCGGGATAATCACCCCGAAGCGACCGACAAACATGCAGAACGCCAGTAAACAGTTc                                                                <  1:54326/251‑1 (MQ=255)
                                                                                                                                                                         aGTGCTCACTCAATATCATCAGGATAGATATTCCGCCTCCGGACCAAGCGCCAGGGCAGGGATAAAGTGCAGTGTGTCAACCAGCAACCCGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTACTCCCCAGCGAACCTGCAATTGCCATCACCGGGTTAATCACCCCGAAGCGACCGACAAACATGCATAACGCCAGTAAACAGTTcc                                                               <  1:2494/251‑1 (MQ=255)
                                                                                                                                                                              tcactcaATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTACTCACCAGCGAACCGGCAATTGCCATCACCGGGATAATCACCCCGAAGCGACCGACAAACATGCAGAACGCCAGTAAACAGTTCCAGAAc                                                          >  2:654938/1‑251 (MQ=255)
                                                                                                                                                                               cactcaATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTACTCACCAGCGAACCGGCAATTGCCATCACCGGGATAATCCCCCCGAAGCGACCGACAAACATGCAGAACTCCAGTAAACAGTTCCAGAAc                                                          <  1:250138/250‑1 (MQ=255)
                                                                                                                                                                                   caATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTACTCACCAGCGAACCGGCAATTGCCATCACCGGGATAATCACCCCGAAGCGACCGACAAACATGCAGAACGCCAGTAAACAGTTCCAGAACgg                                                        >  2:105325/1‑248 (MQ=255)
                                                                                                                                                                                       atcatcAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCTGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTACTCCCCAGCGAACCGGCAATTGCCATCACCGGGATAATCACCCCGAAGCGACCGACAAACATGCAGAACGCCAGTAAACAGTTCCAGAACGGAGAGtt                                                  <  1:127984/250‑1 (MQ=255)
                                                                                                                                                                                       atcatcAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTACTCACCAGCGAACCGGCAATTGCCATCACCGGGATAATCACCCCGAAGCGACCGACAAACATGCAGAACGCCAGTAAACAGTTCCAGAACGGAGAGtt                                                  <  1:487643/250‑1 (MQ=255)
                                                                                                                                                                                                      atTCCGCCCCCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTACTCACCAGCGAACCGGCAATTGCCATCACCGGGATAATCACCCCGAAGCGACCGACAAACATGCAGAACGCCAGTAAACAGTTCCAGAACGGAGAGTTGGCGCTTAATCCGGCa                                  <  2:412549/251‑1 (MQ=255)
                                                                                                                                                                                                             cACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGGTGGCAGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTACTCACCCGCGAACCGGCAATTGCCATCACCGGGATAATCACCCCGAAGCGCCCGACACACATGCAGAACGCCAGTAAACCGTTCCAGAACGGAGAGATGGCGCTTAATCCGGCAAAGgcgc                           >  2:310104/1‑251 (MQ=255)
                                                                                                                                                                                                               cctgACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCTTTGCTCTTTTTACTCACCAGCGAACCGGCAATTGCCATCACCGGGATAATCACCCCGAAGCGACCGACAAACATGCAGAACGCCAGTAAACAGTTCCAGAACGGAGAGTTGGCGCTTAATCCGGCAAAGGCGCTg                         <  1:294113/248‑1 (MQ=255)
                                                                                                                                                                                                                cGGACCAAGCGCCAGGGCTGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGCGCCGGAGCTGGCGGCTTGGCTCTTTTTACTCACCAGCGAACCGGCAATTGCCATCACCGTGAAAATTACCCCGAAGCGACCGACCAACATGCAGAATGCCCGTAAACAGTTCCAGAACGGCGAGTTGGCGCTTAATCCTGCAAAGGCGCTGc                        >  2:223092/1‑251 (MQ=255)
                                                                                                                                                                                                                  gACCAAGCGCCATGGCATGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGAGTGCCGGAGCTGGCGGCTTGGCTCTTTTTACTCACCAGCGAACCGGCAATTGCCATCACCGGGATAATCACCCCGAAGCGACCGACAAACATGCAGAACGCCAGTAAACAGTTCCAGAACGGAGAGTTGGCGCTTAATCCGGCAAAGGCGCTGCCg                      <  2:725/251‑1 (MQ=255)
                                                                                                                                                                                                                      aaGCGCCAGGGCAGGGATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTAATCACCAGCGAACCGGCAATTGCCATCACCGGGATAATCACCCCGAAGCGACCGACAAAAATGCAGAACGCCAGTAAACAGTTCCAGAACGGAGAGTTGGCGCTTAATCCGGCAAAGGCGCTGCCgttgt                  >  2:2116/1‑251 (MQ=255)
                                                                                                                                                                                                                       atggcCAGGGCAGGGATAAAGGGCAGTGTGGCAACCAGCAACGCGGTGCGGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCATCTTGGCTCTTTGTACTCACCAGCGAACCGGCAATTGCCATCACCGGGATAAGCGCCCCGAAGCGGCCGACAAACAGGCAGAACGCCAGTAAACAGTTTCAGAACGGAGAGTTGGCGCTTATTCCGGCAAAGGCGCTGCCgttgtt                 <  2:198974/248‑1 (MQ=255)
                                                                                                                                                                                                                             agggcagggATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTGGGCAGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTACTCACCAGCGAACCGGCAATTGCCATCACCGGGATAATCACCCCGAAGCGACCGACAAACATGCAGAACGCCAGTAAACAGTTCCAGAACGGAGAGTTGGCGCTTAATCCGGCAAAGGCGCTGCCGTTGTTGTTAgc           <  2:416669/251‑1 (MQ=255)
                                                                                                                                                                                                                              gtgcagggATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGGCTACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTCCTCACCAGCGAACCGGCAATTGCCATCACCGGGATAATCTCCCCGAAGCGACCGACAAACATTCAGAACGCCAGTAAACAGTTCCAGAACGGAGAGTTGGCGCTTAATCCGGCAAAGGCGCTGCCGTTGTTGTTAgcg          <  1:604345/249‑1 (MQ=255)
                                                                                                                                                                                                                              gggcagggATAAAGGTCAGTGTGTCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTACTCACCAGCGAACCGGCAATTGCCATCACCGGGATAATCACCCCGAAGCGACCGACAAACATGCAGAACGCCAGTAAACAGTTCCAGAACGGAGAGTTGGCGCTTAATCCGGCAAAGGCGCTGCCGTTGTTGTTAgcg          <  1:276869/251‑1 (MQ=255)
                                                                                                                                                                                                                                      atcaaGGTCAGTGTGTCAACCAGCAACACGGTTCCGCTTCACAGGCCCCCAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGAAGCTGGCGGCTTGGCTCTTTTTAATCACCAGCGAACCGGCAATTGCCCTCACCGGGATAATCACCCCGAAGCGACCGACAAACATGCAGAACGCCAGTAAACAGTTCCAGAACGGAGAGTTGGCGCTTAATCCGGCAAAAGCGCTGCCGTTGTTGTTAACGGCGGATgt  >  2:368554/4‑250 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                           
GCTTGCCATCGCGATGCTAATACAGGTGGTCAGCAGACTGCCGATCCAGACGATAAACATCACCGGATTGCGCCATTGCGCCTGCGGGTTTAATTTTTTCACCGCTTCTTTCAGCGCCTGAACGACAAGTGTTGGTTCGAATAGCGCCAGTTGTTTACGACTCATATTCAGTGCTCACTCAATATCATCAGGAGAGATATTCCGCCACCGGACCAAGCGCCAGGGCAGGGATAAAGGTCAGTGCGCCAACCAGCAACACGGTGCCGATTAACAGGCCAACAAACAGCGGGCCGTGCGTTGGCAGCGTGCCGGAGCTGGCGGCTTGGCTCTTTTTACTCACCAGCGAACCGGCAATTGCCATCACCGGGATAATCACCCCGAAGCGACCGACAAACATGCAGAACGCCAGTAAACAGTTCCAGAACGGAGAGTTGGCGCTTAATCCGGCAAAGGCGCTGCCGTTGTTGTTAGCGGCGGATGA  >  NC_000913/726872‑727352

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: