Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 566,444 T→C 100% pseudogene (69/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,4440TC97.4% 3725.1 / 31.6 1330pseudogene (69/312 nt)peaDpseudogene; DLP12 prophage; replication protein family;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base T (18/18);  new base C (660/633);  total (679/651)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 1332 total aligned reads displayed.

GCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAAC‑TGACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTC  >  NC_000913/566296‑566592
                                                                                                                                                    |                                                                                                                                                     
                                       cTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCCACTTGTACCAGAACATGCGGGc                                                                                                                 >  6:138508/1‑148 (MQ=255)
                                                                         tatggtttACGAGTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAgg                                                                              >  3:356315/8‑149 (MQ=255)
                                                                           tggtttACGAGTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGc                                                                             <  8:221550/143‑1 (MQ=255)
                                                                            ggtttACGAGTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGccg                                                                           >  5:176186/5‑147 (MQ=255)
                                                                            ggtttACGAGTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGccg                                                                           >  1:377001/5‑147 (MQ=255)
                                                                             gtttACGAGTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGc                                                                          <  7:185909/146‑1 (MQ=255)
                                                                              tttACGAGTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACgcc                                                                                    <  1:395455/136‑3 (MQ=255)
                                                                              tttACGAGTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGca                                                                         <  3:44655/147‑2 (MQ=255)
                                                                                  cGAGTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGAt                                                                      >  8:5476/1‑148 (MQ=255)
                                                                                   gAGTATTGCCGGAAGCGATGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTTTATCAGAACATTCGGGACAATGCGCTTCCTGATGCGGAATTCCGCCGTAAGGCCGCAGAt                                                                      >  1:77103/1‑147 (MQ=255)
                                                                                   gAGTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTAc                                                                                                                                                    >  6:358317/1‑69 (MQ=255)
                                                                                   gAGTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGTGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGAt                                                                      >  7:55075/1‑147 (MQ=255)
                                                                                   gAGTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGa                                                                    <  7:90190/149‑1 (MQ=255)
                                                                                   gAGTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGa                                                                    <  3:273890/149‑1 (MQ=255)
                                                                                   gAGTATTGCCGGAAGCGACGCCTGTATCCGGATGCGGTTTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGGTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGa                                                                    <  7:203766/149‑1 (MQ=255)
                                                                                     gTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGagc                                                                  <  5:69475/149‑3 (MQ=255)
                                                                                     gTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGagc                                                                  <  3:181963/149‑3 (MQ=255)
                                                                                      tATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAAC‑t                                                                <  2:86666/149‑1 (MQ=255)
                                                                                       aTTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCCGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑tt                                                               >  4:43350/1‑148 (MQ=255)
                                                                                        ttGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCgctggc                                                                                                     <  4:133277/111‑1 (MQ=255)
                                                                                        ttGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCgctggc                                                                                                     >  3:133277/1‑111 (MQ=255)
                                                                                          gCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTatc                                                                                                                              >  4:243257/1‑82 (MQ=255)
                                                                                          gCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑ttgtc                                                            >  3:228882/1‑145 (MQ=255)
                                                                                          gCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑ttgtc                                                            >  3:158030/1‑145 (MQ=255)
                                                                                          gCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑ttgt                                                             <  1:219625/148‑4 (MQ=255)
                                                                                          gCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑ttgt                                                             <  2:298143/148‑4 (MQ=255)
                                                                                          gCCGGAAGCGAGGCCTGTATCCGGAGGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGGTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑ttgt                                                             <  5:304925/148‑4 (MQ=255)
                                                                                          gCCGGAAGCGAGGCCGGTATCCGGCTGCGGAGTCTTCTCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑ttgt                                                             <  2:120047/148‑4 (MQ=255)
                                                                                           ccGTAAGCGAGGACTGTATCCGTATGCGGAGTCTTATACGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATTCGGAATTACGCCGTAAGGCCGCAGATGAGC‑ttgtcc                                                           <  6:30554/149‑6 (MQ=255)
                                                                                           ccGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGa                                                                       >  6:244060/1‑138 (MQ=255)
                                                                                           ccGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑ttgtcc                                                           <  5:122891/149‑6 (MQ=255)
                                                                                           ccGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑ttgtc                                                            >  2:103796/1‑144 (MQ=255)
                                                                                           ccGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑ttgtc                                                            >  8:283355/1‑144 (MQ=255)
                                                                                            cGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACgcc                                                                                    <  1:240916/124‑3 (MQ=255)
                                                                                             ggAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAAC‑TGACCTg                                                          >  6:160782/1‑148 (MQ=255)
                                                                                             ggAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑ttgtcca                                                          >  7:104980/1‑142 (MQ=255)
                                                                                               aaGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGAGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑ttgtccat                                                         <  4:181529/147‑8 (MQ=255)
                                                                                                aGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑ttgtccat                                                         >  4:164170/1‑139 (MQ=255)
                                                                                                        cTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGc                                                                                                                     >  2:315714/1‑79 (MQ=255)
                                                                                                        cTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATg                                                                     >  6:82191/1‑127 (MQ=255)
                                                                                                           tATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACAAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAAt                                           >  5:373658/1‑149 (MQ=255)
                                                                                                            aTCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCACCCTGTATCAGAACCTGCGGGCCAATGCGCGTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAAt                                           >  4:271399/1‑148 (MQ=255)
                                                                                                            aTCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAAt                                           >  5:34419/1‑148 (MQ=255)
                                                                                                            aTCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAAt                                           >  3:313597/1‑148 (MQ=255)
                                                                                                             tCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAAtt                                          <  6:257923/148‑1 (MQ=255)
                                                                                                             tCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAAtt                                          <  2:258074/148‑1 (MQ=255)
                                                                                                               cGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGagc                                                                  >  4:219442/1‑121 (MQ=255)
                                                                                                               cGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAAc                                       <  2:143350/149‑1 (MQ=255)
                                                                                                               cGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAAc                                       <  2:166421/149‑1 (MQ=255)
                                                                                                               cGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAAc                                       >  2:333201/1‑149 (MQ=255)
                                                                                                                 gATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCg                                     <  1:353699/149‑1 (MQ=255)
                                                                                                                 gATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCg                                     <  7:391075/149‑1 (MQ=255)
                                                                                                                   tGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCCGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCgtg                                   <  2:6973/149‑1 (MQ=255)
                                                                                                                    gCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAAC‑TGACCTGTATGACAGCGCGAATTAACCgtgg                                  >  1:191586/1‑149 (MQ=255)
                                                                                                                     cggagTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCgtggt                                 <  1:118126/147‑1 (MQ=255)
                                                                                                                     cggagTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCgtggt                                 <  6:333598/147‑1 (MQ=255)
                                                                                                                      ggagTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCgtggt                                 >  2:393236/2‑148 (MQ=255)
                                                                                                                      ggagTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACAgtggt                                 >  1:156073/2‑148 (MQ=255)
                                                                                                                       gagTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTaa                                                                                >  5:200132/1‑101 (MQ=255)
                                                                                                                       gagTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTga                               <  7:15639/149‑1 (MQ=255)
                                                                                                                       gagTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGGga                               <  7:295085/149‑1 (MQ=255)
                                                                                                                        agTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTgag                              <  2:236979/149‑1 (MQ=255)
                                                                                                                         gTCTTCTCCGTGGAAATCAACCGCGCACTCCTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTgagg                             <  8:388041/149‑2 (MQ=255)
                                                                                                                         gTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTgagg                             <  8:30789/149‑2 (MQ=255)
                                                                                                                          tCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAAATAACCGTGGTgaggc                            >  4:88887/1‑147 (MQ=255)
                                                                                                                            ttATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGATTTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGa                          >  2:200923/1‑149 (MQ=255)
                                                                                                                            ttATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTgaggc                            >  7:263430/1‑145 (MQ=255)
                                                                                                                            ttATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGa                          >  7:313485/1‑149 (MQ=255)
                                                                                                                            ttATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGa                          >  7:234729/1‑149 (MQ=255)
                                                                                                                            ttATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGa                          >  6:380868/1‑149 (MQ=255)
                                                                                                                             tATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGc                                                                          >  8:195080/1‑101 (MQ=255)
                                                                                                                              aTCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGAtc                        >  8:107220/1‑148 (MQ=255)
                                                                                                                              aTCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGAt                         <  8:7525/148‑1 (MQ=255)
                                                                                                                               tCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑tt                                                               <  5:157946/109‑2 (MQ=255)
                                                                                                                                ccGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGc                                                                          >  3:337499/1‑98 (MQ=255)
                                                                                                                                 cGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATccc                      <  4:347578/148‑1 (MQ=255)
                                                                                                                                 cGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCt                     <  1:317945/149‑1 (MQ=255)
                                                                                                                                 cGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCt                     <  1:299491/149‑1 (MQ=255)
                                                                                                                                  gTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGCTGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCt                     >  3:202659/1‑148 (MQ=255)
                                                                                                                                   tGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTg                    >  3:256140/1‑148 (MQ=255)
                                                                                                                                   tGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGa                   <  2:370007/149‑1 (MQ=255)
                                                                                                                                   tGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATCAGC‑TTGTCCATATTACTGCGAGAATTAACCGTGGTGAGGCGATCCCTg                    >  5:243876/1‑148 (MQ=255)
                                                                                                                                     gAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGa                                                                                                  <  6:17808/69‑1 (MQ=255)
                                                                                                                                     gAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAgg                                                                              <  4:316587/89‑1 (MQ=255)
                                                                                                                                     gAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAAc                 >  8:322017/1‑149 (MQ=255)
                                                                                                                                      aaaTCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAAc                 >  8:351690/1‑148 (MQ=255)
                                                                                                                                        aTCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTACACGTGGTGAGGCGATCCCTGAACCAg              >  8:173224/1‑149 (MQ=255)
                                                                                                                                        aTCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAg              >  7:120635/1‑149 (MQ=255)
                                                                                                                                         tcaaaCGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCATGTCC‑AAATGACTGCGAGAATTAACCGTGGTGAGGCGCTCCCTGAACCAg              >  6:10891/4‑148 (MQ=255)
                                                                                                                                          caaaCGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGt             >  8:141670/3‑148 (MQ=255)
                                                                                                                                           aaaCGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTaa           <  3:7786/148‑1 (MQ=255)
                                                                                                                                           aaaCGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTa            >  7:74188/2‑148 (MQ=255)
                                                                                                                                           aaaCGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTa            >  6:64009/2‑148 (MQ=255)
                                                                                                                                            aaCGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTa                                         >  5:278019/1‑118 (MQ=255)
                                                                                                                                               gcgcACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAaac        >  1:272429/1‑148 (MQ=255)
                                                                                                                                                cgcACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAaacaa      <  7:36544/149‑1 (MQ=255)
                                                                                                                                                 gcACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGCTCCCTGAACCAGTAAaacaac     >  1:101471/1‑149 (MQ=255)
                                                                                                                                                 gcACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTATaaccac     >  4:107119/1‑149 (MQ=255)
                                                                                                                                                 gcACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAaacaac     >  6:211431/1‑149 (MQ=255)
                                                                                                                                                    ttACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAAC‑TGACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTc  >  8:245713/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAAC‑TGACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTC  >  NC_000913/566296‑566592

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