Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,650,147 G→A 5.6% pseudogene (70/765 nt) ydfE → Qin prophage; pseudogene;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,650,1470GA 5.6% 418.8 / 9.0 125pseudogene (70/765 nt)ydfEQin prophage; pseudogene;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base G (66/52);  new base A (4/3);  total (70/55)
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.77e-01

Only 100 of 126 total aligned reads displayed.

AATATCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGA  >  NC_000913/1650001‑1650295
                                                                                                                                                  |                                                                                                                                                    
aataTCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTa                                                                                                                                                    <  4:132560/149‑1 (MQ=255)
aataTCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTa                                                                                                                                                    <  2:291734/149‑1 (MQ=255)
      cGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCcaca                                                                                                                                              >  3:954959/1‑149 (MQ=255)
      cGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCcaca                                                                                                                                              >  3:814606/1‑149 (MQ=255)
         tttCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGcc                                                                                                                                           >  4:505218/1‑149 (MQ=255)
         tttCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGcc                                                                                                                                           >  1:928689/1‑149 (MQ=255)
           tCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCAt                                                                                                                                         <  2:1040271/149‑1 (MQ=255)
            cATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTACAGCCACAGcgcct                                                                                                                                        <  1:999027/149‑5 (MQ=255)
             aTTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTg                                                                                                                                       >  1:465719/1‑149 (MQ=255)
               tCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGaa                                                                                                                                     >  4:576287/1‑149 (MQ=255)
                    tatGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTgg                                                                                                                                >  1:223612/1‑149 (MQ=255)
                               tGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGcc                                                                                                                     <  3:646145/149‑1 (MQ=255)
                                    gggATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCt                                                                                                                <  2:247714/149‑1 (MQ=255)
                                       aTTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGa                                                                                                             <  3:158438/149‑1 (MQ=255)
                                       aTTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGa                                                                                                             >  1:803766/1‑149 (MQ=255)
                                       aTTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGa                                                                                                             >  3:85152/1‑149 (MQ=255)
                                        tttGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAg                                                                                                            >  1:492296/1‑149 (MQ=255)
                                            ttCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCaaa                                                                                                        >  3:419398/1‑149 (MQ=255)
                                            ttCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCaaa                                                                                                        >  3:235638/1‑149 (MQ=255)
                                            ttCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCaaa                                                                                                        <  3:108834/149‑1 (MQ=255)
                                            ttCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCaaa                                                                                                        >  3:815319/1‑149 (MQ=255)
                                            ttCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCGCTGGCTGTAGCCACAGCCATTGAAGCTGTCGACGAACGCCGTGCTCGAGCgaa                                                                                                        >  4:820192/1‑149 (MQ=255)
                                             tCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCaaaa                                                                                                       <  2:689834/149‑1 (MQ=255)
                                              cATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAAt                                                                                                      >  2:824289/1‑149 (MQ=255)
                                                  gttAAATCTGTCATGAGGTTAAAACAAAAGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCa                                                                                                  <  4:930955/148‑1 (MQ=255)
                                                   ttAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCAc                                                                                                 >  3:790731/1‑149 (MQ=255)
                                                   ttAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCAc                                                                                                 >  1:890782/1‑149 (MQ=255)
                                                   ttAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCAc                                                                                                 >  3:666971/1‑149 (MQ=255)
                                                   ttAAATCTGTCATGAGGTTAAAACAAATGAGCAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCAc                                                                                                 >  2:351693/1‑149 (MQ=255)
                                                    tAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACt                                                                                                <  2:223612/149‑1 (MQ=255)
                                                     aaaTCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTg                                                                                               <  3:73302/149‑1 (MQ=255)
                                                      aaTCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTgg                                                                                              >  3:534325/1‑149 (MQ=255)
                                                        tCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGaa                                                                                                                                     <  3:724094/108‑1 (MQ=255)
                                                        tCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGaa                                                                                                                                     >  4:724094/1‑108 (MQ=255)
                                                           gTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAAtt                                                                                         <  1:431580/149‑1 (MQ=255)
                                                           gTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAAtt                                                                                         <  3:131529/149‑1 (MQ=255)
                                                              aTGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCt                                                                                      <  2:437791/149‑1 (MQ=255)
                                                              aTGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCt                                                                                      >  2:527047/1‑149 (MQ=255)
                                                               tGAGGTTAAAAGAAATGAGTAAAGTCTTTGTTTTCGCCGCTATTCCTGACGAACTGGAAACAAGGAAAGAGGGCTCTGTGGCTGTAGCCACAGCCATTGAAGCTCTCGACGAACGCCGTGCTCGAGCAAAATTTCACAGGCAATTCCTg                                                                                     <  3:418285/149‑1 (MQ=255)
                                                               tGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTg                                                                                     >  3:257862/1‑149 (MQ=255)
                                                                    ttAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGa                                                                                                             <  1:1052933/120‑1 (MQ=255)
                                                                    ttAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACa                                                                                >  2:754914/1‑149 (MQ=255)
                                                                      aaaaCAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACAtt                                                                              <  2:1056366/149‑1 (MQ=255)
                                                                         aCAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGATCGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTAta                                                                           <  2:936510/149‑2 (MQ=255)
                                                                             aTGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCa                                                                                                  >  1:736015/1‑122 (MQ=255)
                                                                             aTGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCa                                                                                                  <  2:736015/122‑1 (MQ=255)
                                                                               gAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGTGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGgctg                                                                     >  1:478773/1‑149 (MQ=255)
                                                                                 gTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCAt                                                                                                                                         <  2:386102/79‑1 (MQ=255)
                                                                                    aaGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAg                                                                <  1:1034299/149‑1 (MQ=255)
                                                                                       tCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGAc                                                             >  4:447467/1‑149 (MQ=255)
                                                                                       tCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGAc                                                             >  4:706854/1‑149 (MQ=255)
                                                                                       tCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCCCAGGAc                                                             >  3:631630/1‑149 (MQ=255)
                                                                                             tttGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCt                                                       >  4:432411/1‑149 (MQ=255)
                                                                                                   ccgcTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATaaa                                                 <  3:576287/149‑1 (MQ=255)
                                                                                                    cgcTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAAt                                                >  1:339465/1‑149 (MQ=255)
                                                                                                     gcTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAAtt                                               <  1:85501/149‑1 (MQ=255)
                                                                                                       tATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTa                                             <  1:731706/149‑1 (MQ=255)
                                                                                                       tATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTa                                             <  4:1018921/149‑1 (MQ=255)
                                                                                                         ttCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTAtt                                           <  3:447467/149‑1 (MQ=255)
                                                                                                            cTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTc                                        >  4:280947/1‑149 (MQ=255)
                                                                                                             tGACGAACTGGCAACAAGGTAAGAAGGCGCTATGGCTGTAGCCAAAACCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGTCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGATTATAAATATATTGttt                                       >  4:900102/1‑147 (MQ=255)
                                                                                                             tGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCt                                       >  3:443958/1‑149 (MQ=255)
                                                                                                             tGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCt                                       <  4:815319/149‑1 (MQ=255)
                                                                                                                 gAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGa                                   >  3:850174/1‑149 (MQ=255)
                                                                                                                 gAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGa                                   >  1:1016516/1‑149 (MQ=255)
                                                                                                                 gAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGa                                   <  3:537262/149‑1 (MQ=255)
                                                                                                                 gAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGa                                   >  2:109791/1‑149 (MQ=255)
                                                                                                                 gAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGa                                   <  3:280947/149‑1 (MQ=255)
                                                                                                                 gAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGa                                   <  2:803766/149‑1 (MQ=255)
                                                                                                                  aaCTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAg                                  <  1:602018/149‑1 (MQ=255)
                                                                                                                  aaCTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAg                                  <  4:790731/149‑1 (MQ=255)
                                                                                                                  aaCTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAg                                  <  1:754914/149‑1 (MQ=255)
                                                                                                                     tGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATc                                                                           <  2:384267/105‑1 (MQ=255)
                                                                                                                     tGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATc                                                                           >  1:384267/1‑105 (MQ=255)
                                                                                                                         aacaaGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAAc                           <  2:1016516/149‑1 (MQ=255)
                                                                                                                         aacaaGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAAc                           >  4:637502/1‑149 (MQ=255)
                                                                                                                         aacaaGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAAc                           >  4:565926/1‑149 (MQ=255)
                                                                                                                            aaGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTg                        >  3:613184/1‑149 (MQ=255)
                                                                                                                            aaGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTg                        >  3:356647/1‑149 (MQ=255)
                                                                                                                             aGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTgcgc                                                        >  2:696113/1‑116 (MQ=255)
                                                                                                                             aGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTgcgc                                                        <  1:696113/116‑1 (MQ=255)
                                                                                                                             aGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTgg                       <  1:678149/149‑1 (MQ=255)
                                                                                                                             aGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTgg                       <  1:109791/149‑1 (MQ=255)
                                                                                                                                 aagaagGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATa                   >  4:199115/1‑149 (MQ=255)
                                                                                                                                 aagaagGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATa                   >  2:448617/1‑149 (MQ=255)
                                                                                                                                 aagaagGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATa                   >  2:527101/1‑149 (MQ=255)
                                                                                                                                      ggCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCg              >  4:1016639/1‑149 (MQ=255)
                                                                                                                                      atggcTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTgc                                                          <  4:882532/102‑1 (MQ=255)
                                                                                                                                      atggcTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTgc                                                          >  3:882532/4‑105 (MQ=255)
                                                                                                                                      atggcTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCg              >  4:524100/4‑149 (MQ=255)
                                                                                                                                      atggcTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACCGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCg              >  1:957574/4‑149 (MQ=255)
                                                                                                                                        ggcTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCgcc            >  4:189051/2‑149 (MQ=255)
                                                                                                                                           tgtgGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCgccctg         >  2:144425/1‑149 (MQ=255)
                                                                                                                                             tgGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCgccctgcc       >  3:150345/1‑149 (MQ=255)
                                                                                                                                             tgGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCgccctgcc       <  4:307495/149‑1 (MQ=255)
                                                                                                                                             tgGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCgccctgcc       <  2:615331/149‑1 (MQ=255)
                                                                                                                                               gCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCgccctgccct     >  1:683393/1‑149 (MQ=255)
                                                                                                                                                cTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTc    <  2:53099/149‑1 (MQ=255)
                                                                                                                                                 tGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCg   >  1:214383/1‑149 (MQ=255)
                                                                                                                                                  gTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGa  <  2:339465/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AATATCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGA  >  NC_000913/1650001‑1650295

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: