| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description | 
| RA | NC_000913 | 1,572,928 | C→G | 10.0% | A759P (GCA→CCA) | pqqL ← | putative periplasmic M16 family zinc metalloendopeptidase | 
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,572,928 | 0 | C | G | 10.0% | 342.1 / 22.6 | 100 | A759P (GCA→CCA) | pqqL | putative periplasmic M16 family zinc metalloendopeptidase | 
| Reads supporting (aligned to +/- strand): ref base C (46/44); new base G (5/5); total (51/49) | |||||||||||
| 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.98e-01 | |||||||||||
| AACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACAT  >  NC_000913/1572785‑1573074                                                                                                                                               |                                                                                                                                                  aaCGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGGcaca                                                                                                                                                >  3:810464/1‑148 (MQ=255)  cGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGGCACAgg                                                                                                                                              <  4:810464/148‑1 (MQ=255)                gATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCca                                                                                                                                                  >  1:540396/1‑130 (MQ=255)                gATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCca                                                                                                                                                  <  2:540396/130‑1 (MQ=255)                    ccTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGGCACAGGttt                                                                                                                                           >  1:1035178/1‑132 (MQ=255)                    ccTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGGCACAGGttt                                                                                                                                           <  2:1035178/133‑2 (MQ=255)                         tCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTg                                                                                                                                     <  4:888370/134‑1 (MQ=255)                         tCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTg                                                                                                                                     >  3:888370/1‑134 (MQ=255)                             gAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCcac                                                                                                                                                 >  4:696992/1‑118 (MQ=255)                             gAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCcac                                                                                                                                                 <  3:696992/118‑1 (MQ=255)                              aaTATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTc                                                                                                                                           >  3:701038/1‑123 (MQ=255)                              aaTATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTc                                                                                                                                           <  4:701038/123‑1 (MQ=255)                              aaTATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCtttt                                                                                                                              >  1:267364/1‑136 (MQ=255)                              aaTATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCtttt                                                                                                                              <  2:267364/136‑1 (MQ=255)                                  ttAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCtt                                                                                                                                >  3:526195/1‑130 (MQ=255)                                  ttAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCtt                                                                                                                                <  4:526195/130‑1 (MQ=255)                                   tAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTAc                                                                                                                            <  1:196641/133‑1 (MQ=255)                                   tAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTAc                                                                                                                            >  2:196641/1‑133 (MQ=255)                                          tAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCAtttt                                                                                                                                      <  2:997852/116‑1 (MQ=255)                                          tAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCAtttt                                                                                                                                      >  1:997852/1‑116 (MQ=255)                                          tAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTc                                                                                                      >  3:841057/1‑148 (MQ=255)                                          tAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTc                                                                                                      >  3:840644/1‑148 (MQ=255)                                           aGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCAtttt                                                                                                                                      <  1:643260/115‑1 (MQ=255)                                           aGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCAtttt                                                                                                                                      >  2:643260/1‑115 (MQ=255)                                           aGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTcg                                                                                                     >  4:495780/1‑148 (MQ=255)                                                 tttGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGt                                                                                                       >  4:210913/1‑140 (MQ=255)                                                 tttGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGt                                                                                                       <  3:210913/140‑1 (MQ=255)                                                 tttGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAg                                                                                               <  4:841057/148‑1 (MQ=255)                                                     ccAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTaa                                                                                           <  4:840644/148‑1 (MQ=255)                                                                 tAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCAt                                                                                                                                         >  3:37829/1‑90 (MQ=255)                                                                 tAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCAt                                                                                                                                         <  4:37829/90‑1 (MQ=255)                                                                 tAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGc                                                                               >  4:845136/1‑148 (MQ=255)                                                                   aaaGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCATCGCGCGAGTTAATGGTTTACCTGCgg                                                                             <  3:845136/148‑1 (MQ=255)                                                                      gCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTc                                                                                                                                           <  2:641055/83‑1 (MQ=255)                                                                      gCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTc                                                                                                                                           <  4:555123/83‑1 (MQ=255)                                                                      gCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTc                                                                                                                                           >  3:555123/1‑83 (MQ=255)                                                                      gCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTc                                                                                                                                           >  1:2548/1‑83 (MQ=255)                                                                      gCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTc                                                                                                                                           >  1:641055/1‑83 (MQ=255)                                                                      gCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTc                                                                                                                                           <  2:2548/83‑1 (MQ=255)                                                                         tCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAGACCTGTGCCACAGGTTCATTTTGtt                                                                                                                                   <  1:518329/88‑1 (MQ=255)                                                                         tCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGtt                                                                                                                                   >  2:518329/1‑88 (MQ=255)                                                                          cGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCAt                                                                                                                                         <  3:623206/81‑1 (MQ=255)                                                                          cGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCAt                                                                                                                                         >  4:623206/1‑81 (MQ=255)                                                                            aGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCg                                                                    >  2:668523/1‑148 (MQ=255)                                                                                ccATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCttt                                                                                                                               <  1:169547/85‑1 (MQ=255)                                                                                ccATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCttt                                                                                                                               >  2:169547/1‑85 (MQ=255)                                                                                   tACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCcac                                                                                                                                                 <  2:1062940/64‑1 (MQ=255)                                                                                   tACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCcac                                                                                                                                                 >  1:1062940/1‑64 (MQ=255)                                                                                    aCGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGa                                                                                                                   <  1:165632/93‑1 (MQ=255)                                                                                    aCGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGa                                                                                                                   >  2:165632/1‑93 (MQ=255)                                                                                     cgcgTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTa                                                                                                                             <  3:439350/82‑1 (MQ=255)                                                                                     cgcgTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTa                                                                                                                             >  4:439350/1‑82 (MQ=255)                                                                                       cgTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCag                                                              <  3:640179/143‑1 (MQ=255)                                                                                       cgTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCag                                                              >  4:640179/1‑143 (MQ=255)                                                                                        gTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGtt                                                        <  1:668523/148‑1 (MQ=255)                                                                                          cGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTttg                                                      <  3:495780/148‑1 (MQ=255)                                                                                             cAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATg                                                                                                                 <  2:55760/86‑1 (MQ=255)                                                                                             cAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATg                                                                                                                 >  1:55760/1‑86 (MQ=255)                                                                                             cAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTaa                                                                                           >  1:1046385/1‑108 (MQ=255)                                                                                             cAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTaa                                                                                           <  2:1046385/108‑1 (MQ=255)                                                                                               gATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAgtg                                                          <  4:545307/139‑1 (MQ=255)                                                                                               gATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAgtg                                                          >  3:545307/1‑139 (MQ=255)                                                                                                    aCAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCcaca                                                                                                                                                >  1:1085714/1‑48 (MQ=255)                                                                                                    aCAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCcaca                                                                                                                                                <  2:1085714/48‑1 (MQ=255)                                                                                                       ggTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTttg                                                      <  1:860651/135‑1 (MQ=255)                                                                                                       ggTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTCTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTttg                                                      >  2:860651/1‑135 (MQ=255)                                                                                                         tgtCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCg                                                                                                                    >  3:647010/1‑71 (MQ=255)                                                                                                         tgtCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCg                                                                                                                    <  4:647010/71‑1 (MQ=255)                                                                                                              gggAATCATAACGCTTCCACTGTGAAACCTGTGGCACAGGtt                                                                                                                                            >  4:279430/1‑42 (MQ=255)                                                                                                              gggAATCATAACGCTTCCACTGTGAAACCTGTGGCACAGGtt                                                                                                                                            <  3:279430/42‑1 (MQ=255)                                                                                                                       aaCGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACTGTAACCGATGCGTTGTCCGTCGCGCGAGTTaa                                                                                           <  3:970247/82‑1 (MQ=255)                                                                                                                       aaCGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTaa                                                                                           >  4:970247/1‑82 (MQ=255)                                                                                                                       aaCGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTaa                                                                                           <  2:531162/82‑1 (MQ=255)                                                                                                                       aaCGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTaa                                                                                           >  1:531162/1‑82 (MQ=255)                                                                                                                       aaCGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAAt                                                                 >  1:659159/1‑108 (MQ=255)                                                                                                                       aaCGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAAt                                                                 <  2:659159/108‑1 (MQ=255)                                                                                                                       aaCGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGc                         >  3:599159/1‑148 (MQ=255)                                                                                                                        aCGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGcc                        >  4:157303/1‑148 (MQ=255)                                                                                                                           cTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTa                                            >  1:780673/1‑125 (MQ=255)                                                                                                                           cTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTa                                            <  2:780673/125‑1 (MQ=255)                                                                                                                           cTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCa                       >  4:1108165/1‑146 (MQ=255)                                                                                                                           cTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCa                       <  3:1108165/146‑1 (MQ=255)                                                                                                                             tCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTttgatt                                                   >  2:231624/1‑116 (MQ=255)                                                                                                                             tCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTttgatt                                                   <  1:231624/116‑1 (MQ=255)                                                                                                                                aCTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGttt                >  4:630132/1‑148 (MQ=255)                                                                                                                                 cTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAAt                                                                        >  4:205668/1‑91 (MQ=255)                                                                                                                                 cTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAAt                                                                        <  3:205668/91‑1 (MQ=255)                                                                                                                                  tgtgAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGt              <  4:599159/148‑1 (MQ=255)                                                                                                                                   gtgAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAc                                      <  2:869201/123‑1 (MQ=255)                                                                                                                                   gtgAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAc                                      >  1:869201/1‑123 (MQ=255)                                                                                                                                    atgaaCCTGTGGCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAAt                                                                 <  2:431248/92‑1 (MQ=255)                                                                                                                                    atgaaCCTGTGGCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAAt                                                                 >  1:431248/4‑95 (MQ=255)                                                                                                                                    atgaaCCTGTGGCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTa                                            >  1:1024913/4‑116 (MQ=255)                                                                                                                                    atgaaCCTGTGGCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTa                                            <  2:1024913/113‑1 (MQ=255)                                                                                                                                     gAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGtctt           >  1:1041758/1‑148 (MQ=255)                                                                                                                                        aCCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTg        >  3:700975/1‑148 (MQ=255)                                                                                                                                          cTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCg      <  3:630132/148‑1 (MQ=255)                                                                                                                                             tgCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTcgcgcg                                                                                                 <  2:807052/54‑1 (MQ=255)                                                                                                                                             tgCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTcgcgcg                                                                                                 >  1:807052/1‑54 (MQ=255)                                                                                                                                              gCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACAt  <  4:700975/148‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                  AACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACAT  >  NC_000913/1572785‑1573074 | 
| Alignment Legend | 
|---|
| Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG | 
| Unaligned base: atcgMasked matching base:atcgAlignment gap:‑Deleted base:‑ |