Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,585,687 A→G Y359H (TAT→CAT)  hsdR ← type I restriction enzyme EcoKI endonuclease component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,585,6870AG99.0% 331.4 / NA 96Y359H (TAT→CAT) hsdRtype I restriction enzyme EcoKI endonuclease component
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (42/53);  minor base C (1/0);  total (43/53)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTT  >  NC_000913/4585550‑4585822
                                                                                                                                         |                                                                                                                                       
gCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGct                                                                                                                                       >  1:1312560/1‑140 (MQ=255)
 cGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGctc                                                                                                                                      <  1:1075355/140‑1 (MQ=255)
 cGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGctc                                                                                                                                      <  2:975007/140‑1 (MQ=255)
    gTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTAGTGTCCTGCCAGCTGGTTTCATGCTCTgg                                                                                                                                   >  2:1275052/1‑140 (MQ=255)
     tATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGc                                                                                                                                  >  1:1009348/1‑140 (MQ=255)
      aTCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCa                                                                                                                                 <  1:1279517/140‑1 (MQ=255)
      aTCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGCTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCa                                                                                                                                 <  1:790938/140‑1 (MQ=255)
       tCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCAc                                                                                                                                >  1:1478625/1‑140 (MQ=255)
       tCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCAc                                                                                                                                >  1:1445601/1‑140 (MQ=255)
        cacgcacgTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCa                                                                                                                                 <  1:1265749/138‑1 (MQ=255)
        cacgcacgTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCa                                                                                                                                 >  2:1265749/1‑138 (MQ=255)
        cacgcacgTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACt                                                                                                                               <  1:1320491/140‑1 (MQ=255)
        cacgcacgTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACt                                                                                                                               <  2:1445601/140‑1 (MQ=255)
              cgTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACt                                                                                                                               >  1:269412/1‑134 (MQ=255)
              cgTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACt                                                                                                                               <  2:269412/134‑1 (MQ=255)
               gTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCGGCTAGTTTCATGCACTGGCACTGCTTCAt                                                                                                                        >  1:644559/1‑140 (MQ=255)
               gTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCAt                                                                                                                        >  1:810829/1‑140 (MQ=255)
               gTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCAt                                                                                                                        >  1:1010047/1‑140 (MQ=255)
               gTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTCTGCCGCTGGTGTACCGGCAGCTGGTCTCATCCTCTGGCACTGCTTCAt                                                                                                                        >  2:713841/1‑140 (MQ=255)
                 cgcgATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGc                                                                                                                      >  2:49259/1‑140 (MQ=255)
                 cgcgATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGc                                                                                                                      >  1:313910/1‑140 (MQ=255)
                 cgcgATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGc                                                                                                                      <  2:988205/140‑1 (MQ=255)
                 cgcgATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGc                                                                                                                      >  1:107732/1‑140 (MQ=255)
                        ccAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCAt                                                                                                               <  1:1230804/140‑1 (MQ=255)
                         cAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGGTGGTTTCATGCGCTGGCACTGATTCATGCACTTCATc                                                                                                              >  1:482908/1‑140 (MQ=255)
                           gATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCg                                                                                                            >  1:89323/1‑140 (MQ=255)
                           gATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCg                                                                                                            >  1:945220/1‑140 (MQ=255)
                           gATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCg                                                                                                            >  1:677261/1‑140 (MQ=255)
                           gATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGCCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGCCCTGCCAGCTGGGTGCCGGCTCTGGAACTGCTTCATGCAGTTCATCCg                                                                                                            >  2:743105/1‑140 (MQ=255)
                             tGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGt                                                                                                          <  1:156945/140‑1 (MQ=255)
                             tGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGt                                                                                                          >  1:176917/1‑140 (MQ=255)
                             tGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGt                                                                                                          <  2:1205915/140‑1 (MQ=255)
                               ccGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGa                                                                                                        >  2:976943/1‑140 (MQ=255)
                                cGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTGATGCTCTGGCAGTGGTTCATGGACTTCATCCGGTGag                                                                                                       >  1:396211/1‑140 (MQ=255)
                                cGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGag                                                                                                       <  2:107732/140‑1 (MQ=255)
                                cGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGACAGCTGGTTTCAAGCTATGGCACTGCTTCACGCACTTCATCCGGTGag                                                                                                       >  2:640846/1‑140 (MQ=255)
                                  cTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGagta                                                                                                     <  2:313910/140‑1 (MQ=255)
                                  cTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGagta                                                                                                     >  2:1274180/1‑140 (MQ=255)
                                  cTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTAGCACTGCTGCATGCACTTCATCCGGTGagta                                                                                                     >  2:462105/1‑140 (MQ=255)
                                       ggTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCGGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCt                                                                                                >  1:639827/1‑140 (MQ=255)
                                        gTCTTCATTGTTGCGCGGTATTGGCCCCCATTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTc                                                                                               <  1:556536/140‑1 (MQ=255)
                                        gTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTc                                                                                               <  1:900219/140‑1 (MQ=255)
                                             cATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTcaagca                                                                                          <  2:1009348/140‑1 (MQ=255)
                                                   tgCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTAGCACTGCTGCATCCACTTCATCCGGAGAGTCGTGCTCAAGCAAGGttt                                                                                    >  2:556536/1‑140 (MQ=255)
                                                     cgcgGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTcc                                                                                  <  1:158448/140‑1 (MQ=255)
                                                       cgGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTccct                                                                                >  1:81697/1‑139 (MQ=255)
                                                       cgGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCg                                                                                >  1:1407201/1‑140 (MQ=255)
                                                          tATTCGCGCCCGTTGGTTGAGTCGCAGAAGGCGATCTTACACCGTTGTCTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAg                                                                             <  2:397515/140‑1 (MQ=255)
                                                          tATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAg                                                                             <  2:1312560/140‑1 (MQ=255)
                                                          tATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAg                                                                             >  2:1215651/1‑140 (MQ=255)
                                                          tATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAg                                                                             <  1:271493/140‑1 (MQ=255)
                                                                 gcCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATACACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCa                                                                      <  1:1428525/140‑1 (MQ=255)
                                                                 gcCCGTTGGTCGAGTAGCAGAAGGGCATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCa                                                                      <  1:653880/140‑1 (MQ=255)
                                                                    cGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGGGTAGTGCGCAAGCAGGGTTTCTCGCAGGAAGCCATTa                                                                   >  1:397515/1‑140 (MQ=255)
                                                                    cGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTa                                                                   >  2:824361/1‑140 (MQ=255)
                                                                    cGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTa                                                                   >  2:175088/1‑140 (MQ=255)
                                                                    cGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTa                                                                   >  2:1124106/1‑140 (MQ=255)
                                                                      ttGGTCGAGTAGCAGTAGGGGAATTGAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATc                                                                 <  1:640846/140‑1 (MQ=255)
                                                                      ttGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATc                                                                 >  1:286982/1‑140 (MQ=255)
                                                                      ttGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATc                                                                 >  2:140847/1‑140 (MQ=255)
                                                                      ttGGTCGAGTAGCAGAACGGGAATTTAAACCCTTGCTTGCCGCTGGTGCCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATc                                                                 <  2:396211/140‑1 (MQ=255)
                                                                      ttGGTCGAGTAGCAGAAAGGGAGTTTAAACCGTTGGTTGCCGCTGGTGCCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATc                                                                 <  2:644559/140‑1 (MQ=255)
                                                                       tGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCg                                                                <  1:54934/140‑1 (MQ=255)
                                                                         gTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGaa                                                              <  2:273472/140‑1 (MQ=255)
                                                                           cGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAAc                                                            >  1:995932/1‑140 (MQ=255)
                                                                             aGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACAt                                                          >  2:1406599/1‑140 (MQ=255)
                                                                             aGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACAt                                                          <  1:824361/140‑1 (MQ=255)
                                                                              gTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACAtt                                                         <  2:810829/140‑1 (MQ=255)
                                                                              gTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACAtt                                                         <  2:1010047/140‑1 (MQ=255)
                                                                              gTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACAtt                                                         <  1:976943/140‑1 (MQ=255)
                                                                                aGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTa                                                       <  1:1124106/140‑1 (MQ=255)
                                                                                  cAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACt                                                     >  1:810417/1‑140 (MQ=255)
                                                                                   aGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTg                                                    <  1:135213/140‑1 (MQ=255)
                                                                                    gAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCATTCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGt                                                   <  2:81697/140‑1 (MQ=255)
                                                                                    gAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGt                                                   <  2:677261/140‑1 (MQ=255)
                                                                                       ggggATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGc                                                >  1:182155/1‑140 (MQ=255)
                                                                                          gATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGt                                                   >  1:600436/1‑134 (MQ=255)
                                                                                          gATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGt                                                   <  2:600436/134‑1 (MQ=255)
                                                                                          gATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGAt                                             <  1:140847/140‑1 (MQ=255)
                                                                                            ttttAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATAc                                           >  1:917626/1‑140 (MQ=255)
                                                                                              ttAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGa                                         <  2:846773/140‑1 (MQ=255)
                                                                                              ttAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGa                                         <  2:1478625/140‑1 (MQ=255)
                                                                                                  aCCGTTGTTTGCCGCTAGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCGATTATCGAAACATTTACTGTAGCGATACGACTCa                                     <  1:1275052/140‑1 (MQ=255)
                                                                                                    cGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCAtt                                   <  1:856459/140‑1 (MQ=255)
                                                                                                       tgttTGCCGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGCAACTCCTTCAGCCATTTCTTCCGGGGGTAAGGGCTCAAGAAGGGTTTCCCCCGGGAAGCCTTTTCCAAAACTTTTATTTAACGGTTCGACTCATTGAGc                               >  2:78370/1‑140 (MQ=255)
                                                                                                        gttTGCCGCTGGTGTACTGCCAGCTAGTTCCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGc                               <  1:713841/140‑1 (MQ=255)
                                                                                                         gagacagCTGGTGTCCTGCCAGATGGTTTCATGCTCTCGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGcc                              <  1:495383/134‑1 (MQ=255)
                                                                                                              cGCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCg                         <  2:176917/140‑1 (MQ=255)
                                                                                                               gCTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGAACGACTCATTGAGCCTGGCggg                       <  2:89323/140‑1 (MQ=255)
                                                                                                                cTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCggg                       <  2:1407201/140‑1 (MQ=255)
                                                                                                                cTGGTGTCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCggg                       <  1:1215651/140‑1 (MQ=255)
                                                                                                                  ggTGTCCTGCCAGCTGGTTTCATGCTCTGGCTCTGCTTCTTGCACTTCCTCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGaa                     <  2:482908/140‑1 (MQ=255)
                                                                                                                   gtgtCCTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAAc                    >  1:90457/1‑140 (MQ=255)
                                                                                                                      tCCTGCCAGCTGGTTTCATGCTGTGTCACTGCTTCATCCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTc                 <  1:743105/140‑1 (MQ=255)
                                                                                                                        cTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCTCTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGa               <  1:1406599/140‑1 (MQ=255)
                                                                                                                        cTGCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGa               <  2:182155/140‑1 (MQ=255)
                                                                                                                          gCCAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGata             <  1:49259/140‑1 (MQ=255)
                                                                                                                           ccAGCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAAGATTTACTGGAGCCATACGACTCATTAAGCCCGGCGGGAACgctgtctc            >  2:495383/1‑133 (MQ=255)
                                                                                                                              gCTGGTTTCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGt         <  1:175088/140‑1 (MQ=255)
                                                                                                                                     tCATGCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGttt  <  2:90457/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                       
GCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTT  >  NC_000913/4585550‑4585822

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