Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,128,197 T→C 100% T67A (ACC→GCC)  metJ ← transcriptional repressor, S‑adenosylmethionine‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,128,1970TC89.2% 277.9 / 17.0 100T67A (ACC→GCC) metJtranscriptional repressor, S‑adenosylmethionine‑binding
Reads supporting (aligned to +/- strand):  ref base T (2/9);  new base C (37/52);  total (39/61)
Fisher's exact test for biased strand distribution p-value = 1.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

TAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTT  >  NC_000913/4128101‑4128293
                                                                                                |                                                                                                
tAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCaaa                                                                                               <  2:707759/100‑1 (MQ=255)
 aaTCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAg                                                                                              <  1:1809943/100‑1 (MQ=255)
   tCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGc                                                                                            >  2:436446/1‑100 (MQ=255)
    ccccATCTCACGCATGATCTCTTTTGCCGCTTTCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAg                                                                                              <  1:94074/97‑1 (MQ=255)
    ccccATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCa                                                                                           <  2:1208019/100‑1 (MQ=255)
    ccccATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCa                                                                                           >  1:2827483/1‑100 (MQ=255)
    ccccATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCa                                                                                           >  1:1947676/1‑100 (MQ=255)
     ctcATCTCACGCATTATCTCTTTTGCCGCTTCCGGTATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCa                                                                                           <  1:1169670/97‑1 (MQ=255)
     cccATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCAt                                                                                          >  2:2342165/1‑100 (MQ=255)
     cccATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCAt                                                                                          >  1:392438/1‑100 (MQ=255)
      ccATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATg                                                                                         <  2:65555/100‑1 (MQ=255)
      ccATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATg                                                                                         <  1:451796/100‑1 (MQ=255)
       cATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGc                                                                                        <  2:480601/100‑1 (MQ=255)
        aTCTCTCGCATGATCTCTTTTGCCGCTTCCGGTATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCa                                                                                       <  1:403270/100‑1 (MQ=255)
         tctcACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAg                                                                                      <  2:1600397/100‑1 (MQ=255)
         tctcACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAg                                                                                      >  1:551828/1‑100 (MQ=255)
            cACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGaaa                                                                                   <  1:2342150/100‑1 (MQ=255)
            cACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGaaa                                                                                   >  2:2731871/1‑100 (MQ=255)
             aCGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGTTCGGCATCATCCGGCAAAGGTTACCCGGCAAAGGCATGCAGcaa                                                                                   >  1:1751111/1‑99 (MQ=255)
             aCGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAAc                                                                                  <  2:2556635/100‑1 (MQ=255)
              cGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACg                                                                                 >  2:461885/1‑100 (MQ=255)
                  tGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCtt                                                                              <  2:886500/99‑1 (MQ=255)
                  tGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTc                                                                             >  1:1218360/1‑100 (MQ=255)
                   gATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCg                                                                            >  2:2022764/1‑100 (MQ=255)
                       tctTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGcaca                                                                        >  1:2100624/1‑100 (MQ=255)
                       tctTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGcaca                                                                        <  1:310125/100‑1 (MQ=255)
                        ctTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCAcag                                                                       >  1:2459752/1‑100 (MQ=255)
                         ttttGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCAcagc                                                                      <  2:2443401/100‑1 (MQ=255)
                         ttttGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCAcagc                                                                      <  1:361007/100‑1 (MQ=255)
                         ttttGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCAcagc                                                                      <  2:2580652/100‑1 (MQ=255)
                         ttttGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCAcagc                                                                      <  2:553115/100‑1 (MQ=255)
                         ttttGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCAcagc                                                                      <  1:996564/100‑1 (MQ=255)
                          tttGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCAcagca                                                                     <  1:2769138/100‑1 (MQ=255)
                             gCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCt                                                                  <  1:2821413/100‑1 (MQ=255)
                             gCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCt                                                                  >  2:1562976/1‑100 (MQ=255)
                              ccGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTc                                                                 >  2:130168/1‑100 (MQ=255)
                              ccGCTTCCGGGATATCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTc                                                                 >  1:1936837/1‑100 (MQ=255)
                               cGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCg                                                                >  2:4303/1‑100 (MQ=255)
                                  ttCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTg                                                             >  2:1785234/1‑100 (MQ=255)
                                   tCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCCTCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGt                                                            >  1:1791099/1‑100 (MQ=255)
                                   tCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGGAGAAACGCTTCGCACAGCAGCTCGCTGt                                                            <  2:920940/100‑1 (MQ=255)
                                    ccGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGtt                                                           <  2:1662766/100‑1 (MQ=255)
                                    ccGGGATTTCGTCGCTGCACTCTTTAAGCAGATCGGCCCAATCCGGCAAAGGTTGCACGGTAAAGACATGCAGACACGCTTCGCACAGCAGCTCGCTGtt                                                           <  2:1307923/100‑1 (MQ=255)
                                     cGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTg                                                          >  1:1676693/1‑100 (MQ=255)
                                     cGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTg                                                          >  1:854065/1‑100 (MQ=255)
                                     cGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTg                                                          <  2:2354624/100‑1 (MQ=255)
                                     cGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTg                                                          >  2:2848808/1‑100 (MQ=255)
                                     cGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTg                                                          <  2:656036/100‑1 (MQ=255)
                                      gggATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTgg                                                         <  2:1546907/100‑1 (MQ=255)
                                      gggATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTgg                                                         <  1:348307/100‑1 (MQ=255)
                                      gggATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTgg                                                         <  2:1413848/100‑1 (MQ=255)
                                      gggATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTgg                                                         <  2:834043/100‑1 (MQ=255)
                                       ggATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGt                                                        >  2:1189195/1‑100 (MQ=255)
                                        gATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTgg                                                         <  2:2628656/98‑1 (MQ=255)
                                         aTTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAg                                                      <  1:79442/100‑1 (MQ=255)
                                         aTTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTACCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAg                                                      >  1:2887459/1‑100 (MQ=255)
                                          nttCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGCTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGc                                                     >  1:836356/2‑100 (MQ=255)
                                           ttCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCg                                                    >  2:2204274/1‑100 (MQ=255)
                                             cgtcgCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTg                                                  <  1:478419/100‑1 (MQ=255)
                                             cgtcgCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGTCAGGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTg                                                  >  2:131704/1‑100 (MQ=255)
                                              gtcgCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGa                                                 >  2:1822050/1‑100 (MQ=255)
                                               tcgCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGAc                                                <  1:769058/100‑1 (MQ=255)
                                               tcgCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGAc                                                <  1:638626/100‑1 (MQ=255)
                                                 gCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGc                                              <  1:2587841/100‑1 (MQ=255)
                                                 gCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGc                                              >  2:2284387/1‑100 (MQ=255)
                                                     cgcTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGgt                                          <  1:1360630/100‑1 (MQ=255)
                                                     cgcTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGgt                                          <  1:1306867/100‑1 (MQ=255)
                                                     cgcTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGgt                                          <  1:1085945/100‑1 (MQ=255)
                                                     cgcTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGgt                                          <  1:96507/100‑1 (MQ=255)
                                                     cgcTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGgt                                          <  2:207072/100‑1 (MQ=255)
                                                       ctctTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGgttg                                        >  2:1028850/1‑100 (MQ=255)
                                                        tctTTCCGCAGATCCGCATCATCCGGCAAAGGTTGCCCCGCAAAGGCATGCCGAAACGCTTCGCATAGCAGCTCGCTGTTGGTAGCGTGACGCAGgttgt                                       <  1:1828063/100‑1 (MQ=255)
                                                        tctTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGgttgt                                       <  2:730279/100‑1 (MQ=255)
                                                         ctTTACGCAGATCGGCATCATCCGGCAAAGGTTGTCCGGTAAAGGCATGCAGAAACGCTTTGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGgttgt                                       <  1:1791598/99‑1 (MQ=255)
                                                             aCGCAGATCGTCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCAcc                                  <  2:2284740/100‑1 (MQ=255)
                                                                 aGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTgac                              <  1:103435/100‑1 (MQ=255)
                                                                 aGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTgac                              <  2:1319802/100‑1 (MQ=255)
                                                                  gATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTgacg                             <  2:234103/100‑1 (MQ=255)
                                                                   aTCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTgacga                            >  2:1553546/1‑100 (MQ=255)
                                                                     cGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGacgacg                          >  2:1708323/1‑100 (MQ=255)
                                                                       gCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGAcgcg                        >  1:1350315/1‑100 (MQ=255)
                                                                       gCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGAcgcg                        <  2:517874/100‑1 (MQ=255)
                                                                         atcatcCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACgcgtg                      <  1:2436432/100‑1 (MQ=255)
                                                                          tcatcCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACgcgtgc                     >  2:1634087/1‑100 (MQ=255)
                                                                           catcCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACgcgtgcg                    <  2:1640127/100‑1 (MQ=255)
                                                                           catcCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACgcgtgcg                    >  2:1393491/1‑100 (MQ=255)
                                                                             tcCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGtt                  <  2:836616/100‑1 (MQ=255)
                                                                              ccGGCAAAGGTTGCGCGGCAAAGGCATGCGGAAACGCTTCGCACAGCAGCTCGCGGTTGGTAGCGTGACGCAGGTAGTTTACCTGGCGACgcgtg                      >  1:2314779/1‑95 (MQ=255)
                                                                                 gCAAATGTGGCCGGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCTGGTGGTTTACCTGACGACGCGTGCGTTCATc              <  1:1358323/100‑1 (MQ=255)
                                                                                 gCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATc              <  2:1190468/100‑1 (MQ=255)
                                                                                   aaaGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCgg            <  2:497120/100‑1 (MQ=255)
                                                                                   aaaGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCgg            <  1:2452730/100‑1 (MQ=255)
                                                                                      ggTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTTACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGt           <  2:258594/98‑1 (MQ=255)
                                                                                       gTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAg        <  2:2470674/100‑1 (MQ=255)
                                                                                       gTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAg        <  1:522662/100‑1 (MQ=255)
                                                                                         tGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGa      <  2:1354175/100‑1 (MQ=255)
                                                                                          gCCGGGCAAAGGCATGCAGAAACGAGTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGGGCGTTCATCTGTGAGGa      <  1:728769/99‑1 (MQ=255)
                                                                                          gCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGa      >  2:1858085/1‑99 (MQ=255)
                                                                                           cccGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGAtt    >  2:2346111/1‑100 (MQ=255)
                                                                                             cGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGAtttt  >  2:521512/1‑100 (MQ=255)
                                                                                                |                                                                                                
TAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTT  >  NC_000913/4128101‑4128293

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