Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 568,052 T→C 100% pseudogene (18/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,0520TC95.5% 2929.5 / 57.8 1102pseudogene (18/213 nt)renDpseudogene, DLP12 prophage;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base T (26/25);  new base C (563/483);  total (594/508)
Fisher's exact test for biased strand distribution p-value = 7.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.08e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 1112 total aligned reads displayed.

GTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCG  >  NC_000913/567906‑568194
                                                                                                                                                  |                                                                                                                                              
agaaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGc                                                                                                                                              >  4:184906‑M2/131‑149 (MQ=255)
 gaaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCcgcg                                                                                                                                             <  6:122997‑M2/20‑1 (MQ=255)
 gaaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCcgcg                                                                                                                                             >  7:75677‑M2/130‑149 (MQ=255)
  aaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCgg                                                                                                                                            <  5:21134‑M2/21‑1 (MQ=255)
  aaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCgg                                                                                                                                            <  1:176320‑M2/21‑1 (MQ=255)
  aaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCgg                                                                                                                                            >  5:230851‑M2/129‑149 (MQ=255)
  aaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCgg                                                                                                                                            >  6:151755‑M2/129‑149 (MQ=255)
   tttAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGc                                                                                                                                           >  1:359063/1‑149 (MQ=255)
    atcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagcggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGc                                                                                                                                           >  2:55057‑M2/127‑148 (MQ=255)
     tcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTa                                                                                                                                         <  6:205745‑M2/24‑1 (MQ=255)
      caaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTa                                                                                                                                         >  6:341984‑M2/125‑148 (MQ=255)
       aaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTaaa                                                                                                                                       <  4:48476‑M2/26‑1 (MQ=255)
       aaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTaa                                                                                                                                        >  5:79205‑M2/124‑148 (MQ=255)
       aaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcactaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTaa                                                                                                                                        >  3:158515‑M2/124‑148 (MQ=255)
         agctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAAt                                                                                                                                     >  2:109157‑M2/122‑149 (MQ=255)
         agctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAAt                                                                                                                                     <  4:78476‑M2/28‑1 (MQ=255)
         agctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAAt                                                                                                                                     >  3:244441‑M2/122‑149 (MQ=255)
          gctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATg                                                                                                                                    >  8:107392‑M2/121‑149 (MQ=255)
          gctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATg                                                                                                                                    <  1:302398‑M2/29‑1 (MQ=255)
          gctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATg                                                                                                                                    <  4:216321‑M2/29‑1 (MQ=255)
          atatCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATg                                                                                                                                    >  1:154306/1‑149 (MQ=255)
           ctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATgg                                                                                                                                   >  8:291808‑M2/120‑149 (MQ=255)
           ctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATgg                                                                                                                                   >  1:282869‑M2/120‑149 (MQ=255)
            taagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGgc                                                                                                                                  <  1:136061‑M2/31‑1 (MQ=255)
            taagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgcaggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCCGGCCGCGGCTAAAATgg                                                                                                                                   >  7:272307‑M2/119‑148 (MQ=255)
             aagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATggca                                                                                                                                 >  8:44036‑M2/118‑148 (MQ=255)
              agttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGgc                                                                                                                                  <  7:245993‑M2/31‑1 (MQ=255)
              agttcggactgaaaggagcaagtgtatgacggacaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGgc                                                                                                                                  <  5:115343‑M2/31‑1 (MQ=255)
                ttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattaactggggacgcataatcgcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACgg                                                                                                                              <  3:39646‑M2/35‑1 (MQ=255)
                  cggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACggg                                                                                                                             <  8:361411‑M2/36‑1 (MQ=255)
                  cggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAATCAGGCCGCGGCTAAAATGGCACGGGc                                                                                                                            <  2:19644‑M2/36‑1 (MQ=255)
                   ggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACggg                                                                                                                             >  4:145036‑M2/112‑147 (MQ=255)
                   ggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGc                                                                                                                            <  2:237361‑M2/37‑1 (MQ=255)
                    gactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAg                                                                                                                          <  6:8486‑M2/39‑1 (MQ=255)
                    gactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGAAg                                                                                                                          <  3:309716‑M2/39‑1 (MQ=255)
                     actgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAgg                                                                                                                         >  7:356666‑M2/110‑149 (MQ=255)
                       tgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGcaggtc                                                                                                                       >  3:222287‑M2/108‑147 (MQ=255)
                       tgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgaataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCAGCGGCTAAAATGGCACGGGcaggtc                                                                                                                       >  5:345343‑M2/108‑147 (MQ=255)
                        gaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGgcaggtct                                                                                                                      >  6:229835‑M2/107‑146 (MQ=255)
                         aaaggagcaagtgtatgacgggcaaagaggctattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGggcaggtctt                                                                                                                     >  4:3072‑M2/106‑145 (MQ=255)
                         aaaggagcaagtgtatgacgggcaaagaggccattattcaatacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGcaggtc                                                                                                                       <  8:236729‑M2/42‑3 (MQ=255)
                         aaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGggcaggtctt                                                                                                                     >  2:270769‑M2/106‑145 (MQ=255)
                         aaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGcaggtc                                                                                                                       <  7:318915‑M2/42‑3 (MQ=255)
                          aaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACgggcaggtcttc                                                                                                                    <  6:319519‑M2/45‑6 (MQ=255)
                          aaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACgggcaggtcttc                                                                                                                    <  7:50524‑M2/45‑6 (MQ=255)
                          aaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACgggcaggtcttc                                                                                                                    <  8:70710‑M2/45‑6 (MQ=255)
                          aaggagcaagtgtatgacgggcaaagaggcaattattaattacctggggacgcataatagcttctgtgcgccggacgttgccgcgataagaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGggcaggtctt                                                                                                                     <  6:291881‑M2/44‑5 (MQ=255)
                           aggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCAcgggcaggtcttct                                                                                                                   >  3:182449‑M2/104‑143 (MQ=255)
                           aggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCAcgggcaggtcttct                                                                                                                   >  1:328925‑M2/104‑143 (MQ=255)
                           aggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCAAAAATCAGGCCGCGGCTAATATGGCAcgggcaggtcttct                                                                                                                   <  8:188226‑M2/46‑7 (MQ=255)
                            ggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCacgggcaggtcttctg                                                                                                                  <  2:293478‑M2/47‑8 (MQ=255)
                            ggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCacgggcaggtcttctg                                                                                                                  <  4:98180‑M2/47‑8 (MQ=255)
                            ggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCacgggcaggtcttctg                                                                                                                  >  3:111233‑M2/103‑142 (MQ=255)
                            ggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCacggccaggtcttctg                                                                                                                  <  5:323693‑M2/47‑8 (MQ=255)
                               gcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGt                                                                                                                >  6:271953‑M2/100‑148 (MQ=255)
                               gcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCAGCGGCTAAAATGGCACGGGCAGGTCTTCTGGt                                                                                                                >  5:135226‑M2/100‑148 (MQ=255)
                                caagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTggtta                                                                                                              <  1:213540‑M2/51‑3 (MQ=255)
                                caagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTggtta                                                                                                              >  8:38225‑M2/99‑147 (MQ=255)
                                 aagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaactggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCtggttat                                                                                                             <  3:98128‑M2/52‑4 (MQ=255)
                                 aagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaactggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCtggttat                                                                                                             <  2:147429‑M2/52‑4 (MQ=255)
                                  agtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCGAGCATAAATCAGGCCGCGGCTAAAATGGCAcgggcaggtcttct                                                                                                                   <  5:279559‑M2/46‑7 (MQ=255)
                                  agtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTctggttatc                                                                                                            <  2:312058‑M2/53‑5 (MQ=255)
                                  agtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgacgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTctggttatc                                                                                                            <  2:357695‑M2/53‑5 (MQ=255)
                                   gtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAg                                                                                                                          >  5:173020‑M2/96‑134 (MQ=255)
                                   gtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTctggttatc                                                                                                            >  3:28094‑M2/96‑144 (MQ=255)
                                     gtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTtctggttatcg                                                                                                           >  7:14962‑M2/94‑142 (MQ=255)
                                      tatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTcttctggttatcgaa                                                                                                         <  6:180984‑M2/56‑8 (MQ=255)
                                        tgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGcaggtc                                                                                                                       <  1:42684‑M2/42‑3 (MQ=255)
                                         gacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATgg                                                                                                                                   <  6:38796‑M2/30‑1 (MQ=255)
                                         gacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGgtt                                                                                                               <  8:122379‑M2/50‑2 (MQ=255)
                                          acgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTaa                                                                                                    <  7:185138‑M2/61‑1 (MQ=255)
                                           cgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAg                                                                                                   <  6:105328‑M2/62‑1 (MQ=255)
                                           cgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAg                                                                                                   <  5:37174‑M2/62‑1 (MQ=255)
                                           cgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAg                                                                                                   >  8:215291‑M2/88‑149 (MQ=255)
                                            gggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAgg                                                                                                  >  8:253660‑M2/87‑149 (MQ=255)
                                             ggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGaaggtaaggt                                                                                                 >  8:325944‑M2/86‑149 (MQ=255)
                                              gcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTc                                                                                                <  8:160084‑M2/65‑1 (MQ=255)
                                              gcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTc                                                                                                >  7:124486‑M2/85‑149 (MQ=255)
                                               caaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCt                                                                                               <  7:222829‑M2/66‑1 (MQ=255)
                                                aaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTg                                                                                              >  4:265212‑M2/83‑149 (MQ=255)
                                                aaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTg                                                                                              >  7:295410‑M2/83‑149 (MQ=255)
                                                aaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTg                                                                                              >  5:25370‑M2/83‑149 (MQ=255)
                                                 aagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTgg                                                                                             <  7:108215‑M2/68‑1 (MQ=255)
                                                  gtgtGGCGGCCTCAGTCCGGAAGAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTTCGTGGTTAGGATCAAGGAGATTTAATCAGGCTGCGGCTAAAATGGCGCGTGCAGGATTCCTGGTCGTTGATGGTAAGGTCTgg                                                                                             <  3:206338/148‑1 (MQ=21)
                                                  agaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGc                                                                                            <  2:163494‑M2/69‑1 (MQ=255)
                                                  agaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGc                                                                                            <  2:387830‑M2/69‑1 (MQ=255)
                                                   gaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGc                                                                                            >  6:68311‑M2/80‑148 (MQ=255)
                                                    aggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCg                                                                                           >  7:385007‑M2/79‑148 (MQ=255)
                                                                    cctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGt                                                                          <  4:122232‑M2/87‑1 (MQ=255)
                                                                      tggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGgttt                                                                        >  1:244208‑M2/61‑148 (MQ=255)
                                                                       ggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCggtttg                                                                       >  4:351718‑M2/60‑147 (MQ=255)
                                                                                                   gacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAaagat                                           >  3:345680‑M2/32‑147 (MQ=255)
                                                                                                      gttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTa                                                                    <  1:389259‑M2/93‑1 (MQ=255)
                                                                                                           cgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCGCGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACg                                    >  4:329098‑M2/24‑148 (MQ=255)
                                                                                                             atatTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAAt                                 >  1:106493/1‑149 (MQ=255)
                                                                                                                              gtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGttt                                                                        <  1:73879/91‑2 (MQ=255)
                                                                                                                                   cAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTa                                                                               <  3:72832/81‑1 (MQ=255)
                                                                                                                                   cAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGa                                         >  4:236636/1‑119 (MQ=255)
                                                                                                                                       aTAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACggg                                                     >  2:77896/1‑103 (MQ=255)
                                                                                                                                            tCAGGCCGCGGCTAAAATGACACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCa                                      <  5:54778/113‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
GTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCG  >  NC_000913/567906‑568194

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: