Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,305,424:1 +G 7.8% coding (1175/1986 nt) yjcS ← metallo‑beta‑lactamase superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,305,4221.G 7.8% 178.6 / 8.7 51coding (1177/1986 nt)yjcSmetallo‑beta‑lactamase superfamily protein
Reads supporting (aligned to +/- strand):  ref base . (20/27);  new base G (2/2);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

GTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTA‑ATCATGTCGCCGATTTCATTCAT‑GGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGCGTCGTTACCCCACATATCCAGCGTTTCGTTCAGATACTCGGTCCACTTGCTGGTGTCGCGGGTTTTCGCGCCGCGCAGAGTGTAGAAG  >  NC_000913/4305204‑4305645
                                                                                                                                                                                                                            |                                                                                                                                                                                                                               
gTTGATGTCCCGGGCGTAGCCGGCCAGCGCCTGGACGTAACGTTTACCCATCTCCACCTGACCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAGCCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTA‑ATCATGTCGCCGATTTCATTCAT‑GGTGTAGcc                                                                                                                                                                                                                        <  2:270339/228‑1 (MQ=255)
         acGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGTTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGACAGTGAAGACGGCAGCTAA‑AGTATGCTGCAGATTTCTATTAT‑CGGTTAGACCTGGTACGc                                                                                                                                                                                                               >  1:99136/1‑228 (MQ=255)
         acGGGCAGAACCGCCCAGCGCCTGCACGTAACGGTTACCCATCTCCACCTGACCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCCCGGGCGTTGTGGCTGACAGAACCGTAAAAGCCGCGGCTGGCCCAGTTATTGGCCACGGCAGCGGGCAGCTTA‑ACCAAGACGCCCATTTCATTCAT‑GGTGAAACCCCGGTTccc                                                                                                                                                                                                               >  1:176877/1‑226 (MQ=255)
           gggCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGCCAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCCAGTGCAGGCGGCAGCTGA‑AGCATGTCCCCGTCTTCATTCAT‑GGTGTGGCCCTGGTCAGCCa                                                                                                                                                                                                             >  2:237130/1‑228 (MQ=255)
            ggCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTA‑ATCATGTCGCCGATTTCATTCAT‑GGTGTAGCCCTGGTTCGcc                                                                                                                                                                                                              <  2:176877/226‑1 (MQ=255)
            ggCAGAACCGCCCAGCGCCGGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTA‑ATCATGTCGCCGATTTCATTCAT‑GGTGTAGCCCTGGTTCGcc                                                                                                                                                                                                              <  2:123311/226‑1 (MQ=255)
                      cccAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGAAGGCGGCAGCTTA‑ATCATGTCGCCGATTTCATTCAT‑GGTGTAGCCCTGGTTTGCCCGGGGCAGGGTc                                                                                                                                                                                                  >  2:372614/1‑228 (MQ=255)
                      cccAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGAATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACATCGCGCGGGCGCTGTGGCTGACAGAAACGTAATAGCCGCGCCTAGCCCAGTTCTTGGCAAGTGCAGGCGGCAGCTGACTTCACG‑GGCCGACTCCATCTAA‑GGTGTAGCCCCGGTTCCCCAGGtcacgggta                                                                                                                                                                                                  >  1:261910/1‑220 (MQ=255)
                           cgcCTGCCCGTACCGTTTACCCATCCCCAACTGACCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTA‑ATCATGTCGCCGATTTCATTCAT‑GGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTc                                                                                                                                                                                             <  2:352510/228‑1 (MQ=255)
                                 cACGTAACGTTTACCCATCTCCACCTGACCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCCGTTATTGGCAAGTGCAGGCGGCAGCTTA‑ATCAGGCCGCCCATTTCATTCTTGGGGGTTGCCCTGGTTCCCCAGGGGCAgggcccgggcggcaa                                                                                                                                                                                        >  1:71374/1‑216 (MQ=255)
                                 cACGTAACGTTTACCCATCTCCACCTGACCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCCCGTCCGGCCCAGGTATTGGCAAGTGCAGGCGGCAGCTTA‑ATCACGCGCCCGATTTCTTTCATGGGTGCAGCCCGGGTGCGCCAGGGGGCGGGCCTGTTCGTGaa                                                                                                                                                                                        >  1:644903/1‑228 (MQ=255)
                                              cccATCTCCACCTGACCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTA‑ATCATGTCGCCGCTTTCCTTCAC‑GGTGTATGCCCGGGTCGCCAGGTGCAGGGTCTGGTCGTAGATGAAATTGAAGGTa                                                                                                                                                                          >  1:56862/1‑228 (MQ=255)
                                                 aTCTCCATCTGCCCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTA‑ATCATGTCGCCGATTTCATTCAT‑GGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATcg                                                                                                                                                                       <  2:239094/228‑1 (MQ=255)
                                                  tctcCACCTGACCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGACGGCCGCTTA‑ATCATGTCGCCGATTTCATTCAT‑GGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCCTGaaatgagttttaggaattgg                                                                                                                                                                      >  1:351582/1‑210 (MQ=255)
                                                           ggccATGCGGATGCAGGTTGGCCGGGTTACCGTCAGAATTGCCAAGATAGAAGTTATCCCCCGCGCGGGCGTTGTTGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTA‑ATCATGTCGCCGATTTCATTCAT‑GGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTAc                                                                                                                                                             <  2:99136/226‑1 (MQ=255)
                                                           gACCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACGCCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTA‑ATCATGTCGCCGATTTCATTCAT‑GGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTAc                                                                                                                                                             <  2:62812/228‑1 (MQ=255)
                                                                tACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGGCGCAGCTTA‑ATCCTGTCGCCGATTTCATTCAT‑GGTGTAGCCCCGGTTCGCCAGGCGCACGGCCTTGGCGTGCATGTATTTGATgggtttgttcccttttaaaaat                                                                                                                                                        >  1:90288/1‑208 (MQ=255)
                                                                 aCGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTA‑ATCATGTCGCCGATTTCATTCAT‑GGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTa                                                                                                                                                                                    <  2:90288/199‑1 (MQ=255)
                                                                  cGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGCTAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTA‑ATCATGTCGCCGATTTCATTCATGGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAatat                                                                                                                                                       <  2:644903/228‑1 (MQ=255)
                                                                  cGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTA‑ATCATGTCGCCGATTTCATTCATGGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAatat                                                                                                                                                       <  2:71374/228‑1 (MQ=255)
                                                                              cgCCGGGTTTCCGTGCTAATAGCCAAGATAGAAGTTATTCAGCGCGCGGGGGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTA‑ATCATGTCGCCGATTTCATTCAT‑GGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATg                                                                                                                                          <  1:683085/227‑1 (MQ=255)
                                                                                  gggTTCCCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTA‑ACCATGTCGCCGATTTCATTCAT‑GGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTa                                                                                                                                      <  2:649413/228‑1 (MQ=255)
                                                                                  gggTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGTGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTA‑ATCATGTCGCCGATTTCATTCAT‑GGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTa                                                                                                                                      <  2:261910/228‑1 (MQ=255)
                                                                                   ggTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGGTATTGGCACGTGCAGGCGGCAGATTA‑ATCATGGCGCCGATTTTATTGAT‑GGTGTAGCCCCGGTTCGCCAGGGGCAGGGCCTGGTGGGGAATGTACTTGGTGGGAATGCGGGATTTACCAAAATAAACAATGAAAAGCTTAt                                                                                                                                     >  1:65713/1‑228 (MQ=255)
                                                                                     ttACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTA‑ATCATGTCGCCGATTTCATTCAT‑GGGGGAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGGGAATGTACTTGATGGGATTCCGGAATTTCACAATATAATAATTGTTAATCTTAttt                                                                                                                                   >  1:415367/1‑227 (MQ=255)
                                                                                      tACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTA‑ATCATGTCGCCGATTTCATTCAT‑GGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGatat                                                                                                                                           <  2:415367/219‑1 (MQ=255)
                                                                                            gatAATAGCCAAGATAGAAGTTACCCCCCGCGGGGGGGTTGTGGCTGACAGACCCGGAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTA‑ATCATGTCGCCGATTGCATTCAT‑GGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGa                                                                                                                            <  1:372614/226‑1 (MQ=255)
                                                                                            cgtAAGTGCCAAGATAGAAGTTATGCCCCGGGTGGGCGTTGTGGCGGACAGAACCGTAATGGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTA‑AGCATGTCGCCGATTTCATTCAT‑GGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGa                                                                                                                            <  1:202023/228‑1 (MQ=255)
                                                                                              tAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTA‑ATCATGTCGCCGATTTCATTCAT‑GGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGGATCGCGGGATTTACCAATATAATCATTGATATGCTTATTGCCCCAGAcc                                                                                                                          >  2:485239/1‑228 (MQ=255)
                                                                                               aaTAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTA‑ATCATGTCGCCGATTTCATTCAT‑GGTGTAGACCTGGTTCGCCAGGGGCAAGGTCTGGTCGTTAATGTACTTGGTGGTATCGCGGATTTTACCAATAAAATCATTGATATGCTTATTGCCCACGACCg                                                                                                                         >  1:76424/1‑228 (MQ=255)
                                                                                                     cAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTA‑ATCATGTCGCCGATTTCATTCAT‑GGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATTTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAgg                                                                                                                   >  1:444308/1‑228 (MQ=255)
                                                                                                        gattggagttgtcctccggcagggTGTTGGGGCTTACAGAACCTTAATACCCGAGGTTGGCCAAGTTATTGGAAAGTGCAGGCGGCCGCTTA‑ATCATGTCGCCGATTTCATTCAT‑GGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGgtgt                                                                                                                <  1:575559/207‑1 (MQ=255)
                                                                                                        gATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTA‑ATCATGTCGCCGATTTCATTCAT‑GGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGgtgt                                                                                                                <  2:351582/228‑1 (MQ=255)
                                                                                                        gATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTA‑ATCATGTCGCCGATTTCATTCAT‑GGTGTAGCCCTGGTTCGCCAGGTGCAGGGGCTGGTCGTGAATGTACTTGATGGTATCCCGGCATTTACCAATATAATCATTGATATGCTTATTGCCCAAGCCCGGctcggcgg                                                                                                                >  2:647949/1‑221 (MQ=255)
                                                                                                                ttATACCCCGCGCGGGGGTTGTGGCTGGCAGACCCGGAATAGCCGCGGCTGGCCCATTTAGTGGCAAGTGCGGGCGGCAGCTTA‑ATCATGTCGCCGATTTCATTCAT‑GGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATa                                                                                                        <  1:483761/228‑1 (MQ=255)
                                                                                                                   taCACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTA‑ATCATGTCGCCGATTCCATTCAT‑GGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGGATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATaaa                                                                                                      >  2:81938/1‑227 (MQ=255)
                                                                                                                         gcgcgGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTA‑ATCATGTCGCCGATTTCATTCAT‑GGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAata                                                                                                                                                        <  2:371071/171‑1 (MQ=255)
                                                                                                                          cgcgGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTA‑ATCATGTCGCCGATTTCATTCAT‑GGTGTAGCCCTGGTTCGCCAGGTGCAGGGCCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAGTCATTGATATGCTAATTGCCCCAGACAGGCCAGGAGTGCGCAATAAacaacaat                                                                                               >  2:170771/1‑225 (MQ=255)
                                                                                                                                  ttGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTA‑ATCATGTCGCCGATTTCATTCAT‑GGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGGATCGCGGTATTTACCAATATAATCCTTGATATGCTTATTGCCCCAGACCGGCCAGGGGTGGGGCATAAACGACACTTCCGcgtcg                                                                                      >  2:256122/1‑228 (MQ=255)
                                                                                                                                   tggggcgtgcAGAACCGTAATTGCCCGGGGTGGCCCATGTATTGGCAAGTGCAGGCGGCAGCTTA‑ATCATGGCGCCGATTCCATTCAT‑GGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGcgtcgt                                                                                     <  1:81938/219‑1 (MQ=255)
                                                                                                                                           aCAGAACCGTAATAGCCGCGGCTTGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTA‑ATCATGTCGCCGATTTCAGTCAT‑GGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGCGTCGTTACCCca                                                                              <  2:65713/227‑1 (MQ=255)
                                                                                                                                                          ccGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTA‑ATCATGTCGCCGATTTCATTCAT‑GGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGGATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGCGTCGTTACCCCCCATATCCAGCGTTTCg                                                              >  2:32961/1‑228 (MQ=255)
                                                                                                                                                                        tttATTGGGATGTCGAGTCGGCAGTTGA‑ATCATGTCGGCGATTTCATTAAT‑GGTGTAGCCCGGGTTCGCCAGGTTCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGCGTCGTTCCCCCACATATCCAGCGTTTCGTTCAGAttagaggc                                                <  1:237130/227‑8 (MQ=255)
                                                                                                                                                                                   gTGCAGGCGGCAGCTTA‑ATCATGTCGCCGATTTCATTCAT‑GGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGAATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTTCGGCATAAACAGCACTTCCGCGTCTTTGCCCCACATATCCAGCGTTTCGTTCAGATACCCGTGCACCTTGCTgg                                     >  2:266010/1‑228 (MQ=255)
                                                                                                                                                                                              aGCTTA‑ATCATGTCGCCGCTTTCATTCAT‑GGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACCGCCCTTCCGCGTCGTTACCCCCCATATTCAGCGTTTCGTTCAGATACTCGGTCCACTTGCTGGTGTCGCGGGt                           >  2:407780/1‑227 (MQ=255)
                                                                                                                                                                                                      tCATGTCGGCGATTTCATTCAT‑GGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCCCTTCCGCGTCGTTACCCCACATATCCAGCGTTTCGTTCAGATACTCGGTCCACTTGCTGGTGTCGCGGGTTTTcgcg                    <  2:85044/227‑1 (MQ=255)
                                                                                                                                                                                                          gacGCCGATTTCATTCAT‑GGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGCGTCGTTACCCCACATATCCAGCGTTTCGTTCAGATACTCGGTCCACTTGCTGGTGTCGCGGGTTTTcgcgccgcg               <  2:524371/226‑1 (MQ=255)
                                                                                                                                                                                                           ttgCCGTTTTCGTTCAT‑GGTGTAGCCGTGGGTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTAGTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGCGTCGTTACCCCACATATCCAGCGTTTCGTTCAGATACTCGGTCCACTTGCTGGTGTCGCGGGTTTTcgcgccgcg               <  2:338783/225‑1 (MQ=255)
                                                                                                                                                                                                           cggctgttttCATTCAT‑GGTGGAGCTCGTGTTCGCTAAGTGCAGGGTCTTGTCGTGAATGTACTTGATGGTATCGCTGTATTTACCACGATAAGTATTGATATGCTGAGTACCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGCGTCGTCGCCCCACATGTCCAGCGTTGCGTTCAGATTCTCGGTCCACTTGCTGGTGTCGCGGGTTTTcgcgccgcgc              <  1:220772/221‑1 (MQ=255)
                                                                                                                                                                                                            cgcggttttCATTGAT‑GGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGCGTTGAATGTACTTGAGGGTCTCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGTCCAGGCGTGCGGCATAAACAGCACTTCCGCGGCGTTACCCCACATATCCAGCGTTTCGTTCAGATACTCGGTCCACTTGCTGGTGTCGCGGGTTTTCGCGCCGCGCa             <  1:210118/222‑1 (MQ=255)
                                                                                                                                                                                                             gCCGATTTCATTCAT‑GGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGCGTCGTTACCCCACATATCCAGCGTTTCGTTCAGATACTCGGTCCACTTGCTGGTGTCGCGGGTTTTCGCGCCGCGCa             <  2:76424/227‑1 (MQ=255)
                                                                                                                                                                                                                        tcat‑GGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGCGTCGTTACCCCACATATCCAGCGTTTCGTTCAGATACTCGGTCCACTTGCTGGTGTCGCGGGTTTTCGCGCCGCGCAGAGTGTAGAAg  <  2:444308/227‑1 (MQ=255)
                                                                                                                                                                                                                            |                                                                                                                                                                                                                               
GTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTA‑ATCATGTCGCCGATTTCATTCAT‑GGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGCGTCGTTACCCCACATATCCAGCGTTTCGTTCAGATACTCGGTCCACTTGCTGGTGTCGCGGGTTTTCGCGCCGCGCAGAGTGTAGAAG  >  NC_000913/4305204‑4305645

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