Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,304,821 T→C 6.6% N593D (AAC→GAC)  yjcS ← metallo‑beta‑lactamase superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,304,8210TC 6.6% 177.1 / 4.8 61N593D (AAC→GAC) yjcSmetallo‑beta‑lactamase superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (38/19);  new base C (2/2);  total (40/21)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-02

CAAAAGCGAGTGACATTAATTTGGGGTTACGATATTCACCCACAAATCGAAATTATCCAGACAGGCGATAATTTCTTCCAGTTTCGCGCCATTGCCAATAATTTTGGCTTTCTTCGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGT  >  NC_000913/4304597‑4305039
                                                                                                                                                                                                                                |                                                                                                                                                                                                                           
cAAAGGCGAGTGACATTAATTTGGGGTTACGATATTCACCCACAAATCGAAATTATCCAGACGGGCGATAATTTCTTCCAGTTTCGCGCCATTGCCAATAATTTTGGCTTTCTTCGTTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAg                                                                                                                                                                                                                          <  1:321295/228‑1 (MQ=255)
       gAGTGACATTAATTTGGGGTTACGATATTCACCCACAAATCGAAATTATCCAGACAGGCGATAATTTCTTCCAGTTTCGCGCCATTGCCAATAATTTTGGCTTTCTTCGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGCTCAGCACGCTCTGCCGATCTTCCCGGCCGTTGTAGACATAGGCGTCCGCTTGCGGTTGCGCCGTTTTCCGGTAGTTAATCCCGCTa                                                                                                                                                                                                                   >  1:67242/1‑228 (MQ=255)
         gTGACATTAATTTGGGGTTACGATATTCACCCACAAATCGAAATTATCCAGACAGGCGATAATTTCTTCCAGTTTCGCGCCATTGCCAATAATTTTGGCTTTCTTCGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCGGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGGTGCAGCGGTTTCCGGTGGTTAAGCACGCTATc                                                                                                                                                                                                                 >  1:525766/1‑228 (MQ=255)
                       gggTTACGATATTCACCCACAAATCGAAATTATCCAGACAGGCGATAATTTCTTCCAGTTTCGCGCCATTGCCAATAATTTTGGCTTTCTTCGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAATAGGCGTCGGCTTGCGGTTGAAGCGTTTTCCGGTAGTTAAGCACGCCATCGTTCAGCGTCCgg                                                                                                                                                                                                    >  1:100436/1‑227 (MQ=255)
                        gtttattattttCCCCCACAATGCAAAATTTTCCGTACGGGGGATAGTTTTTTCCGTGTTCGCGCCATTGCCATTATTTTTGGTTTTCTCCGCTTTTACCAGTCCCGCCCTTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGtt                                                                                                                                                                                                  <  2:263968/218‑1 (MQ=255)
                        ggTTACGATATTCACCCACAAATCGAAATTATCCAGACAGGCGATAATTTCTTCCAGTTTCGCGCCATTGCCAATAATTTTGGCTTTCTTCGCTTTTACCAGATCCGCCATTTTGGCTTTTCCGGTCAGCCCGGCGTGCAGATCTTCACGGGTGTTGGAGAAAGAGACGCCGGCTTGCGGGTGCAGCGTTTTCCCGTAGTTAAGAACGCCACCGTTCAGCGTCCgggg                                                                                                                                                                                                  >  2:459856/1‑226 (MQ=255)
                            aCGATATTCACCCACAAATCGAAATTATCCAGACAGGCGATAATTTCTTCCAGTTTCGCGCCATTGCCAATAATTTTGGCTTTCTTCGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGCgg                                                                                                                                                                                              >  1:353044/1‑228 (MQ=255)
                            aCGATATTCACCCACAAATCGAAATTATCCAGACAGGCGATAATTTCTTCCAGTTTCGCGCCATTGCCAATAATTTTGGCTTTCTTCGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAgg                                                                                                                                                                                              >  1:263968/1‑228 (MQ=255)
                             cGATATTCACCCACAAATCGAAATTATCCAGACAGGCGATAATTTCTTCCAGTTTCGCGCCATTGCCAATAATTTTGGCTTTCTTCGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGGAGCGTTTTCTGGTAGTTAAGCATGCTATCGTTCAGCGTCTGGGTTAgg                                                                                                                                                                                              >  1:270142/1‑227 (MQ=255)
                                       ccACAAATCGAAATTATCCAGACAGGCGATAATTTCTTCCAGTTTCGCGCCATTGCCAATAATTTTGGCTTTCTTCGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAGGGTATTGCCGtt                                                                                                                                                                                  >  1:205608/1‑228 (MQ=255)
                                          caAATCGAAATTATCCAGACAGGCGATAATTTCTTCCAGTTTCGCGCCATTGCCAATAATTTTGGCTTTCTTCGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGTATTAGAAAGAGGCGTCGGCTTTCGGTTGCACCGGTTTCCGGTAGTTAAGACCGCTATCTGTCAACGCAAGGTTGGGGGTTTCCCCTTTggc                                                                                                                                                                               >  2:388690/1‑226 (MQ=255)
                                          caAATCGAAATTATCCAGACAGGCGATAATTTCTTCCAGTTTCGCGCCATTGCCAATAATTTTGGCTTTCTTCGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACg                                                                                                                                                                                                                      >  2:409202/1‑190 (MQ=255)
                                          caAATCGAAATTATCCAGACAGGCGATAATTTCTTCCAGTTTCGCGCCATTGCCAATAATTTTGGCTTTCTTCGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACg                                                                                                                                                                                                                      <  1:409202/190‑1 (MQ=255)
                                          caAATCGAAATTATCCAGACAGGCGATAATTTCTTCCAGTTTCGCGCCATTGCCAATAATTTTGGCTTTCTTCGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGCCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAACACGCTATCGTTCAGCGGCAGGGTGGGGTTGTCGcccattct                                                                                                                                                                               >  2:205515/1‑222 (MQ=255)
                                            aaTCGAAATTATCCAGACAGGCGATAATTTCTTCCAGTTTCGCGCCATTGCCAATAATTTTGGCTTTCTTCGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCCGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCCGGGTGAGGTTATCGCCGTTGCTcc                                                                                                                                                                             >  1:379123/1‑227 (MQ=255)
                                            aaTCAAATTTTTCCAGATCGGCGATTATTTTTTCCAGTTTCGCGCCATTGCCAATAATTGTGGCTTTCTCCGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCa                                                                                                                                                                             <  2:387730/228‑1 (MQ=255)
                                                 aaaTTATCCAGACAGGCGATAATTTCTTCCAGTTTCGCGCCATTGCCAATAATTTTGGCTTTCTTCGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCCCGGCTGATGTAGAAAGAGGCGTCGGCGTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCTGGTTGGGGGTGTGCCCGTTGCTCagt                                                                                                                                                                           >  1:93324/1‑223 (MQ=255)
                                                       tCCAGACAGGCGATAATTTCTTCCAGTTTCGCGCCATTGCCAATAATTTTGGCTTTCTTCGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCCGATACTCACGGCTGGTGTAGAAAGATGCGTCGGCGTGCCGTTGCAGCGTTTTCCGGTAGGTAAGGCCGCTAGCGTTCAGCGTCAGGGTGAGGTTTACCCCCGTTGCCATTTTGAAAttg                                                                                                                                                                  >  2:436208/1‑227 (MQ=255)
                                                            aCAGGCGATAATTTCTTCCAGTTTCGCGCCATTGCCAATAATTTTGGCTTTCTTCGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGc                                                                                                                                                              >  1:197148/1‑227 (MQ=255)
                                                                gCGATAATTTCTTCCAGTTTCGCGCCATTGCCAATAATTTTGGCTTTCTTCGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGata                                                                                                                                                         >  1:196992/1‑228 (MQ=255)
                                                                gCGATAATTTCTTCCAGTTTCGCGCCATTGCCAATAATTTTGGCTTTCTTCGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTTAGGTTATCGCCGTTGGTCATATTGCAGTTCCGGCTGata                                                                                                                                                         >  1:542972/1‑228 (MQ=255)
                                                                    tttttttttCCGGTTTCCCGCCATTGGCAATTATTTTGGCTTTCTTCGCTTTGACCAGGTCCGCCATTTTGGGTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATAtttt                                                                                                                                                     <  2:205608/225‑1 (MQ=255)
                                                                          cttcCAGTTTCGCGCCATTGCCAATAATTTTGGCTTTCTTCGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGCAGTTAAGCACGCAATCGTTCACCGTCAGGTTGTGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTAccg                                                                                                                                                >  1:244198/1‑226 (MQ=255)
                                                                                     gcgcCATTGCCAATAATTTTGGCTTTCTTCGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAGGGTATCGCCGTTGGTCATATTGAAGTTCAGGCTGATATTTTTAACCGCAGCTTTcggg                                                                                                                                    >  1:278496/1‑226 (MQ=255)
                                                                                                   tgttttGTCTTTCTTCGTTTTTACATGGTCCGACATTTGGGCTGGCCCGGGCGGCGCGGGGTGCAGATGTTTACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGGAGCGTTTTCCGGTAGTTAAGCACGTTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAAc                                                                                                                      <  2:222731/226‑1 (MQ=255)
                                                                                                       ttGGCTTTCTTCGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTTAAGTTCAGGCTGATATTTTTACCCGCACCTTTTGCGCTATCGAGGCGACCgcc                                                                                                                   >  1:222731/1‑225 (MQ=255)
                                                                                                        tGGCTTTCTTCGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGCAGTTCCGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGccc                                                                                                                 >  1:29838/1‑227 (MQ=255)
                                                                                                                        tACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGGTGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGGCAGGTGGGGGTTATCGCCGTTCCTCAAATTTAAGTTCAGGCTGATATTTTTACCCGCGCCTTTCGCGCTATCGGGGGGACCGGCCATAAAGTGGacagcca                                                                                                  >  2:438159/1‑221 (MQ=255)
                                                                                                                                  gCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGAAATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGCACAGCATTTcgcccgac                                                                                       >  1:588582/1‑228 (MQ=255)
                                                                                                                                   ccATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGAGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGTACGGCTCTAATGGCGAAAAGCATTTCGACCGCCa                                                                                      >  1:102346/1‑228 (MQ=255)
                                                                                                                                    taTTTTGGCTTGTCCGGTCCGCACGGCGTGCAGGTCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACAt                                                                                     <  1:41004/227‑1 (MQ=255)
                                                                                                                                            ctttgcgggTCAGCACTGCGTTCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGaa                                                                             <  1:618634/221‑1 (MQ=255)
                                                                                                                                                tCCGGTCAGCACGGCGTGCAGCTCTTCAAGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGGTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCATACAGCATTTCGACCACCATCCAGCGAAAgg                                                                          >  1:160781/1‑227 (MQ=255)
                                                                                                                                                tCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGGAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTCAAGCAc                                                                                                                                                                                                                       >  2:268981/1‑87 (MQ=255)
                                                                                                                                                 ccGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTCAAGCACg                                                                                                                                                                                                                      <  1:268981/87‑1 (MQ=255)
                                                                                                                                                     tCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACACCATTTCGACCGCCATCCCGCGAATGGTGGCCg                                                                    >  1:482360/1‑228 (MQ=255)
                                                                                                                                                      tcgcccGGCGTGCAGGTTTTCACGGGTTATGTAAAAAGGGGCGTCGGCTTGGGGTTGCAGGGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTCCgg                                                                   <  2:379123/223‑1 (MQ=255)
                                                                                                                                                              cGTGCAGATCTTCCAGGCTGATGTAGAAAGGGGCGTCGGCTTGCGGTTGCAGCTTTTTCCGGTAGTTAAGCCCGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTCCGGGGAACCgg                                                           <  1:545358/228‑1 (MQ=255)
                                                                                                                                                                     atcgtccgtgCTGGTGGAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGcc                                                    <  1:210292/219‑1 (MQ=255)
                                                                                                                                                                            tggttGATGTAGAAAGGGGCGTCGGTTTGGGGTTGGAGCGTTTTCCGGTAGTTAAGCACGTTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTg                                             <  2:542972/224‑1 (MQ=255)
                                                                                                                                                                            gggCTGATGTAGAAGGGGGGGTCGGTTTGCGGTTGGAGCGTTTTCCGGTAGTTAAGCACGCTATCGTGCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTg                                             <  2:197148/227‑1 (MQ=255)
                                                                                                                                                                               cTGATGGAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCCTAAAGTCGAACAGCATTTCCCCCGACATCCCGCGAATGGGGTCCGGGCAAACCGGGGAGCCGTGGCGGaaa                                          >  2:210292/1‑227 (MQ=255)
                                                                                                                                                                                  aTGTAGGAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGGTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTGGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTTGCCCCACATCCCCGGAATGGTGTCCGGGGCACCGTGGGGGCGCGGGCGGAACTTa                                       >  2:542358/1‑228 (MQ=255)
                                                                                                                                                                                  aTGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTCCGGGGAAc                                                              <  1:367497/205‑1 (MQ=255)
                                                                                                                                                                                  aTGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGGGTCCGGGGAACCGTGGGTGCCTTGGCTGAACTTa                                       >  1:390211/1‑228 (MQ=255)
                                                                                                                                                                                           ggggtgtgggCTTGGGGTTGCAGCGTTTTCCGGTCTTTAAGCCCGCTATCGTTTAGCGTCAGGTTGAGGTTATCGCCTTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCt                              <  2:353044/220‑1 (MQ=255)
                                                                                                                                                                                            agGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGGGTCCGCGGAACCGGTGTTGCGTTGGCTGAACTTCTCCACCCCtc                             >  1:467786/1‑227 (MQ=255)
                                                                                                                                                                                            agGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCACGCGAATGGTGTCCGGGGAACCGGTGGGGCCGTGGCTGACCTTATGCACCCCtt                             >  1:631695/1‑228 (MQ=255)
                                                                                                                                                                                              gCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGGGTCCGGTGAACCGAGGTCGCCGTGGCTGAACTTATGCACCCCTTc                            >  1:682187/1‑227 (MQ=255)
                                                                                                                                                                                               cGTCAGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGAACCCCtcagc                          >  1:561407/1‑224 (MQ=255)
                                                                                                                                                                                                 tCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTCCGGGGAACCGGCGGTGCCGTGGCTGAACTTATTCACCCCTCcgcg                         >  1:168298/1‑227 (MQ=255)
                                                                                                                                                                                                 tCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGGGTCCGGGGAACCGGTGGTGCCGTGGCGGAAATTATGCACCCCTTggcg                         >  1:155970/1‑227 (MQ=255)
                                                                                                                                                                                                              gCAGCGTTTTCCGGTAGTTAAGCACGCTAGCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCgcgc           <  1:35822/228‑1 (MQ=255)
                                                                                                                                                                                                                  gtttttCCGGTTTTTCAGCCGGTTGTCGTTACGCGTCAGGTTGAGGTTATCGCTGTTGCTCATATGGAAGCTCGGGGTCATATTTTTCCCCGGAGCTTTCGCGGTAACGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGt       <  2:168298/226‑1 (MQ=255)
                                                                                                                                                                                                                    ttttCCGGTAGTCAAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTc                                                                                                                                       >  2:205504/1‑98 (MQ=255)
                                                                                                                                                                                                                     tttCCGGTAGTCAAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTcg                                                                                                                                      <  1:205504/98‑1 (MQ=255)
                                                                                                                                                                                                                      ttCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGACCCCCTTCCCGCAGCTCTTTCGCGCCGGGCAg    >  1:166089/1‑227 (MQ=255)
                                                                                                                                                                                                                      ttCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCCCGAATGGTGTCCGGGGAACCGGTGGGGCCGTGGCTGACCTTCTGACCCCCTTCGCGCGGCCCTTCCGCCCCGGCCAgg   >  2:35822/1‑228 (MQ=255)
                                                                                                                                                                                                                      ttCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGAATTTCGACCGACATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTGGCGCCGGGCAg    >  1:681951/1‑227 (MQ=255)
                                                                                                                                                                                                                        ccGGTTGTTAAGCACGCTATCGGTGAGCGTTAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGt  <  1:205515/227‑1 (MQ=255)
                                                                                                                                                                                                                        ccGGTAGTTAAGCCGGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTAATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTGCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGt  <  1:438159/227‑1 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                           
CAAAAGCGAGTGACATTAATTTGGGGTTACGATATTCACCCACAAATCGAAATTATCCAGACAGGCGATAATTTCTTCCAGTTTCGCGCCATTGCCAATAATTTTGGCTTTCTTCGCTTTTACCAGATCCGCCATTTTGGCTTGTCCGGTCAGCACGGCGTGCAGATCTTCACGGCTGATGTAGAAAGAGGCGTCGGCTTGCGGTTGCAGCGTTTTCCGGTAGTTAAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGT  >  NC_000913/4304597‑4305039

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