Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 568,067 G→A 94.3% pseudogene (33/213 nt) renD → pseudogene, DLP12 prophage;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913568,0670GA94.3% 2102.4 / 54.8 806pseudogene (33/213 nt)renDpseudogene, DLP12 prophage;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base G (22/25);  new base A (400/354);  total (426/380)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

Only 100 of 855 total aligned reads displayed.

GAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGG  >  NC_000913/567921‑568205
                                                                                                                                                  |                                                                                                                                          
gttcggactgaaaggagcaagtgtatgccgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGcacg                                                                                                                                          >  4:193074‑M2/116‑149 (MQ=255)
gttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGcacg                                                                                                                                          <  2:277240‑M2/34‑1 (MQ=255)
 ttcggactgaaaggagcaagtgtatgccgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCGAAAATGGCACgg                                                                                                                                         >  7:127911‑M2/115‑149 (MQ=255)
 ttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattaattacctggggacgcataatagcttctgtgcgccggacgttgccgagctaacaggcccaacagtacACAGCATAAATCAGGCAGCGGCTAAAATGGCACgg                                                                                                                                         >  6:296450‑M2/115‑149 (MQ=255)
  tcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACggg                                                                                                                                        >  4:149435‑M2/114‑149 (MQ=255)
   cggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACggg                                                                                                                                        <  4:333552‑M2/36‑1 (MQ=255)
   cggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACggg                                                                                                                                        <  6:112223‑M2/36‑1 (MQ=255)
    ggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAACTCCGGCCGCGGCTAAAATGGCACggg                                                                                                                                        >  2:302655‑M2/112‑147 (MQ=255)
    ggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGc                                                                                                                                       <  1:40403‑M2/37‑1 (MQ=255)
     gactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCCGGCCGCGGCTAAAATGGCACGGGCAg                                                                                                                                     >  6:217713‑M2/111‑149 (MQ=255)
     gactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAg                                                                                                                                     >  5:102900‑M2/111‑149 (MQ=255)
     gactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAg                                                                                                                                     >  6:233532‑M2/111‑149 (MQ=255)
      actgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAgg                                                                                                                                    >  3:217762‑M2/110‑149 (MQ=255)
      actgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAgg                                                                                                                                    >  7:356666‑M2/110‑149 (MQ=255)
       ctgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCaggt                                                                                                                                   >  5:190023‑M2/109‑148 (MQ=255)
        ttaaaggagcaagtgtatgacgtgcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGcaggtc                                                                                                                                  <  5:279470‑M2/42‑3 (MQ=255)
        tgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGcaggtc                                                                                                                                  >  8:327834‑M2/108‑147 (MQ=255)
        tgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgaataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCAGCGGCTAAAATGGCACGGGcaggtc                                                                                                                                  >  5:345343‑M2/108‑147 (MQ=255)
          aaaggagcaagtgtatgacgggcaaagaggccattattcaatacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGcaggtc                                                                                                                                  <  8:236729‑M2/42‑3 (MQ=255)
          aaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggatgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGggcaggtctt                                                                                                                                >  4:312632‑M2/106‑145 (MQ=255)
          aaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCCCGggcaggtctt                                                                                                                                >  8:203417‑M2/106‑145 (MQ=255)
          aaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGggcaggtctt                                                                                                                                >  2:18897‑M2/106‑145 (MQ=255)
          aaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGggcaggtctt                                                                                                                                >  6:12984‑M2/106‑145 (MQ=255)
          aaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGggcaggtctt                                                                                                                                >  6:279118‑M2/106‑145 (MQ=255)
          aaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGggcaggtctt                                                                                                                                >  7:116177‑M2/106‑145 (MQ=255)
          aaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATAAGCCCGCGGCTAAAATGGCACGGgcaggtca                                                                                                                                 >  2:321865‑M2/106‑145 (MQ=255)
          aaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgcagcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGtgcaggtctt                                                                                                                                >  3:51536‑M2/106‑145 (MQ=255)
          aaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgacgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGggcaggtctt                                                                                                                                >  5:161595‑M2/106‑145 (MQ=255)
          aaaggagcaagtgtaggacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcacaacagtaaCCAGCATAAATCACGCCGCGGCTAAAATGGCACGggcaggtctt                                                                                                                                >  1:341366‑M2/106‑145 (MQ=255)
           aaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcatattagcttctctgcgcgggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACgggcaggtcttc                                                                                                                               <  1:145407‑M2/45‑6 (MQ=255)
           aaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACgggcaggtcttc                                                                                                                               >  7:65716‑M2/105‑144 (MQ=255)
           aaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACgggcaggtcttc                                                                                                                               <  8:70710‑M2/45‑6 (MQ=255)
           aaggaccaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataaaagcttctgtgcgacggacgttgccgcgctaacaggcccaacagaaaCCAGCAAAAATCAGGCCGCGGCTAAAATGGCACgggcaggtcttc                                                                                                                               >  4:362626‑M2/105‑144 (MQ=255)
            aggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCAcgggcaggtcttct                                                                                                                              <  3:329354‑M2/46‑7 (MQ=255)
             ggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCacgggcaggtcttctg                                                                                                                             <  6:196944‑M2/47‑8 (MQ=255)
             ggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCacgggcaggtcttctg                                                                                                                             >  3:111233‑M2/103‑142 (MQ=255)
             ggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCacggccaggtcttctg                                                                                                                             <  5:323693‑M2/47‑8 (MQ=255)
                gcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAGTGGCACGGGCAGGTCTTCTGgtt                                                                                                                          >  7:86771‑M2/100‑148 (MQ=255)
                gcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGcaggtc                                                                                                                                  <  6:382629‑M2/42‑3 (MQ=255)
                gcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGt                                                                                                                           >  8:25347‑M2/100‑148 (MQ=255)
                gcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGt                                                                                                                           >  1:363480‑M2/100‑148 (MQ=255)
                gcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGt                                                                                                                           >  4:129713‑M2/100‑148 (MQ=255)
                 caagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTggtta                                                                                                                         >  8:38225‑M2/99‑147 (MQ=255)
                  aagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaactggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCtggttat                                                                                                                        <  2:147429‑M2/52‑4 (MQ=255)
                  aagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCtggttat                                                                                                                        >  7:327716‑M2/98‑146 (MQ=255)
                   tgtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTctggttatc                                                                                                                       <  3:155115‑M2/53‑5 (MQ=255)
                   agtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctatcaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTctggttatc                                                                                                                       >  3:190263‑M2/97‑145 (MQ=255)
                   agtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCacgggcaggtcttctg                                                                                                                             <  7:48681‑M2/47‑8 (MQ=255)
                   agtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgacgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTctggttatc                                                                                                                       <  2:357695‑M2/53‑5 (MQ=255)
                    gtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTctggttatc                                                                                                                       >  8:34604‑M2/96‑144 (MQ=255)
                     tgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTtctggttatcg                                                                                                                      >  1:304243‑M2/95‑143 (MQ=255)
                      gtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTtctggttatcg                                                                                                                      >  5:232352‑M2/94‑142 (MQ=255)
                       tatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTcttctggttagcgaa                                                                                                                    <  4:102199‑M2/56‑8 (MQ=255)
                       tatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCtggttat                                                                                                                        >  1:18047‑M2/93‑141 (MQ=255)
                         tgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGcaggtc                                                                                                                                  >  2:42684‑M2/91‑130 (MQ=255)
                          gacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGgtt                                                                                                                          <  8:122379‑M2/50‑2 (MQ=255)
                           acgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTaa                                                                                                               <  2:384496‑M2/61‑1 (MQ=255)
                            cgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAg                                                                                                              <  7:246805‑M2/62‑1 (MQ=255)
                             gggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaagaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTcttctggttatcgaa                                                                                                                    <  5:37073‑M2/56‑8 (MQ=255)
                             gggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAgg                                                                                                             >  4:86165‑M2/87‑149 (MQ=255)
                             gggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAgg                                                                                                             >  8:39832‑M2/87‑149 (MQ=255)
                              ggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGaaggtaaggt                                                                                                            <  1:160418‑M2/64‑1 (MQ=255)
                              ggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGaaggtaaggt                                                                                                            >  6:87930‑M2/86‑149 (MQ=255)
                               gcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTc                                                                                                           <  8:160084‑M2/65‑1 (MQ=255)
                               gcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTc                                                                                                           <  3:273554‑M2/65‑1 (MQ=255)
                                caaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCt                                                                                                          <  7:303024‑M2/66‑1 (MQ=255)
                                caaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCt                                                                                                          <  7:222829‑M2/66‑1 (MQ=255)
                                caaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCt                                                                                                          >  3:206656‑M2/84‑149 (MQ=255)
                                 aaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTg                                                                                                         <  2:112002‑M2/67‑1 (MQ=255)
                                 aaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTg                                                                                                         <  1:234330‑M2/67‑1 (MQ=255)
                                 aaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTg                                                                                                         >  5:25370‑M2/83‑149 (MQ=255)
                                  aagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTgg                                                                                                        >  6:170168‑M2/82‑149 (MQ=255)
                                  aagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTgg                                                                                                        <  3:359722‑M2/68‑1 (MQ=255)
                                   agaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTc                                                                                                           >  1:264488‑M2/81‑145 (MQ=255)
                                   agaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGc                                                                                                       <  2:163494‑M2/69‑1 (MQ=255)
                                   agaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGc                                                                                                       >  2:227624‑M2/81‑149 (MQ=255)
                                   agaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGc                                                                                                       <  5:7217‑M2/69‑1 (MQ=255)
                                   agaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGc                                                                                                       <  8:217513‑M2/69‑1 (MQ=255)
                                    gaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGc                                                                                                       >  6:68311‑M2/80‑148 (MQ=255)
                                    gaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGc                                                                                                       >  3:300987‑M2/80‑148 (MQ=255)
                                     aggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCg                                                                                                      >  7:385007‑M2/79‑148 (MQ=255)
                                                        ggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCggtttg                                                                                  >  1:288601‑M2/60‑147 (MQ=255)
                                                        ggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCggtttg                                                                                  >  8:29138‑M2/60‑147 (MQ=255)
                                                                      aGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAATGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAgaagaa                                                                    <  5:123993/149‑1 (MQ=255)
                                                                       gAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAAt                                                                   >  4:52625/1‑149 (MQ=255)
                                                                                 ccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGaaa                                                         <  7:188655‑M2/115‑1 (MQ=255)
                                                                                  cggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAg                                                        <  6:327137‑M2/116‑1 (MQ=255)
                                                                                    gacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAaagat                                                      >  4:252453‑M2/32‑147 (MQ=255)
                                                                                       gttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTa                                                                               <  1:389259‑M2/93‑1 (MQ=255)
                                                                                              atatTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAAt                                            >  1:106493/1‑149 (MQ=255)
                                                                                                               tCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAgtg                           <  3:90538/149‑1 (MQ=255)
                                                                                                               gtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGttt                                                                                   <  1:73879/91‑2 (MQ=255)
                                                                                                                 aaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAggga                                                                        >  5:54497/1‑102 (MQ=255)
                                                                                                                    cAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGa                                                    <  3:236636/119‑1 (MQ=255)
                                                                                                                      gCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGgt                                                                                                >  4:266361/1‑73 (MQ=255)
                                                                                                                      gCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCag                                                                           <  4:219058/94‑1 (MQ=255)
                                                                                                                           aaTCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCaggg                                                                         <  2:74518/91‑3 (MQ=255)
                                                                                                                            aTCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGggaagga                                                           >  7:155290/1‑104 (MQ=255)
                                                                                                                               aGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGa                                                    >  7:278153/1‑108 (MQ=255)
                                                                                                                                  ccgCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGgtgt                                                                                              >  5:214480/3‑63 (MQ=255)
                                                                                                                                                  |                                                                                                                                          
GAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGG  >  NC_000913/567921‑568205

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: