Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 568,140 T→C 70.4% pseudogene (106/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,1400TC70.4% 149.2 / 75.1 128pseudogene (106/213 nt)renDpseudogene, DLP12 prophage;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base T (17/21);  new base C (38/52);  total (55/73)
Fisher's exact test for biased strand distribution p-value = 8.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 131 total aligned reads displayed.

GAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTC  >  NC_000913/567992‑568285
                                                                                                                                                    |                                                                                                                                                 
gAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAAt                                                                                                                                                   <  1:245830/149‑1 (MQ=255)
gAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGGTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAAt                                                                                                                                                   <  2:302802/149‑1 (MQ=255)
     tCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATggg                                                                                                                                                >  7:203872/1‑147 (MQ=255)
          gcggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGaaa                                                                                                                                         <  5:61865‑M2/115‑1 (MQ=255)
          ccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGaaa                                                                                                                                         <  4:99109‑M2/115‑1 (MQ=255)
          ccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGaaa                                                                                                                                         <  5:280938‑M2/115‑1 (MQ=255)
          ccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGaaa                                                                                                                                         <  5:65307‑M2/115‑1 (MQ=255)
          ccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGaaa                                                                                                                                         <  3:224050‑M2/115‑1 (MQ=255)
          ccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGaaa                                                                                                                                         <  2:173204‑M2/115‑1 (MQ=255)
          ccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGaaa                                                                                                                                         <  7:148924‑M2/115‑1 (MQ=255)
          ccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGaaa                                                                                                                                         <  7:188655‑M2/115‑1 (MQ=255)
          ccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGaaa                                                                                                                                         <  7:253660‑M2/115‑1 (MQ=255)
          ccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGaaa                                                                                                                                         <  8:11972‑M2/115‑1 (MQ=255)
          ccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGaaa                                                                                                                                         <  1:319510‑M2/115‑1 (MQ=255)
          ccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGaaa                                                                                                                                         <  8:395655‑M2/115‑1 (MQ=255)
          ccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGaaa                                                                                                                                         <  4:17886‑M2/115‑1 (MQ=255)
          ccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGaaa                                                                                                                                         <  1:244536‑M2/115‑1 (MQ=255)
          ccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGaaa                                                                                                                                         <  1:150419‑M2/115‑1 (MQ=255)
           cggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGaaa                                                                                                                                         <  3:204623‑M2/115‑1 (MQ=255)
           cggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAg                                                                                                                                        <  4:216632‑M2/116‑1 (MQ=255)
           cggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAg                                                                                                                                        <  6:327137‑M2/116‑1 (MQ=255)
           cggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAg                                                                                                                                        <  8:217164‑M2/116‑1 (MQ=255)
           cggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAg                                                                                                                                        <  1:270769‑M2/116‑1 (MQ=255)
            ggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGggaagga                                                                                                                                           >  4:170711‑M2/33‑145 (MQ=255)
             gacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGACGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAaagat                                                                                                                                      >  5:218477‑M2/32‑147 (MQ=255)
             gacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAaagat                                                                                                                                      >  4:252453‑M2/32‑147 (MQ=255)
             gacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAaagat                                                                                                                                      >  8:64574‑M2/32‑147 (MQ=255)
             gacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAaagat                                                                                                                                      >  3:345680‑M2/32‑147 (MQ=255)
             gacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAaagat                                                                                                                                      >  3:186872‑M2/32‑147 (MQ=255)
             gCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGGCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGt                                                                                                                                      <  7:60394/149‑1 (MQ=255)
                     cgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCGCGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACg                                                                                                                               >  4:329098‑M2/24‑148 (MQ=255)
                     ccATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACg                                                                                                                               >  4:90538/1‑148 (MQ=255)
                       atatTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAAt                                                                                                                            >  1:106493/1‑149 (MQ=255)
                          taacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCa                                                                                                                                 >  3:235245‑M2/19‑141 (MQ=255)
                                  gTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTaa                                                                                                                 <  7:330550/149‑1 (MQ=255)
                                    aacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGa                                                                                                                        >  7:121247/7‑140 (MQ=255)
                                    aacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGa                                                                                                                        <  8:121247/134‑1 (MQ=255)
                                       aTCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAgt                                                                                                            >  7:391951/1‑149 (MQ=255)
                                       aTCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAgt                                                                                                            >  7:386640/1‑149 (MQ=255)
                                        tCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAgtg                                                                                                           <  3:90538/149‑1 (MQ=255)
                                         taaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAgat                                                                                                                                      <  5:198644/120‑3 (MQ=255)
                                         taaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAgat                                                                                                                                      >  6:198644/2‑119 (MQ=255)
                                            ccAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGAtgaac                                                                                                                   >  6:82422/1‑133 (MQ=255)
                                            ccAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGAtgaac                                                                                                                   <  5:82422/137‑5 (MQ=255)
                                             cAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGa                                                                                                                                    >  4:236636/1‑119 (MQ=255)
                                             cAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGa                                                                                                                                    <  3:236636/119‑1 (MQ=255)
                                             cAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGaac                                                                                                                            <  4:8965/127‑2 (MQ=255)
                                             cAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGCTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGaac                                                                                                                            <  4:211570/127‑2 (MQ=255)
                                              aGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGAGTATTACCGGTTTGCTACCAGGGAAGAACggg                                                                                                                                                >  3:259429/1‑106 (MQ=255)
                                                 aTAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACggg                                                                                                                                                <  1:77896/103‑1 (MQ=255)
                                                 aTAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACggg                                                                                                                                                >  2:77896/1‑103 (MQ=255)
                                                 aTAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGAtga                                                                                                                     >  5:341955/1‑128 (MQ=255)
                                                 aTAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGAtga                                                                                                                     <  6:341955/130‑3 (MQ=255)
                                                 aTAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGAt                                                                                                                       <  5:321363/128‑1 (MQ=255)
                                                     aTCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCAc                                                                                                                                <  8:24685/115‑1 (MQ=255)
                                                     aTCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCAc                                                                                                                                >  7:24685/1‑115 (MQ=255)
                                                     aTCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGggaagga                                                                                                                                           >  6:107678/1‑104 (MQ=255)
                                                     aTCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGggaagga                                                                                                                                           >  7:155290/1‑104 (MQ=255)
                                                     aTCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGggaagga                                                                                                                                           <  5:107678/104‑1 (MQ=255)
                                                     aTCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCa                                                                                                                                 <  7:315501/114‑1 (MQ=255)
                                                      tCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCcg                                                                                             <  1:324612/149‑1 (MQ=255)
                                                      tCAGGCCGCGGCTAAAATGACACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCa                                                                                                                                 >  6:54778/1‑113 (MQ=255)
                                                        aGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGa                                                                                                                                    <  8:278153/108‑1 (MQ=255)
                                                             gCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAgat                                                                                                                                      >  8:337598/1‑99 (MQ=255)
                                                               ggCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGAtg                                                                                                                      <  5:230930/115‑2 (MQ=255)
                                                               ggCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGAtg                                                                                                                      >  6:230930/1‑114 (MQ=255)
                                                                gCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGTGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACg                                                                                   <  3:36282/149‑1 (MQ=255)
                                                                 cTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACgg                                                                                  <  8:203872/149‑1 (MQ=255)
                                                                 cTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGAt                                                                                                                       <  2:350024/112‑1 (MQ=255)
                                                                 cTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGAt                                                                                                                       >  1:350024/1‑112 (MQ=255)
                                                                  tAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGACCGGGAAGGAAAGATGAGCACGAACCTGa                                                                                                                        <  5:2046/110‑1 (MQ=255)
                                                                  tAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACgg                                                                                                                                                 >  2:271841/1‑85 (MQ=255)
                                                                  tAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACgg                                                                                                                                                 <  1:271841/85‑1 (MQ=255)
                                                                  tAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGggaa                                                                                                                                              <  8:338131/88‑1 (MQ=255)
                                                                  tAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGggaa                                                                                                                                              >  7:338131/1‑88 (MQ=255)
                                                                  tAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGa                                                                                                                        >  6:2046/1‑110 (MQ=255)
                                                                    aaaTGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGAtgaac                                                                                                                   <  8:57112/113‑5 (MQ=255)
                                                                    aaaTGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGAtgaac                                                                                                                   >  7:57112/1‑109 (MQ=255)
                                                                      aTGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACggg                                                                                                                                                >  6:204458/1‑82 (MQ=255)
                                                                      aTGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACggg                                                                                                                                                <  5:204458/82‑1 (MQ=255)
                                                                           acgGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAgat                                                                                                                                      >  5:372231/2‑85 (MQ=255)
                                                                            cgGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGt                                                                       <  4:205542/149‑1 (MQ=255)
                                                                             gggCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACg                                                                                                                               <  7:240360/92‑1 (MQ=255)
                                                                             gggCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACg                                                                                                                               >  8:240360/1‑92 (MQ=255)
                                                                                    aaTCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATaaaa                                                               <  6:34575/149‑1 (MQ=255)
                                                                                           gTCGTTGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTTATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATc                                                        <  3:317294/149‑1 (MQ=255)
                                                                                                  aTGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATggg                                                  <  2:106493/148‑1 (MQ=255)
                                                                                                          gTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAga                                                                                                                                       >  2:198310/1‑54 (MQ=255)
                                                                                                               gCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAatgatg                                      <  3:204262/147‑1 (MQ=255)
                                                                                                                     ggTGTATTACCTGTTCGCTACCAGAGAAGAATGGGAAGGAACTGTGAGCACGACTCTGATTTTTAAGGAGTGTCTCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAATAACATCAATGGGAACACAATGATGAAACAgtt                              <  4:287114/149‑3 (MQ=255)
                                                                                                                                 gTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTc                   >  4:50956/1‑148 (MQ=255)
                                                                                                                                       tACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGa                                                                                                              >  6:29758/1‑51 (MQ=255)
                                                                                                                                         ccAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTgtag          >  6:234568/1‑149 (MQ=255)
                                                                                                                                          cAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCaaa                <  7:381583/142‑1 (MQ=255)
                                                                                                                                          cAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCaaa                >  8:381583/1‑142 (MQ=255)
                                                                                                                                            gagaAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAACCGGGTAGTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGCGTTCAAACTgtagta        >  1:324151/1‑148 (MQ=255)
                                                                                                                                              gaagaaTGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTgtag          >  1:308802/1‑144 (MQ=255)
                                                                                                                                              gaagaaTGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTgtag          <  2:308802/144‑1 (MQ=255)
                                                                                                                                                 gaaTGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATtc  >  8:158995/1‑149 (MQ=255)
                                                                                                                                                  aaTGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGtgagttga                       >  4:239225/1‑127 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTC  >  NC_000913/567992‑568285

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