Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 566,468 C→T 100% pseudogene (93/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,4680CT98.3% 5285.2 / 21.9 1851pseudogene (93/312 nt)peaDpseudogene; DLP12 prophage; replication protein family;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (17/16);  new base T (963/848);  total (983/868)
Fisher's exact test for biased strand distribution p-value = 8.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 1852 total aligned reads displayed.

GCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAAC‑TGACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTC  >  NC_000913/566321‑566613
                                                                                                                                                   |                                                                                                                                                  
                                                tatggtttACGAGTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTTCTGATGCGGAATTACGCCGTAAgg                                                                                                   >  1:273941/8‑149 (MQ=255)
                                                  tggtttACGAGTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGcc                                                                                                 >  5:248931/6‑148 (MQ=255)
                                                    gtttACGAGTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAATTCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTTAGGcag                                                                                                >  6:247131/4‑146 (MQ=255)
                                                    gtttACGAGTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGccg                                                                                                >  6:59632/4‑146 (MQ=255)
                                                    gtttACGAGTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGccg                                                                                                <  1:352418/145‑3 (MQ=255)
                                                    gtttACGAGTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGc                                                                                               <  7:160875/146‑1 (MQ=255)
                                                     tttACGAGTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGca                                                                                              <  3:346635/147‑2 (MQ=255)
                                                         cGAGTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGAtgcg                                                                                                                   <  4:89194/124‑1 (MQ=255)
                                                         cGAGTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGAt                                                                                           >  6:296022/1‑148 (MQ=255)
                                                          gAGTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGa                                                                                         <  7:216100/149‑1 (MQ=255)
                                                          gAGTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGa                                                                                         <  6:128271/149‑1 (MQ=255)
                                                          gAGTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGa                                                                                         <  3:87540/149‑1 (MQ=255)
                                                            gTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGagc                                                                                       >  5:374789/1‑147 (MQ=255)
                                                             tATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑t                                                                                     >  1:325085/1‑149 (MQ=255)
                                                               ttGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGAtgcg                                                                                                                   <  7:249997/118‑1 (MQ=255)
                                                               ttGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑ttg                                                                                   >  1:154712/1‑147 (MQ=255)
                                                               ttGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑tt                                                                                    <  6:213198/148‑2 (MQ=255)
                                                                tGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGTTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGAtgc                                                                                                                    <  8:385107/116‑1 (MQ=255)
                                                                tGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑ttg                                                                                   <  8:12203/148‑3 (MQ=255)
                                                                 gCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑ttgt                                                                                  <  1:219625/148‑4 (MQ=255)
                                                                  ccGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACGAACCTGTTTCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTTAGGCCGCAGCTGagcgtgtc                                                                                  >  1:241570/1‑141 (MQ=255)
                                                                  ccGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGc                                                                                                  <  2:120851/132‑1 (MQ=255)
                                                                  ccGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑ttgtcc                                                                                <  5:122891/149‑6 (MQ=255)
                                                                  ccGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑ttgtc                                                                                 <  6:81011/148‑5 (MQ=255)
                                                                  ccGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑ttgtc                                                                                 >  2:103796/1‑144 (MQ=255)
                                                                   cGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑ttgtcca                                                                               >  3:286026/1‑143 (MQ=255)
                                                                    ggAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAAC‑TGACCTg                                                                               >  6:160782/1‑148 (MQ=255)
                                                                    ggAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATgc                                                                                                                                <  4:50354/100‑1 (MQ=255)
                                                                       aGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑ttgtccat                                                                              <  3:164170/146‑8 (MQ=255)
                                                                                   aTCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCACCCTGTATCAGAACCTGCGGGCCAATGCGCGTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAAt                                                                >  4:271399/1‑148 (MQ=255)
                                                                                    tCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAAtt                                                               <  6:257923/148‑1 (MQ=255)
                                                                                    tCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAAtt                                                               >  6:329322/1‑148 (MQ=255)
                                                                                    tCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAAtt                                                               <  4:144741/148‑1 (MQ=255)
                                                                                       ggATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAAcc                                                           <  8:318982/149‑1 (MQ=255)
                                                                                          ttcggagTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCgtg                                                        <  5:176702/145‑1 (MQ=255)
                                                                                          tGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCgtg                                                        <  6:357140/149‑1 (MQ=255)
                                                                                           gCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAAC‑TGACCTGTATGACAGCGCGAATTAACCgtgg                                                       >  1:191586/1‑149 (MQ=255)
                                                                                            cggagTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCgtgg                                                       <  3:233063/146‑1 (MQ=255)
                                                                                             ggagTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCgtggt                                                      >  2:21253/2‑148 (MQ=255)
                                                                                              gagTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGAtgcg                                                                                                                   <  4:222231/87‑1 (MQ=255)
                                                                                               agTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTgag                                                   <  8:54574/149‑1 (MQ=255)
                                                                                                 tCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑tt                                                                                    >  7:6385/1‑113 (MQ=255)
                                                                                                   ttATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTTACCGTGGTGAGGCGa                                               >  2:149329/1‑149 (MQ=255)
                                                                                                   ttATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTgaggc                                                 <  8:263430/147‑3 (MQ=255)
                                                                                                   ttATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGa                                               >  6:380868/1‑149 (MQ=255)
                                                                                                      tCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGAtcc                                            >  5:281615/1‑147 (MQ=255)
                                                                                                        cGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCgtg                                                        >  1:115508/1‑135 (MQ=255)
                                                                                                        cGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCt                                          <  1:299491/149‑1 (MQ=255)
                                                                                                        cGTGGAAATCAAACGAGCACTAGTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATccc                                           <  5:89886/148‑1 (MQ=255)
                                                                                                          tGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAAtt                                                               <  4:320882/126‑1 (MQ=255)
                                                                                                          tGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAAtt                                                               >  3:320882/1‑126 (MQ=255)
                                                                                                            gAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAGTTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAAc                                      <  8:207586/149‑1 (MQ=255)
                                                                                                            gAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGaa                                       <  7:333543/148‑1 (MQ=255)
                                                                                                            gAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAAc                                      <  2:299335/149‑1 (MQ=255)
                                                                                                              aaTCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAAc                                      >  8:179151/1‑147 (MQ=255)
                                                                                                              aaTCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAAc                                      >  7:144965/1‑147 (MQ=255)
                                                                                                              aaTCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAAc                                      >  4:11270/1‑147 (MQ=255)
                                                                                                              aaTCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCa                                    <  8:234729/149‑1 (MQ=255)
                                                                                                                tcaaaCGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGt                                  <  8:302502/146‑1 (MQ=255)
                                                                                                                  aaaCGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCa                                    <  4:195801/144‑1 (MQ=255)
                                                                                                                  aaaCGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTa                                 >  7:271914/2‑148 (MQ=255)
                                                                                                                   aaCGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTaaa                               <  7:9607/149‑1 (MQ=255)
                                                                                                                    aCGCGCACTACTGGCTGGTTACCACCCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTaaaa                              <  2:70747/149‑1 (MQ=255)
                                                                                                                    aCGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTaaaa                              <  7:19752/149‑1 (MQ=255)
                                                                                                                     cgcgCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTaaaa                              >  7:85218/1‑148 (MQ=255)
                                                                                                                     cgcgCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAaac                             <  5:207984/149‑1 (MQ=255)
                                                                                                                      gcgcACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑ttg                                                                                   >  5:310218/1‑92 (MQ=255)
                                                                                                                        gcACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGCTCCCTGAACCAGTAAaacaac                          >  1:101471/1‑149 (MQ=255)
                                                                                                                        gcACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGaa                                       <  8:75621/136‑1 (MQ=255)
                                                                                                                         cactACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTaa                                                                                                     >  2:367184/4‑74 (MQ=255)
                                                                                                                         cactACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGAtcc                                            <  3:143197/127‑3 (MQ=255)
                                                                                                                           ctACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTc                       <  4:251725/148‑1 (MQ=255)
                                                                                                                           ctACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTc                       <  8:349450/148‑1 (MQ=255)
                                                                                                                           ctACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTc                       <  1:373765/148‑1 (MQ=255)
                                                                                                                            tACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCAG‑ATGACTGCGAGAATTAACCGTGGTGCGGCGATCCCTGAACCAGTAAAACAACTTcc                      >  5:154417/1‑149 (MQ=255)
                                                                                                                            tACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTcc                      >  7:377907/1‑149 (MQ=255)
                                                                                                                              ctggctggTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCt                     <  7:307405/148‑1 (MQ=255)
                                                                                                                              ctggctggTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTg                    >  3:52493/1‑149 (MQ=255)
                                                                                                                              ctggctggTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTg                    >  8:146145/1‑149 (MQ=255)
                                                                                                                               tggctggtTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGt                   <  2:392820/149‑1 (MQ=255)
                                                                                                                                 gctggttACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTaaaa                              <  3:208360/134‑1 (MQ=255)
                                                                                                                                  ctggttACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCGTATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCAt                <  3:183591/148‑1 (MQ=255)
                                                                                                                                  ctggttACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCAt                <  2:170397/148‑1 (MQ=255)
                                                                                                                                   tggttACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGAGTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACAAGTAAAACAACTTCCTGTCATg               >  5:33687/1‑149 (MQ=255)
                                                                                                                                    ggttACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATgg              <  2:268635/149‑1 (MQ=255)
                                                                                                                                     gttACCAACCTGTATCAGAACATGCGGGCCAGTGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATggg             <  7:44664/149‑1 (MQ=255)
                                                                                                                                        aCCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATggg             <  2:241105/146‑1 (MQ=255)
                                                                                                                                        aCCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCgg          >  1:73589/1‑149 (MQ=255)
                                                                                                                                          cAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCgg          >  7:112801/1‑147 (MQ=255)
                                                                                                                                          cAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTa        <  3:366045/149‑1 (MQ=255)
                                                                                                                                          cAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGAGCCCTGACCCAGTAAACCAACTGCCTGTCATGGGCgg          >  2:318820/1‑147 (MQ=255)
                                                                                                                                           accctGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCgg          >  1:389334/5‑146 (MQ=255)
                                                                                                                                           aaCCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCgg          >  3:259008/1‑146 (MQ=255)
                                                                                                                                            acctGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTGCATATGACTGCGAGAATTCACCGTGGTGAGGCGATCCCGGAACCAGTAAAACAACGTCCGGTCATGGGCGGTAGa      >  6:17866/4‑149 (MQ=255)
                                                                                                                                             cctGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGCCTGCGAGAATTAACCGTGGGGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGAc     <  4:170547/147‑1 (MQ=255)
                                                                                                                                             cctGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCgg          >  4:358322/3‑144 (MQ=255)
                                                                                                                                             cctGTATCAGAACATGCGGGCCAATGCGCTTACAGATGCGGAATTACTCCGTAAGGCCGCAGATGAGCATGTCC‑ATATAACTGAGAGAATAAAACGTAGTGATGCGATCCATGAACCACTAAAACAACTTCCAGtc                  >  6:96363/3‑135 (MQ=255)
                                                                                                                                               tGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGAt                                              <  2:336136/106‑1 (MQ=255)
                                                                                                                                               tGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGCTCCAGTAAAACAACATCCTGTCATGGGCGGTAGCCct   >  2:97937/1‑149 (MQ=255)
                                                                                                                                               tGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACct   >  5:64878/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAAC‑TGACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTC  >  NC_000913/566321‑566613

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: ‑