Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,128,262 T→C 35.3% Q45R (CAG→CGG)  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,2620TC35.3% 70.7 / 80.2 85Q45R (CAG→CGG) metJtranscriptional repressor, S‑adenosylmethionine‑binding
Reads supporting (aligned to +/- strand):  ref base T (30/25);  new base C (16/14);  total (46/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

CGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGT  >  NC_000913/4128163‑4128361
                                                                                                   |                                                                                                   
cGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCt                                                                                                     >  1:1881237/1‑100 (MQ=255)
  cAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCCga                                                                                                   >  2:26634/1‑100 (MQ=255)
   aGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTgac                                                                                                  >  2:2221601/1‑100 (MQ=255)
   aGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTgac                                                                                                  <  1:1770575/100‑1 (MQ=255)
   aGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTgac                                                                                                  <  1:744094/100‑1 (MQ=255)
    gATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTgacg                                                                                                 <  2:1522826/100‑1 (MQ=255)
    gATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTgacg                                                                                                 >  1:1463131/1‑100 (MQ=255)
      tCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTgacgac                                                                                               >  2:1299957/1‑100 (MQ=255)
        ggCATCATCCGGCAAAGGTTGACCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGAGTGACGCAGATTGTTCACCCGACGAcgc                                                                                             >  2:1894957/1‑100 (MQ=255)
         gCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGAcgcg                                                                                            <  1:1041121/100‑1 (MQ=255)
           atcatcCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACgcgtg                                                                                          <  2:397056/100‑1 (MQ=255)
             catcCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACgcgtgcg                                                                                        >  2:339712/1‑100 (MQ=255)
              atcCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCCGACGACgcgtgcgt                                                                                       >  1:1436211/1‑100 (MQ=255)
              atcCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACgcgtgcgt                                                                                       >  2:1241303/1‑100 (MQ=255)
                 cGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCCGACGACGCGTGCGTTCa                                                                                    >  2:34872/1‑100 (MQ=255)
                 cGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCCGACGACGCGTGCGTTCa                                                                                    >  1:592122/1‑100 (MQ=255)
                   gCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATc                                                                                  <  2:858219/100‑1 (MQ=255)
                   gCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATc                                                                                  <  1:719444/100‑1 (MQ=255)
                     aaaGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCCCGCTGTTGGTAGCGTTACGCAGGTTGTTAACCTGACGACACGTGCGTTCATCgg                                                                                >  1:33105/1‑100 (MQ=255)
                      aagtttttccGGCAATGGCTTGCAGAAACGCTTCGCACAGCAGCTTGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCg                                                                                 <  1:45540/90‑1 (MQ=255)
                       aGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTg                                                                              <  1:1299963/100‑1 (MQ=255)
                          ttGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAgg                                                                           >  2:337301/1‑100 (MQ=255)
                           tGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCCGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAgg                                                                           >  2:1267614/1‑99 (MQ=255)
                            gCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTCCAACCGACGACGCGTGCGTTCATCGGTGAGGAt                                                                         <  2:1475098/100‑1 (MQ=255)
                             cccGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCCCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGAtt                                                                        >  2:12292/1‑100 (MQ=255)
                              ccGGCACAGGCATGCAGAAACGCTTCGCACAACAGCTAGCTGTTGGTAGCGTGACGAAGGTTCTTACCCTGAAGACCCGTGCGTTCATCGGTGa                                                                             >  2:714305/1‑94 (MQ=255)
                               cGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCCGACGACGCGTGCGTTCATCGGTGAGGAtttt                                                                      >  1:1297098/1‑100 (MQ=255)
                                ggTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCCGACGACGCGTGCGGTCATCGGTGAGGAttttt                                                                     >  1:1571187/1‑100 (MQ=255)
                                 gcaaaGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGAtttt                                                                      <  1:654586/96‑1 (MQ=255)
                                 gcaaaGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTa                                                                    <  2:2035894/98‑1 (MQ=255)
                                 gTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCCGACGACGCGTGCGTTCATCGGTGAGGATTTTTa                                                                    <  1:1378730/100‑1 (MQ=255)
                                  tAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCCGACGACGCGTGCGTTCATCGGTGAGGATTTTTa                                                                    <  1:1953686/99‑1 (MQ=255)
                                  caaaGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTaa                                                                   >  2:30396/2‑100 (MQ=255)
                                     tccgatctaGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCCGACGACGCGTGCGTTCATCGGTGAGGATTTTTAacac                                                                <  1:348330/92‑1 (MQ=255)
                                      ggCATGTAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCCGACGACGCGTGCGTTCATCGGTGAGGATTTTTAac                                                                  <  2:778402/97‑1 (MQ=255)
                                      ggCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCCGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACAcc                                                               >  1:114524/1‑100 (MQ=255)
                                       gCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCCGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACAcc                                                               <  1:1752896/99‑1 (MQ=255)
                                            cAGAAACGCTTCGCACCGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTAATCGGTGCGGATTTTTAAAACCTTaa                                                           >  2:581528/1‑98 (MQ=255)
                                             aGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCCGACGACGCGTGCGTTCATCGGTGAGGATTTTTAacacagatcgga                                                        >  2:348331/1‑92 (MQ=255)
                                              gAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCAGCCGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAagagn                                                       <  2:843716/100‑2 (MQ=255)
                                              gAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCCGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAgg                                                       <  1:1143211/100‑1 (MQ=255)
                                                 aCGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAAt                                                    <  2:1260949/100‑1 (MQ=255)
                                                  cGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCCGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATg                                                   <  1:2232222/100‑1 (MQ=255)
                                                   gCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATgg                                                  >  1:1795389/1‑100 (MQ=255)
                                                   gCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATgg                                                  <  2:1264265/100‑1 (MQ=255)
                                                   gCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCCGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATgg                                                  <  1:752573/100‑1 (MQ=255)
                                                    cTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGa                                                 >  1:924469/1‑100 (MQ=255)
                                                    cTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGa                                                 <  2:155679/100‑1 (MQ=255)
                                                     ttCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGaa                                                >  1:461478/1‑100 (MQ=255)
                                                       cGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAAc                                              <  2:1494106/100‑1 (MQ=255)
                                                         cacaGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCCGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCg                                            <  1:2171797/100‑1 (MQ=255)
                                                         cacaGCAGCTCGCTGTTGGTAGCGTAACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCg                                            >  2:594758/1‑100 (MQ=255)
                                                          acaGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGt                                           >  1:2121067/1‑100 (MQ=255)
                                                           cagcagCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTa                                          >  2:1490849/1‑100 (MQ=255)
                                                              cagcTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACTTGACGACGCGGGCGTGCATCGGTGAGGATTTTTAACACCTTAAGTGGAATGGAAACCGTaa                                         <  2:980051/98‑1 (MQ=255)
                                                                gcTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCtt                                     >  1:641212/1‑100 (MQ=255)
                                                                gcTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGAATTTTAACACCTTAAGAGGAATGGAAACCGTAATCtt                                     >  1:1573223/1‑100 (MQ=255)
                                                                gcTCGCTGTTAGTCGCGTGACGCAGGTTGTTCACCTGACCCCGCGTGCGATCATCGGTGAGTATACTTAACACATTACGAGCAATGGAAATCGCAACCtt                                     >  2:1883243/1‑100 (MQ=255)
                                                                 cTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCttt                                    >  2:326129/1‑100 (MQ=255)
                                                                   cGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCtttt                                   >  1:524411/1‑99 (MQ=255)
                                                                    gCTGTTGGTAGCGTGACGCAGGTTGTTCACCCGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTg                                 <  1:1087460/100‑1 (MQ=255)
                                                                     cTGTGGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGa                                <  1:866646/100‑1 (MQ=255)
                                                                       gTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACt                              >  1:1254442/1‑100 (MQ=255)
                                                                       gTTGGTAGCGTGACGCAGGTTGTTCACCCGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACCCCTTAAGAGGAATGGAAACCGTAATCTTTTTGa                                >  1:1959872/1‑98 (MQ=255)
                                                                          ggTAGCGTGACGCAGGTTGTTCACCCGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGt                           >  1:1868203/1‑100 (MQ=255)
                                                                            tAGCGTGACGCAGGTTTTTCACCCGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAAAAATGGAAACCGTAATCTTTTTGACTTGTTc                         <  2:90031/100‑1 (MQ=255)
                                                                              gCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCAc                       <  1:1508442/100‑1 (MQ=255)
                                                                              gCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCAc                       <  2:459534/100‑1 (MQ=255)
                                                                              gCGTGACGCAGGTTGTTCACCCGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCAc                       >  1:1705319/1‑100 (MQ=255)
                                                                              gCGTGACGCAGGTTGTTCACCCGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCAc                       >  1:429003/1‑100 (MQ=255)
                                                                              gCGTGACGCAGGTTGTTCACCCGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCAc                       >  1:428935/1‑100 (MQ=255)
                                                                                gTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCActc                     >  1:1758944/1‑100 (MQ=255)
                                                                                gTGACGCAGGTTGTTCACCCGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCActc                     >  2:82533/1‑100 (MQ=255)
                                                                                 tGACGCAGGTTGTTCACCCGACGACGCGTGCGTTCATCGGTGAGGATTTTTAATCCCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCActct                    <  2:1089034/100‑1 (MQ=255)
                                                                                    cGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACtcttct                 >  2:1956686/1‑100 (MQ=255)
                                                                                    cGCAGGTTGTTCACCCGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACtcttct                 >  1:1718482/1‑100 (MQ=255)
                                                                                     gCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTcttctt                >  1:1820279/1‑100 (MQ=255)
                                                                                       aGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAGACCGTAATCTTTTTGACTTGTTCACTCTTCTTGc              >  2:884263/1‑100 (MQ=255)
                                                                                       aGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGc              <  2:1544838/100‑1 (MQ=255)
                                                                                        ggTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGcc             <  1:732094/100‑1 (MQ=255)
                                                                                        ggTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGcc             >  1:950146/1‑100 (MQ=255)
                                                                                         gttgttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCg            <  2:1683297/100‑1 (MQ=255)
                                                                                           tgttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTg          <  2:282959/100‑1 (MQ=255)
                                                                                            gttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTg          >  1:1962134/1‑99 (MQ=255)
                                                                                              tCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGc         <  2:530400/98‑1 (MQ=255)
                                                                                                  cTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCg   <  1:1707937/100‑1 (MQ=255)
                                                                                                   tGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGt  >  2:1879976/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGT  >  NC_000913/4128163‑4128361

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