Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 566,519 T→C 100% pseudogene (144/312 nt) peaD → pseudogene; DLP12 prophage; replication protein family;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913566,5190TC98.4% 7091.4 / 29.8 2496pseudogene (144/312 nt)peaDpseudogene; DLP12 prophage; replication protein family;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base T (23/20);  new base C (1283/1162);  total (1310/1186)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.70e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 2497 total aligned reads displayed.

TACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAACTGACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGAC  >  NC_000913/566376‑566667
                                                                                                                                               |                                                                                                                                                    
   gAGTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGAt                                                                                                                                                >  5:355434/1‑147 (MQ=255)
   gAGTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGa                                                                                                                                              <  5:298668/149‑1 (MQ=255)
    aGTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGag                                                                                                                                             <  7:176625/149‑2 (MQ=255)
     gTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGagc                                                                                                                                            >  3:319360/1‑147 (MQ=255)
     gTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGagc                                                                                                                                            <  5:69475/149‑3 (MQ=255)
        ttGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCttg                                                                                                                                         >  4:177977/1‑147 (MQ=255)
        ttGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCtt                                                                                                                                          <  8:255136/148‑2 (MQ=255)
        ttGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCtag                                                                                                                                         >  5:271347/1‑147 (MQ=255)
          gCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCttgt                                                                                                                                        <  6:212446/148‑4 (MQ=255)
           ccGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCCCTACTGGCTGGTTACCAACATGTATCAGAACATGCGTGCCTATCCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCttgtc                                                                                                                                       >  1:182684/1‑144 (MQ=255)
           ccGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGTTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCttgtc                                                                                                                                       <  6:135998/148‑5 (MQ=255)
            cGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCttgtcc                                                                                                                                      >  8:227107/1‑143 (MQ=255)
             ggAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAACTGACCTg                                                                                                                                     >  6:160782/1‑148 (MQ=255)
             ggAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCttgtccat                                                                                                                                    <  1:207493/149‑8 (MQ=255)
                             tCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTAACAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAAtt                                                                                                                     >  1:366243/1‑148 (MQ=255)
                              ccGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTaa                                                                                                                   >  6:106711/1‑149 (MQ=255)
                               cGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTaa                                                                                                                   <  3:224600/148‑1 (MQ=255)
                                   tGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCgtg                                                                                                              >  2:322748/1‑149 (MQ=255)
                                   tGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCgtg                                                                                                              <  4:117577/149‑1 (MQ=255)
                                     cggagTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCgtgg                                                                                                             <  3:256663/146‑1 (MQ=255)
                                          tCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCtt                                                                                                                                          >  7:6385/1‑113 (MQ=255)
                                           cTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTgaggc                                                                                                       >  8:333543/1‑146 (MQ=255)
                                            ttATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCtt                                                                                                                                          >  1:88731/1‑111 (MQ=255)
                                            ttATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGa                                                                                                     <  2:81622/149‑1 (MQ=255)
                                              aTCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGAt                                                                                                    <  8:7525/148‑1 (MQ=255)
                                                 cGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATccc                                                                                                 <  6:295642/148‑1 (MQ=255)
                                                 cGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGGGAGGCGATccc                                                                                                 <  3:90306/148‑1 (MQ=255)
                                                     gAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGaa                                                                                             <  2:301128/148‑1 (MQ=255)
                                                     gAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAAc                                                                                            >  1:154018/1‑149 (MQ=255)
                                                     gAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAAc                                                                                            <  7:43458/149‑1 (MQ=255)
                                                     gAAATCAAACGCGCACTACTGGCTGGTTACAAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAAc                                                                                            >  1:233816/1‑149 (MQ=255)
                                                        aTCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAg                                                                                         >  7:120635/1‑149 (MQ=255)
                                                        aTCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAg                                                                                         <  4:319360/149‑1 (MQ=255)
                                                        aTCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAg                                                                                         >  1:59385/1‑149 (MQ=255)
                                                           aaaCGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTaa                                                                                      <  5:1256/148‑1 (MQ=255)
                                                            aaCGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTACCCGTGGTGAGGCGATCCCTGAACCAGTaa                                                                                      >  4:56985/1‑148 (MQ=255)
                                                              cgcgCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAaac                                                                                   <  2:13985/149‑1 (MQ=255)
                                                                 gcACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGCTCCCTGAACCAGTAAaacaac                                                                                >  1:101471/1‑149 (MQ=255)
                                                                 gcACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTaaa                                                                                     <  5:240837/144‑1 (MQ=255)
                                                                 gcACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAaacaac                                                                                <  2:225433/149‑1 (MQ=255)
                                                                     tACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTcc                                                                            >  6:75961/1‑149 (MQ=255)
                                                                     tACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTcc                                                                            >  7:377907/1‑149 (MQ=255)
                                                                        tggctggtTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGt                                                                         >  7:205868/1‑149 (MQ=255)
                                                                        tggctggtTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGt                                                                         >  7:318347/1‑149 (MQ=255)
                                                                          gctggttACCAACGTGTATCAGAACTTGCGGGCCATTGCGCTTACTGATGCGGTATTCCGCCGTAAGGCCAAAGATGAGCTTGTCCATATGACTGCGAGATTTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCATGtca                                                                       <  8:498/147‑3 (MQ=255)
                                                                          gctggttACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGtca                                                                       >  2:115849/3‑147 (MQ=255)
                                                                          gctggttACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGtca                                                                       >  6:321796/3‑147 (MQ=255)
                                                                             ggttACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATgg                                                                    >  5:45392/1‑149 (MQ=255)
                                                                              gttACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAAc                                                                                            <  2:343012/124‑1 (MQ=255)
                                                                              gttACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTaaaa                                                                                    <  2:373802/132‑1 (MQ=255)
                                                                                  ccAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCgg                                                                <  8:138226/148‑1 (MQ=255)
                                                                                   cAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTc                                                                             <  2:270199/134‑1 (MQ=255)
                                                                                   cAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCgg                                                                >  7:232418/1‑147 (MQ=255)
                                                                                   cAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCgg                                                                >  1:387808/1‑147 (MQ=255)
                                                                                   cAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCgg                                                                >  4:150448/1‑147 (MQ=255)
                                                                                   cAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCgg                                                                >  1:258048/1‑147 (MQ=255)
                                                                                   cAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTa                                                              <  3:366045/149‑1 (MQ=255)
                                                                                   cAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTa                                                              <  7:350720/149‑1 (MQ=255)
                                                                                   cAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTa                                                              <  3:341792/149‑1 (MQ=255)
                                                                                   cAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTa                                                              <  4:209335/149‑1 (MQ=255)
                                                                                   cAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATCAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTa                                                              <  7:152368/149‑1 (MQ=255)
                                                                                     acctGTATCAGAACATGCGGGCCAATGAGCTTACTGATGCGGAATTACGCAGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGa                                                            >  6:97818/4‑149 (MQ=255)
                                                                                      cctGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCgg                                                                >  4:358322/3‑144 (MQ=255)
                                                                                      cctGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGAc                                                           <  6:260262/147‑1 (MQ=255)
                                                                                      cctGTATCAGAACATCCGGGCCAATGCGCTTACTGAGGCGGAATTACGCCGGACGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGAc                                                           <  8:163063/147‑1 (MQ=255)
                                                                                      cctGTAGCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGAc                                                           <  6:53547/147‑1 (MQ=255)
                                                                                       ctGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGAcc                                                          >  3:173907/2‑149 (MQ=255)
                                                                                       ctGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGAcc                                                          <  7:260527/148‑1 (MQ=255)
                                                                                       ctGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGAcc                                                          >  4:93277/2‑149 (MQ=255)
                                                                                        tGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGAt                                                                                                    >  1:336136/1‑106 (MQ=255)
                                                                                         gTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGACCCAGTAAAACAACTTCCTGTCATGGGCGGTAGACctc                                                        >  8:312067/1‑149 (MQ=255)
                                                                                            tcAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGCTCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTa                                                      >  1:306488/2‑148 (MQ=255)
                                                                                            tcAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTa                                                      >  7:88515/2‑148 (MQ=255)
                                                                                            tcAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATCACTGCGAGAATTAACCGTGGGGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTa                                                      >  2:159517/2‑148 (MQ=255)
                                                                                             cAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTaaa                                                    >  3:285701/1‑149 (MQ=255)
                                                                                             cAGAACATGCGGGCAAAAGCGCTGGCTGACGCGCAATTACGGCGCAAGCCTGCCGATGAACTGACCTGTATGACAGCGAGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTaaa                                                    >  4:22804/1‑149 (MQ=255)
                                                                                               gAACTTGCGGGCCAATGCGCTTACTGATGCGGAATTTCGCCGTAAGGCCGCAGATGAGCTTGTCCATTTGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATc                                                  <  5:156629/149‑1 (MQ=255)
                                                                                               gAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATc                                                  <  6:374789/149‑1 (MQ=255)
                                                                                               gAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATc                                                  <  1:281404/149‑1 (MQ=255)
                                                                                               gAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATc                                                  <  6:242142/149‑1 (MQ=255)
                                                                                               gAACATGCGGGCCAATGCGCTAACTGATGCGGAATTACGCCGTAACGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAAATTCCTGTCATGGGCGGTAGACCTCTAAATc                                                  >  4:284981/1‑149 (MQ=255)
                                                                                                aaCATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCg                                                 >  8:36585/1‑149 (MQ=255)
                                                                                                 aCATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGt                                                <  2:51371/149‑1 (MQ=255)
                                                                                                  cATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTTGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTg                                               <  2:250297/149‑1 (MQ=255)
                                                                                                  cATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGACCTTGTCCATATGACTGCGAGGATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTg                                               <  3:198347/149‑1 (MQ=255)
                                                                                                    tGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGca                                             <  4:259008/149‑1 (MQ=255)
                                                                                                       gggCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAg                                          >  5:317849/1‑149 (MQ=255)
                                                                                                        ggCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATCACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAgg                                         >  8:286292/1‑149 (MQ=255)
                                                                                                                 gcttactgATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCg                                 >  7:289738/6‑148 (MQ=255)
                                                                                                                  cttactgATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCTATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGaa                               <  6:186881/145‑1 (MQ=255)
                                                                                                                   ttactgATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCg                                 >  3:146789/4‑146 (MQ=255)
                                                                                                                     actgATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCg                                 >  5:173778/2‑144 (MQ=255)
                                                                                                                     actgATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGa                             <  3:53887/147‑1 (MQ=255)
                                                                                                                      ctgATGCGGAATTACGCAGTTATGCCGCTGATGAGCTTGTCCATATGACTGAGAGATATAACCGTGGGGAGGCGATCCCTTAACCCGTAAAACAACTTCCTGTCATGGGCGGTAGAGCTCTAAATAGTGCACAGGCTCTGGCGAAGATc                           <  6:350963/149‑1 (MQ=255)
                                                                                                                       tgATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCg                          >  5:206759/1‑149 (MQ=255)
                                                                                                                        gatgcgGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGCCCTCTAAATCGTGCACAGGCTCTGGCGAAGATCg                          >  1:117898/4‑148 (MQ=255)
                                                                                                                         atgcgGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCa                        <  8:196269/147‑1 (MQ=255)
                                                                                                                            cgGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGaa                     <  8:111879/149‑1 (MQ=255)
                                                                                                                              gAATTACGCCGTAAGTCCGCAGATGAGCTAGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAAt                   <  3:11270/149‑1 (MQ=255)
                                                                                                                                              cTGCCGATGAACTGACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGa   <  2:191586/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
TACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAACTGACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGAC  >  NC_000913/566376‑566667

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: