Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,900,542 A→T 6.7% I12I (ATT→ATA yieL ← putative xylanase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,900,5420AT 6.7% 309.2 / 8.6 88I12I (ATT→ATAyieLputative xylanase
Reads supporting (aligned to +/- strand):  ref base A (38/44);  new base T (3/3);  total (41/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

CATCAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAAC  >  NC_000913/3900268‑3900785
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   
cATCAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCa                                                                                                                                                                                                                                                     >  2:348069/1‑275 (MQ=255)
   cAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATa                                                                                                                                                                                                                                                  >  2:97974/1‑275 (MQ=255)
                    gcGTATAGGGTGCCTTTAAGTAGATGTTTGCGCCAAGCCCAGATCACTGCTTGTTCTTTGGTCATCGGGTGAACATGGTCCGGAACCGGAAAACCAACCAATACAGAGACATTTTTTGCCCCAGGGTCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAagcagc                                                                                                                                                                                                                                 <  4:366671/274‑1 (MQ=255)
                           ggTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCGGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAAttttt                                                                                                                                                                                                                          >  4:356017/1‑275 (MQ=255)
                                  ttttcggttgggTGGGCGCCCCGACCAGACCCCTGCTCCGTCTTTGGTCATCTGGGGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGGCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCa                                                                                                                                                                                                                   <  4:519974/265‑1 (MQ=255)
                                  tttCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCa                                                                                                                                                                                                                   <  4:348068/275‑1 (MQ=255)
                                  tttCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCa                                                                                                                                                                                                                   <  2:483856/275‑1 (MQ=255)
                                          tAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTaa                                                                                                                                                                                                             <  4:448647/273‑1 (MQ=255)
                                                  cgcCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTTCAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTTACCTGAGTCACATACTGTTTCACCGGGCTAGTGCGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCttt                                                                                                                                                                                                    >  2:343137/1‑274 (MQ=255)
                                                       cGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACGGAGACATTTTTTGCCCCAGGGGCAAAGTAGCCAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTaa                                                                                                                                                                                              >  2:95324/1‑275 (MQ=255)
                                                                          tCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTtata                                                                                                                                                                                                                       >  4:129368/1‑231 (MQ=255)
                                                                              tGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCa                                                                                                                                                                                                                   <  2:129369/231‑1 (MQ=255)
                                                                                           aaTATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACCTACTGTTTCACCGGGATAGTGTGTGCCGGGCTGGCTGGCCTATCCGCTGCTATGGCCCCCTGAGCAGAAATACCGCTGGCAATAACCAGCGTTAAAGCAGCAATTTTTATATTCATCGTACTATTACTTTCATTAAGTTCCATTACCACCAGATTTCAGACTGGGCACCGa                                                                                                                                                           >  1:145098/1‑274 (MQ=255)
                                                                                                 gTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCaaa                                                                                                                                                    <  4:97973/275‑1 (MQ=255)
                                                                                                  tCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCTCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGCCAGCAAAc                                                                                                                                                   >  2:150702/1‑275 (MQ=255)
                                                                                                  tCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAAc                                                                                                                                                   >  2:46581/1‑275 (MQ=255)
                                                                                                       aaccggaacACCCACCACTAACGAGACATTTTTTGCCCCAGGGGCAAATTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGgt                                                                                                                                               <  4:67949/274‑1 (MQ=255)
                                                                                                       aaccggaaaACCCACCACTACAGGGCAATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGgt                                                                                                                                               <  4:299604/274‑1 (MQ=255)
                                                                                                                 cccaccACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAAtttt                                                                                                                                                                                                                           >  4:225364/1‑188 (MQ=255)
                                                                                                                     ccacTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTtata                                                                                                                                                                                                                       <  2:225365/188‑1 (MQ=255)
                                                                                                                            agagACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAACGCAGCAATTTTTATATTCATCGTACTATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAATgg                                                                                                                          >  1:36410/1‑274 (MQ=255)
                                                                                                                                tatatttttGCCCCGGGGGCAAAGTATCANNNGGTCACGCTGTTATCGGCATTGACCTGAGTCACCTACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGTCCCACTGGGCAGAAATACCGCTGGTAATAACCAGCGTTAAATCAGCAATGTTTATATTCATCGTAATATTCCTTTCATTAAGTACCATTACCACCAGATTTCAGCCTGAGCACCGACAGCAAACTGGTCATCTTTATTGTCTTCAAAGGTGAAg                                                                                                                     <  3:145098/271‑1 (MQ=255)
                                                                                                                                aCATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAAt                                                                                                                                                                                                                                                   >  4:564051/1‑149 (MQ=255)
                                                                                                                                aCATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTAAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTTTCTTCAAAGGTGAAg                                                                                                                     >  2:1091/1‑275 (MQ=255)
                                                                                                                                   ttttttGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGcc                                                                                                                   >  2:173056/1‑274 (MQ=255)
                                                                                                                                    tttttGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGcc                                                                                                                                                                                                                                               <  2:564052/149‑1 (MQ=255)
                                                                                                                                     ttttGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATc                                                                                                                >  2:255176/1‑275 (MQ=255)
                                                                                                                                             cAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATTGGCATTAACCGGAGTACCATACGGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATGTCCGCAGCTAAGGCCCACTGAGCAGAAATACCGCTGGCAATAACCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTATCATTAAGTTCCATTACCACCAGATGTCAGCCTGGGCACCGACAGCAAACTAGTCGTCGTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCg                                                                                                         <  4:343136/274‑1 (MQ=255)
                                                                                                                                             cAGGGGAAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCg                                                                                                         <  4:46580/274‑1 (MQ=255)
                                                                                                                                                     aaGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCaa                                                                                                                                                                                                                               >  8:1086697/1‑148 (MQ=255)
                                                                                                                                                        tAGGGAAAGGTCACGCTGGTATCGGCAGTGACCTGTGGCACATACGTTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCTGCTGTTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCt                                                                                             <  4:1091/275‑1 (MQ=255)
                                                                                                                                                        tAGGGAAAGGGCACGCTGTTAGCGGCATTGACAGGAGTCCAATCCTGTTTAACCGGGATATTGCGTGCCGGGCTGGCTGGCATATCCGCTGCGATGGCCCACTGAGCAAAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGGTTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGc                                                                                              <  4:458255/274‑1 (MQ=255)
                                                                                                                                                         atcGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCt                                                                                             <  1:113126/272‑1 (MQ=255)
                                                                                                                                                           cGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTaa                                                                                                                                                                                                                                       >  7:1052848/1‑134 (MQ=255)
                                                                                                                                                           cGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTaa                                                                                                                                                                                                                                       <  8:1052848/134‑1 (MQ=255)
                                                                                                                                                             aaaGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAAt                                                                                                                                                                                                                                                   <  8:440954/120‑1 (MQ=255)
                                                                                                                                                             aaaGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAAt                                                                                                                                                                                                                                                   >  7:440954/1‑120 (MQ=255)
                                                                                                                                                             aaaGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGc                                                                                                                                                                                                                                                <  7:762855/123‑1 (MQ=255)
                                                                                                                                                             aaaGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGc                                                                                                                                                                                                                                                >  8:762855/1‑123 (MQ=255)
                                                                                                                                                                      gCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTa                                                                                                                                                                                                              <  7:1086697/148‑1 (MQ=255)
                                                                                                                                                                         gTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTa                                                                                                                                                                                                                                        <  8:1011318/119‑1 (MQ=255)
                                                                                                                                                                         gTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTa                                                                                                                                                                                                                                        >  7:1011318/1‑119 (MQ=255)
                                                                                                                                                                              cGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAAt                                                                                                                                                                                                                                                   <  8:883594/103‑1 (MQ=255)
                                                                                                                                                                              cGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAAt                                                                                                                                                                                                                                                   >  7:883594/1‑103 (MQ=255)
                                                                                                                                                                                    tGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTgg                                                                                                                                                                  >  4:455751/1‑178 (MQ=255)
                                                                                                                                                                                        cTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCAc                                                                                                                                                              <  2:455752/178‑1 (MQ=255)
                                                                                                                                                                                            gTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGAcgc                                                         <  4:255175/275‑1 (MQ=255)
                                                                                                                                                                                              cacaTACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAAtt                                                                                                                                                                                                                             <  7:1133033/109‑1 (MQ=255)
                                                                                                                                                                                              cacaTACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAAtt                                                                                                                                                                                                                             >  8:1133033/1‑109 (MQ=255)
                                                                                                                                                                                                caTACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAg                                                                                                                                                                                                                                                 <  5:778212/87‑1 (MQ=255)
                                                                                                                                                                                                caTACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAg                                                                                                                                                                                                                                                 >  6:778212/1‑87 (MQ=255)
                                                                                                                                                                                                caTACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGcc                                                                                                                                                                                                                                               >  8:720866/1‑89 (MQ=255)
                                                                                                                                                                                                caTACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAAtttt                                                                                                                                                                                                                           >  7:885777/1‑109 (MQ=255)
                                                                                                                                                                                                caTACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAAtttt                                                                                                                                                                                                                           <  8:885777/109‑1 (MQ=255)
                                                                                                                                                                                                caTACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGACCCACTGAGCAGAAATACCGCTGGCAATAGcc                                                                                                                                                                                                                                               <  7:720866/89‑1 (MQ=255)
                                                                                                                                                                                                  tACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTaaa                                                                                                                                                                                                                                      >  5:779078/1‑96 (MQ=255)
                                                                                                                                                                                                  tACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTaaa                                                                                                                                                                                                                                      <  6:779078/96‑1 (MQ=255)
                                                                                                                                                                                                  tACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAAttt                                                                                                                                                                                                                            >  5:332899/1‑106 (MQ=255)
                                                                                                                                                                                                  tACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCACACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAAttt                                                                                                                                                                                                                            <  6:332899/106‑1 (MQ=255)
                                                                                                                                                                                                      gTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATAACGCTGGCTATTGCCAGCg                                                                                                                                                                                                                                           >  6:839784/1‑87 (MQ=255)
                                                                                                                                                                                                      gTTTCAACGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATAACGCTGGCTATTGCCAGCg                                                                                                                                                                                                                                           <  5:839784/87‑1 (MQ=255)
                                                                                                                                                                                                            ccGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCTATTGCCAGCGGTattt                                                                                                                                                                                                                                     <  8:148394/87‑4 (MQ=38)
                                                                                                                                                                                                            ccGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCTATTGCCAGCGGTattt                                                                                                                                                                                                                                     >  7:148394/1‑84 (MQ=38)
                                                                                                                                                                                                            ccGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAAtttt                                                                                                                                                                                                                           <  6:729714/97‑1 (MQ=255)
                                                                                                                                                                                                            ccGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAAtttt                                                                                                                                                                                                                           >  5:729714/1‑97 (MQ=255)
                                                                                                                                                                                                            ccGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAAttt                                                                                                                                                                                                                            >  8:148908/1‑96 (MQ=255)
                                                                                                                                                                                                            ccGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAAttt                                                                                                                                                                                                                            <  7:148908/96‑1 (MQ=255)
                                                                                                                                                                                                             cGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCCGGGTCACCATACTGGt                                        >  4:504242/1‑275 (MQ=255)
                                                                                                                                                                                                                gATAGTTGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTAtt                                     <  4:95323/275‑1 (MQ=255)
                                                                                                                                                                                                                       ggTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGGGCAGAAATACCGCTGGCAATa                                                                                                                                                                                                                                                  <  8:700664/63‑1 (MQ=255)
                                                                                                                                                                                                                       ggTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATa                                                                                                                                                                                                                                                  >  7:700664/1‑63 (MQ=255)
                                                                                                                                                                                                                        gTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTTGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTaa                             <  2:58597/275‑1 (MQ=255)
                                                                                                                                                                                                                        gTGCCGGGATGGCTGGCATATCCGCAGCAATGGGCCACTGGGCAGTAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAGGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTaa                             <  4:150701/275‑1 (MQ=255)
                                                                                                                                                                                                                                         atatCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGt             <  2:356018/274‑1 (MQ=255)
                                                                                                                                                                                                                                                 tgctATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGttt                                                                                                      >  4:98186/1‑177 (MQ=255)
                                                                                                                                                                                                                                                     aTGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCa                                                                                                  <  2:98187/177‑1 (MQ=255)
                                                                                                                                                                                                                                                         cccACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAagc                                                                                                                                                                                                                                    >  5:965537/1‑43 (MQ=255)
                                                                                                                                                                                                                                                         cccACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAagc                                                                                                                                                                                                                                    <  6:965537/43‑1 (MQ=255)
                                                                                                                                                                                                                                                         cccACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTc                                                                                                                                                                                          >  7:400991/1‑85 (MQ=255)
                                                                                                                                                                                                                                                         cccACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTc                                                                                                                                                                                          <  8:400991/85‑1 (MQ=255)
                                                                                                                                                                                                                                                                          ccGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTta      >  4:147096/1‑248 (MQ=255)
                                                                                                                                                                                                                                                                            gCTGGCTATTGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGAttt                                                                                                                                                                          >  7:818598/1‑82 (MQ=255)
                                                                                                                                                                                                                                                                            gCTGGCTATTGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGAttt                                                                                                                                                                          <  8:818598/82‑1 (MQ=255)
                                                                                                                                                                                                                                                                            gCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGAtt                                                                                                                                                                           <  5:654238/81‑1 (MQ=255)
                                                                                                                                                                                                                                                                            gCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGAtt                                                                                                                                                                           >  6:654238/1‑81 (MQ=255)
                                                                                                                                                                                                                                                                              tGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAAc  <  2:147097/248‑1 (MQ=255)
                                                                                                                                                                                                                                                                               ggCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAAc                                                                                                                                                   >  4:533656/1‑102 (MQ=255)
                                                                                                                                                                                                                                                                               ggCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGa                                                                                                                       >  5:355979/1‑130 (MQ=255)
                                                                                                                                                                                                                                                                               ggCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGa                                                                                                                       <  6:355979/130‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   
CATCAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAAC  >  NC_000913/3900268‑3900785

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