Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 566,513 C→T 100% pseudogene (138/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,5130CT98.6% 7508.5 / 24.0 2604pseudogene (138/312 nt)peaDpseudogene; DLP12 prophage; replication protein family;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (22/18);  new base T (1358/1196);  total (1388/1216)
Fisher's exact test for biased strand distribution p-value = 8.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 2606 total aligned reads displayed.

CCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGA‑AC‑TGACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGT  >  NC_000913/566367‑566661
                                                                                                                                                  |                                                                                                                                                      
  tatggtttACGAGTATTGCCGGAAGCGAGGCCTGTATCCCGATGCCCAGACTTATACGTGGAAATCAAACGAGCGCTACTGGCTGGTTACCCACCTGTATCAGAACATGCGGGACAAATCGCTTACTGATGCGGAATTACGCCGTAAgg                                                                                                                                                    >  4:178519/8‑149 (MQ=255)
       tttACGAGTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGccg                                                                                                                                                 >  3:347578/3‑145 (MQ=255)
           cGAGTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGAt                                                                                                                                            <  4:57351/148‑1 (MQ=255)
            gAGTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCGGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACTCCGTAAGGCCGCAGATGa                                                                                                                                          <  8:108082/149‑1 (MQ=255)
            gAGTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGa                                                                                                                                          <  8:309719/149‑1 (MQ=255)
            gAGTATTGCCGGAAGCGACGCCTGTATCCGGATGCGGTTTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGGTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGa                                                                                                                                          <  7:203766/149‑1 (MQ=255)
               tATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑t                                                                                                                                     <  1:234598/149‑1 (MQ=255)
                 ttGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGa                                                                                                                                          >  3:2614/1‑144 (MQ=255)
                 ttGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑ttg                                                                                                                                   >  1:154712/1‑147 (MQ=255)
                    ccGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑ttgtcc                                                                                                                                <  2:19618/149‑6 (MQ=255)
                      ggAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑ttgtccat                                                                                                                              <  1:207493/149‑8 (MQ=255)
                      ggAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑ttgtccat                                                                                                                              <  3:32478/149‑8 (MQ=255)
                             aggCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGAt                                                                                                                                            >  1:308851/2‑130 (MQ=255)
                               gCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGc                                                                                                                                                <  8:190582/124‑1 (MQ=255)
                                    aatccggatccggAGTCTTATCCGTGGAAAGCACACGCGCACTACTGGCGGGTTACCAACCTGTAGCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAAt                                                                                                                <  2:375704/148‑1 (MQ=255)
                                        cGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTaa                                                                                                             <  5:48965/148‑1 (MQ=255)
                                            tGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCgtg                                                                                                        <  7:63518/149‑1 (MQ=255)
                                              cggagTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCgtggt                                                                                                      <  2:325085/147‑1 (MQ=255)
                                               ggagTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCgtggtg                                                                                                     <  5:131674/148‑1 (MQ=255)
                                               ggagTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCgtg                                                                                                        <  3:23002/145‑1 (MQ=255)
                                               ggagTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGAAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCgtggt                                                                                                      >  6:299234/2‑148 (MQ=255)
                                                 agTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTgag                                                                                                   <  6:213937/149‑1 (MQ=255)
                                                  gTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTgagg                                                                                                  <  7:119401/149‑2 (MQ=255)
                                                     ttATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCATCCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGa                                                                                               <  4:153434/149‑1 (MQ=255)
                                                     ttATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑tt                                                                                                                                    >  1:88731/1‑111 (MQ=255)
                                                     ttATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTgag                                                                                                   >  1:59318/1‑145 (MQ=255)
                                                     ttATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGa                                                                                               >  3:367958/1‑149 (MQ=255)
                                                     ttATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGGGGTGAGGCGa                                                                                               >  1:118796/1‑149 (MQ=255)
                                                       aTCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGca                                                                                                                                               >  4:291541/1‑100 (MQ=255)
                                                       aTCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCg                                                                                                <  8:156960/146‑1 (MQ=255)
                                                       aTCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGa                                                                                               <  2:305022/147‑1 (MQ=255)
                                                       aTCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGAGAGGCGAtc                                                                                             >  8:136568/1‑148 (MQ=255)
                                                         ccGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGA‑AC‑TGACCTGTATGACAGCGCGAATTAAc                                                                                                            <  6:354394/132‑1 (MQ=255)
                                                           gTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGCTGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCt                                                                                          >  3:202659/1‑148 (MQ=255)
                                                           gTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGAtcc                                                                                            <  3:55458/146‑3 (MQ=255)
                                                           gTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCt                                                                                          >  1:336432/1‑148 (MQ=255)
                                                           gTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCt                                                                                          >  2:276195/1‑148 (MQ=255)
                                                            tGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTg                                                                                         >  8:219808/1‑148 (MQ=255)
                                                              gAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGaa                                                                                       <  3:287027/148‑1 (MQ=255)
                                                                 aTCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGA‑AC‑TGACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAg                                                                                   <  6:108298/149‑1 (MQ=255)
                                                                  tcaaaCGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGt                                                                                  <  8:74637/146‑1 (MQ=255)
                                                                    aaaCGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCTCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTa                                                                                 >  2:139749/2‑148 (MQ=255)
                                                                    aaaCGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTa                                                                                 >  8:36544/2‑148 (MQ=255)
                                                                     aaCGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTa                                                                                                              >  5:278019/1‑118 (MQ=255)
                                                                          gcACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAAc                                                                                                            <  6:5492/115‑1 (MQ=255)
                                                                          gcACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGCTCCCTGAACCAGTAAaacaac                                                                          >  1:101471/1‑149 (MQ=255)
                                                                           cactACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTaa                                                                                                                                                      <  1:367184/71‑1 (MQ=255)
                                                                           cactACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACt                                                                         <  3:174936/146‑1 (MQ=255)
                                                                             ctACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAaacaa                                                                           >  6:37848/2‑145 (MQ=255)
                                                                              tACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTcc                                                                      >  3:118241/1‑149 (MQ=255)
                                                                              tACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTcc                                                                      >  5:211312/1‑149 (MQ=255)
                                                                              tACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTcc                                                                      >  7:377907/1‑149 (MQ=255)
                                                                                cttgctggttACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGGGAGAATTAACCGTGTTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTg                                                                    <  5:161778/144‑1 (MQ=255)
                                                                                   gctggttACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGtca                                                                 <  1:184602/147‑3 (MQ=255)
                                                                                   gctggttACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGtca                                                                 >  2:115849/3‑147 (MQ=255)
                                                                                   gctggttACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGtca                                                                 >  6:161199/3‑147 (MQ=255)
                                                                                   gctggttACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGt                                                                   <  7:149972/145‑1 (MQ=255)
                                                                                       gttACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATggg                                                             >  7:88025/1‑149 (MQ=255)
                                                                                       gttACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATgg                                                              <  7:35055/148‑1 (MQ=255)
                                                                                          aCCAACCTGTATCAGAACATGCGGGCCAAGGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCATGAACCAGTAAAACAACTTCCTGACATGGGCgg                                                          >  2:347809/1‑149 (MQ=255)
                                                                                           ccAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCC‑TGTCC‑ATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCgg                                                          <  3:148559/148‑1 (MQ=255)
                                                                                             aaCCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAg                                                       >  5:323489/1‑149 (MQ=255)
                                                                                              acctGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCg                                                           <  4:98953/141‑1 (MQ=255)
                                                                                              acctGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGa                                                      >  5:302933/4‑149 (MQ=255)
                                                                                              acctGTATCAGAACATGCGGGCCAATGAGCTTACTGATGCGGAATTACGCAGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGa                                                      >  6:97818/4‑149 (MQ=255)
                                                                                               cctGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGAc                                                     >  5:188237/3‑149 (MQ=255)
                                                                                               cctGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGAc                                                     <  7:231735/147‑1 (MQ=255)
                                                                                                  ttataaGAACATGCCGCCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGCCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCTTGGTGAGTCGATCCCTTAACCAGTAAAACAACTTCCTGACATGCGCGGTAGACctc                                                  >  6:136386/6‑149 (MQ=255)
                                                                                                     tcAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATCACTGCGAGAATTAACCGTGGGGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTa                                                >  2:159517/2‑148 (MQ=255)
                                                                                                     tcAGAACATGAGGGCCAATGCGCTTACTGATGCGGAATTACTCCGTAAGGCCGCACATGA‑GCATGTCC‑ATATGACTGCGAGAATTAACCGTGGTGAGCCGATCCCTGCACCAGTAAAACAACTTCCTGTCATCGGCGGTAGACCTCTa                                                >  6:187935/2‑148 (MQ=255)
                                                                                                      cAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTaaa                                              >  7:43098/1‑149 (MQ=255)
                                                                                                        gAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACTCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAAATAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATc                                            >  6:327974/1‑149 (MQ=255)
                                                                                                          aCATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGt                                          <  4:23298/149‑1 (MQ=255)
                                                                                                          aCATGCGGGCCAATGCGCTTACTGATGCGGAATTACGACGTAAGGCCGCAGATGA‑GC‑TTGTCCATTTGACTGCGAGAATTAAGCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTGATGGGGGGTAGACCTCTAAATCGt                                          <  5:359395/149‑1 (MQ=255)
                                                                                                           cATGCGGGCCAATGCGCTTACTGATGCGGACTAACGCCGTAAGGCCGCAGATGA‑GC‑TGGTCCATATGACTGCGAGACTTACCCGTGGTGAGGCGATCCCTGAACCAGTAAAACACCTTCCTGTCATGGGCGGTAGACCTCTAAATCgag                                         <  4:211577/149‑3 (MQ=255)
                                                                                                           cATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTg                                         <  7:317049/149‑1 (MQ=255)
                                                                                                               cGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGca                                       >  3:93063/1‑147 (MQ=255)
                                                                                                                  gCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGt                                                         >  8:332988/1‑126 (MQ=255)
                                                                                                                    cAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGctc                                >  2:251691/1‑149 (MQ=255)
                                                                                                                     aaTGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGctct                               >  7:173283/1‑149 (MQ=255)
                                                                                                                      aTGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTg                              >  5:245226/1‑149 (MQ=255)
                                                                                                                      aTGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTg                              <  4:264065/149‑1 (MQ=255)
                                                                                                                       tGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCTATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGGAGACCTCTAAATCGTGCACAGGCTCTgg                             <  1:283767/149‑1 (MQ=255)
                                                                                                                       tGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTgg                             <  7:288870/149‑1 (MQ=255)
                                                                                                                       tGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATAAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTgg                             >  6:24495/1‑149 (MQ=255)
                                                                                                                         cgcttactgATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGc                            <  6:76540/142‑1 (MQ=255)
                                                                                                                           cttactgATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATTGGCGTTAGACCTCTAAATCGTGCACAGTCTCTGGCGaa                         <  2:233816/145‑1 (MQ=255)
                                                                                                                           cttactgATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCg                           >  1:258878/5‑147 (MQ=255)
                                                                                                                           cttactgATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGaa                         <  4:321455/145‑1 (MQ=255)
                                                                                                                            ttactgATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCg                           >  7:259754/4‑146 (MQ=255)
                                                                                                                              actgATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTcc                                                                      >  8:283003/2‑101 (MQ=255)
                                                                                                                               ctgATGCGGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATc                     <  2:52506/149‑1 (MQ=255)
                                                                                                                                  atgcgGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCg                    <  8:48287/145‑1 (MQ=255)
                                                                                                                                   tgcgGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGaa                         <  5:75460/140‑1 (MQ=255)
                                                                                                                                   tgcgGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCg                    >  2:83356/2‑146 (MQ=255)
                                                                                                                                   tgcgGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAg                 <  4:36116/148‑1 (MQ=255)
                                                                                                                                    gcgGAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGt                                          >  6:264197/1‑123 (MQ=255)
                                                                                                                                      ggAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAATACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGaaa              <  1:38727/149‑1 (MQ=255)
                                                                                                                                       gAATTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAAt             >  4:125756/1‑149 (MQ=255)
                                                                                                                                        aaTTACGCCGTAAGGCCGCAGATGA‑GC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAAt             >  8:244922/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
CCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGA‑AC‑TGACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGT  >  NC_000913/566367‑566661

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