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

GTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGCGTCGTTACCCCACATATCCAGCGTTTCGTTCAGATACTCGGTCCACTTGCTGGTGTCGCGGGTTTTCGCGCCGCGCAGAGTGTAGAAG  >  NC_000913/4305204‑4305645
                                                                                                                                                                                                                            |                                                                                                                                                                                                                               
gTTGATGTCCCGGGCGTAGCCGGCCAGCGCCTGGACGTAACGTTTACCCATCTCCACCTGACCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAGCCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGcc                                                                                                                                                                                                                        <  2:270339/228‑1 (MQ=255)
         acGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGTTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGACAGTGAAGACGGCAGCTAAAGTATGCTGCAGATTTCTATTATCGGTTAGACCTGGTACGc                                                                                                                                                                                                               >  1:99136/1‑228 (MQ=255)
         acGGGCAGAACCGCCCAGCGCCTGCACGTAACGGTTACCCATCTCCACCTGACCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCCCGGGCGTTGTGGCTGACAGAACCGTAAAAGCCGCGGCTGGCCCAGTTATTGGCCACGGCAGCGGGCAGCTTAACCAAGACGCCCATTTCATTCATGGTGAAACCCCGGTTccc                                                                                                                                                                                                               >  1:176877/1‑226 (MQ=255)
           gggCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGCCAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCCAGTGCAGGCGGCAGCTGAAGCATGTCCCCGTCTTCATTCATGGTGTGGCCCTGGTCAGCCa                                                                                                                                                                                                             >  2:237130/1‑228 (MQ=255)
            ggCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGcc                                                                                                                                                                                                              <  2:176877/226‑1 (MQ=255)
            ggCAGAACCGCCCAGCGCCGGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGcc                                                                                                                                                                                                              <  2:123311/226‑1 (MQ=255)
                      cccAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGAAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTTGCCCGGGGCAGGGTc                                                                                                                                                                                                  >  2:372614/1‑228 (MQ=255)
                      cccAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGAATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACATCGCGCGGGCGCTGTGGCTGACAGAAACGTAATAGCCGCGCCTAGCCCAGTTCTTGGCAAGTGCAGGCGGCAGCTGACTTCACGGGCCGACTCCATCTAAGGTGTAGCCCCGGTTCCCCAGGtcacgggta                                                                                                                                                                                                  >  1:261910/1‑220 (MQ=255)
                           cgcCTGCCCGTACCGTTTACCCATCCCCAACTGACCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTc                                                                                                                                                                                             <  2:352510/228‑1 (MQ=255)
                                 cACGTAACGTTTACCCATCTCCACCTGACCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCCGTTATTGGCAAGTGCAGGCGGCAGCTTAATCAGGCCGCCCATTTCATTCTTGGGGGTTGCCCTGGTTCCCCAGGGGCAgggcccgggcggcaa                                                                                                                                                                                        >  1:71374/1‑216 (MQ=255)
                                 cACGTAACGTTTACCCATCTCCACCTGACCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCCCGTCCGGCCCAGGTATTGGCAAGTGCAGGCGGCAGCTTAATCACGCGCCCGATTTCTTTCATGGGTGCAGCCCGGGTGCGCCAGGGGGCGGGCCTGTTCGTGaa                                                                                                                                                                                        >  1:644903/1‑228 (MQ=255)
                                              cccATCTCCACCTGACCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGCTTTCCTTCACGGTGTATGCCCGGGTCGCCAGGTGCAGGGTCTGGTCGTAGATGAAATTGAAGGTa                                                                                                                                                                          >  1:56862/1‑228 (MQ=255)
                                                 aTCTCCATCTGCCCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATcg                                                                                                                                                                       <  2:239094/228‑1 (MQ=255)
                                                  tctcCACCTGACCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGACGGCCGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCCTGaaatgagttttaggaattgg                                                                                                                                                                      >  1:351582/1‑210 (MQ=255)
                                                           ggccATGCGGATGCAGGTTGGCCGGGTTACCGTCAGAATTGCCAAGATAGAAGTTATCCCCCGCGCGGGCGTTGTTGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTAc                                                                                                                                                             <  2:99136/226‑1 (MQ=255)
                                                           gACCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACGCCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTAc                                                                                                                                                             <  2:62812/228‑1 (MQ=255)
                                                                tACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGGCGCAGCTTAATCCTGTCGCCGATTTCATTCATGGTGTAGCCCCGGTTCGCCAGGCGCACGGCCTTGGCGTGCATGTATTTGATgggtttgttcccttttaaaaat                                                                                                                                                        >  1:90288/1‑208 (MQ=255)
                                                                 aCGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTa                                                                                                                                                                                    <  2:90288/199‑1 (MQ=255)
                                                                  cGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGCTAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAatat                                                                                                                                                       <  2:644903/228‑1 (MQ=255)
                                                                  cGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAatat                                                                                                                                                       <  2:71374/228‑1 (MQ=255)
                                                                              cgCCGGGTTTCCGTGCTAATAGCCAAGATAGAAGTTATTCAGCGCGCGGGGGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATg                                                                                                                                          <  1:683085/227‑1 (MQ=255)
                                                                                  gggTTCCCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAACCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTa                                                                                                                                      <  2:649413/228‑1 (MQ=255)
                                                                                  gggTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGTGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTa                                                                                                                                      <  2:261910/228‑1 (MQ=255)
                                                                                   ggTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGGTATTGGCACGTGCAGGCGGCAGATTAATCATGGCGCCGATTTTATTGATGGTGTAGCCCCGGTTCGCCAGGGGCAGGGCCTGGTGGGGAATGTACTTGGTGGGAATGCGGGATTTACCAAAATAAACAATGAAAAGCTTAt                                                                                                                                     >  1:65713/1‑228 (MQ=255)
                                                                                     ttACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGGGGAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGGGAATGTACTTGATGGGATTCCGGAATTTCACAATATAATAATTGTTAATCTTAttt                                                                                                                                   >  1:415367/1‑227 (MQ=255)
                                                                                      tACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGatat                                                                                                                                           <  2:415367/219‑1 (MQ=255)
                                                                                            gatAATAGCCAAGATAGAAGTTACCCCCCGCGGGGGGGTTGTGGCTGACAGACCCGGAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTGCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGa                                                                                                                            <  1:372614/226‑1 (MQ=255)
                                                                                            cgtAAGTGCCAAGATAGAAGTTATGCCCCGGGTGGGCGTTGTGGCGGACAGAACCGTAATGGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAAGCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGa                                                                                                                            <  1:202023/228‑1 (MQ=255)
                                                                                              tAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGGATCGCGGGATTTACCAATATAATCATTGATATGCTTATTGCCCCAGAcc                                                                                                                          >  2:485239/1‑228 (MQ=255)
                                                                                               aaTAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGACCTGGTTCGCCAGGGGCAAGGTCTGGTCGTTAATGTACTTGGTGGTATCGCGGATTTTACCAATAAAATCATTGATATGCTTATTGCCCACGACCg                                                                                                                         >  1:76424/1‑228 (MQ=255)
                                                                                                     cAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATTTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAgg                                                                                                                   >  1:444308/1‑228 (MQ=255)
                                                                                                        gattggagttgtcctccggcagggTGTTGGGGCTTACAGAACCTTAATACCCGAGGTTGGCCAAGTTATTGGAAAGTGCAGGCGGCCGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGgtgt                                                                                                                <  1:575559/207‑1 (MQ=255)
                                                                                                        gATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGgtgt                                                                                                                <  2:351582/228‑1 (MQ=255)
                                                                                                        gATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGGCTGGTCGTGAATGTACTTGATGGTATCCCGGCATTTACCAATATAATCATTGATATGCTTATTGCCCAAGCCCGGctcggcgg                                                                                                                >  2:647949/1‑221 (MQ=255)
                                                                                                                ttATACCCCGCGCGGGGGTTGTGGCTGGCAGACCCGGAATAGCCGCGGCTGGCCCATTTAGTGGCAAGTGCGGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATa                                                                                                        <  1:483761/228‑1 (MQ=255)
                                                                                                                   taCACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTCCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGGATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATaaa                                                                                                      >  2:81938/1‑227 (MQ=255)
                                                                                                                         gcgcgGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAata                                                                                                                                                        <  2:371071/171‑1 (MQ=255)
                                                                                                                          cgcgGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGCCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAGTCATTGATATGCTAATTGCCCCAGACAGGCCAGGAGTGCGCAATAAacaacaat                                                                                               >  2:170771/1‑225 (MQ=255)
                                                                                                                                  ttGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGGATCGCGGTATTTACCAATATAATCCTTGATATGCTTATTGCCCCAGACCGGCCAGGGGTGGGGCATAAACGACACTTCCGcgtcg                                                                                      >  2:256122/1‑228 (MQ=255)
                                                                                                                                   tggggcgtgcAGAACCGTAATTGCCCGGGGTGGCCCATGTATTGGCAAGTGCAGGCGGCAGCTTAATCATGGCGCCGATTCCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGcgtcgt                                                                                     <  1:81938/219‑1 (MQ=255)
                                                                                                                                           aCAGAACCGTAATAGCCGCGGCTTGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCAGTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGCGTCGTTACCCca                                                                              <  2:65713/227‑1 (MQ=255)
                                                                                                                                                          ccGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGGATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGCGTCGTTACCCCCCATATCCAGCGTTTCg                                                              >  2:32961/1‑228 (MQ=255)
                                                                                                                                                                        tttATTGGGATGTCGAGTCGGCAGTTGAATCATGTCGGCGATTTCATTAATGGTGTAGCCCGGGTTCGCCAGGTTCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGCGTCGTTCCCCCACATATCCAGCGTTTCGTTCAGAttagaggc                                                <  1:237130/227‑8 (MQ=255)
                                                                                                                                                                                   gTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGAATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTTCGGCATAAACAGCACTTCCGCGTCTTTGCCCCACATATCCAGCGTTTCGTTCAGATACCCGTGCACCTTGCTgg                                     >  2:266010/1‑228 (MQ=255)
                                                                                                                                                                                              aGCTTAATCATGTCGCCGCTTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACCGCCCTTCCGCGTCGTTACCCCCCATATTCAGCGTTTCGTTCAGATACTCGGTCCACTTGCTGGTGTCGCGGGt                           >  2:407780/1‑227 (MQ=255)
                                                                                                                                                                                                      tCATGTCGGCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCCCTTCCGCGTCGTTACCCCACATATCCAGCGTTTCGTTCAGATACTCGGTCCACTTGCTGGTGTCGCGGGTTTTcgcg                    <  2:85044/227‑1 (MQ=255)
                                                                                                                                                                                                          gacGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGCGTCGTTACCCCACATATCCAGCGTTTCGTTCAGATACTCGGTCCACTTGCTGGTGTCGCGGGTTTTcgcgccgcg               <  2:524371/226‑1 (MQ=255)
                                                                                                                                                                                                           ttgCCGTTTTCGTTCATGGTGTAGCCGTGGGTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTAGTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGCGTCGTTACCCCACATATCCAGCGTTTCGTTCAGATACTCGGTCCACTTGCTGGTGTCGCGGGTTTTcgcgccgcg               <  2:338783/225‑1 (MQ=255)
                                                                                                                                                                                                           cggctgttttCATTCATGGTGGAGCTCGTGTTCGCTAAGTGCAGGGTCTTGTCGTGAATGTACTTGATGGTATCGCTGTATTTACCACGATAAGTATTGATATGCTGAGTACCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGCGTCGTCGCCCCACATGTCCAGCGTTGCGTTCAGATTCTCGGTCCACTTGCTGGTGTCGCGGGTTTTcgcgccgcgc              <  1:220772/221‑1 (MQ=255)
                                                                                                                                                                                                            cgcggttttCATTGATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGCGTTGAATGTACTTGAGGGTCTCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGTCCAGGCGTGCGGCATAAACAGCACTTCCGCGGCGTTACCCCACATATCCAGCGTTTCGTTCAGATACTCGGTCCACTTGCTGGTGTCGCGGGTTTTCGCGCCGCGCa             <  1:210118/222‑1 (MQ=255)
                                                                                                                                                                                                             gCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGCGTCGTTACCCCACATATCCAGCGTTTCGTTCAGATACTCGGTCCACTTGCTGGTGTCGCGGGTTTTCGCGCCGCGCa             <  2:76424/227‑1 (MQ=255)
                                                                                                                                                                                                                        tcatGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGCGTCGTTACCCCACATATCCAGCGTTTCGTTCAGATACTCGGTCCACTTGCTGGTGTCGCGGGTTTTCGCGCCGCGCAGAGTGTAGAAg  <  2:444308/227‑1 (MQ=255)
                                                                                                                                                                                                                            |                                                                                                                                                                                                                               
GTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGCGTCGTTACCCCACATATCCAGCGTTTCGTTCAGATACTCGGTCCACTTGCTGGTGTCGCGGGTTTTCGCGCCGCGCAGAGTGTAGAAG  >  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: