Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,900,545 A→T 7.1% A11A (GCT→GCA yieL ← putative xylanase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,900,5450AT 7.1% 288.3 / 8.9 83A11A (GCT→GCAyieLputative xylanase
Reads supporting (aligned to +/- strand):  ref base A (34/43);  new base T (3/3);  total (37/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-02

CAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAAC  >  NC_000913/3900271‑3900785
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                
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)
                                                                                                                             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)
                                                                                                                                                          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)
                                                                                                                                                                                 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                                                                                                                                                                                                                                                 >  6:778212/1‑87 (MQ=255)
                                                                                                                                                                                             caTACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAg                                                                                                                                                                                                                                                 <  5:778212/87‑1 (MQ=255)
                                                                                                                                                                                             caTACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGcc                                                                                                                                                                                                                                               >  8:720866/1‑89 (MQ=255)
                                                                                                                                                                                             caTACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAAtttt                                                                                                                                                                                                                           <  8:885777/109‑1 (MQ=255)
                                                                                                                                                                                             caTACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAAtttt                                                                                                                                                                                                                           >  7:885777/1‑109 (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                                                                                                                                                                                                                            <  7:148908/96‑1 (MQ=255)
                                                                                                                                                                                                         ccGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAAttt                                                                                                                                                                                                                            >  8:148908/1‑96 (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                                                                                                                                                                                                                                    <  6:965537/43‑1 (MQ=255)
                                                                                                                                                                                                                                                      cccACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAagc                                                                                                                                                                                                                                    >  5:965537/1‑43 (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                                                                                                                                                                           >  6:654238/1‑81 (MQ=255)
                                                                                                                                                                                                                                                                         gCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGAtt                                                                                                                                                                           <  5:654238/81‑1 (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)
                                                                                                                                                                                                                                                                                aTAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGgt                                                                                                                                               <  2:533657/102‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 tAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTg                                                                                                                                                 <  8:999700/99‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 tAGCCAGCGTGAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTg                                                                                                                                                 >  7:999700/1‑99 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                
CAACATTGAAAAAATACTCGTACAGGTTGCCTTTCAGGATAGGTGTGCGCCACGACCAGACTCCTGCTTCGTCTTTGGTCATCGGGTGAATATTGTCCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAAC  >  NC_000913/3900271‑3900785

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