Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,427,556 A→G T39A (ACG→GCG)  cycA → transport of D‑alanine, D‑serine, and glycine

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,427,5560AG99.6% 770.3 / NA 241T39A (ACG→GCG) cycAtransport of D‑alanine, D‑serine, and glycine
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (100/140);  minor base C (0/1);  total (100/141)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 248 total aligned reads displayed.

GCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTT  >  NC_000913/4427466‑4427645
                                                                                          |                                                                                         
    atgatCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGgt                                                                                     >  2:4567746/1‑93 (MQ=255)
     tgatCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGgtt                                                                                    >  1:4223008/1‑93 (MQ=255)
     tgatCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGgtt                                                                                    >  1:2578826/1‑93 (MQ=255)
      gatCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGgttg                                                                                   >  1:4954135/1‑93 (MQ=255)
       atCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGgttgt                                                                                  >  2:3611554/1‑93 (MQ=255)
       atCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGgttgt                                                                                  >  1:810551/1‑93 (MQ=255)
        tCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGgttgtt                                                                                 <  1:4530259/93‑1 (MQ=255)
          aGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTa                                                                               <  1:3871229/93‑1 (MQ=255)
          aGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTa                                                                               <  2:3243570/93‑1 (MQ=255)
           ggCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTAt                                                                              <  1:475996/93‑1 (MQ=255)
             cTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTTGTGCGGGGTTGTTTATgg                                                                            >  2:3756734/1‑93 (MQ=255)
                 ggCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTc                                                                        <  1:4311869/93‑1 (MQ=255)
                   cTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTg                                                                      <  1:3706585/93‑1 (MQ=255)
                    tGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTgg                                                                     <  1:113364/93‑1 (MQ=255)
                     gAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGc                                                                    >  2:2259706/1‑93 (MQ=255)
                     gAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGc                                                                    >  2:6303391/1‑93 (MQ=255)
                      aaCAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCa                                                                   >  2:378774/1‑93 (MQ=255)
                       aCAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCaa                                                                  >  2:5498278/1‑93 (MQ=255)
                       aCAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCaa                                                                  >  2:5605229/1‑93 (MQ=255)
                            cGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGa                                                             <  1:5114001/93‑1 (MQ=255)
                            cGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGa                                                             <  2:1386852/93‑1 (MQ=255)
                            cGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGa                                                             <  1:4036063/93‑1 (MQ=255)
                              cTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGAtt                                                           <  1:693281/93‑1 (MQ=255)
                              cTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGGTGTTTATGGGGTCTGGCAAAACGAtt                                                           >  2:127020/1‑93 (MQ=255)
                               tACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTa                                                          >  2:262663/1‑93 (MQ=255)
                                aCGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGGAAAACGATTAg                                                         >  2:2984659/1‑93 (MQ=255)
                                 cGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGc                                                        <  2:2044097/93‑1 (MQ=255)
                                  ggCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGcc                                                       >  2:4692878/1‑93 (MQ=255)
                                   gcgcAATCTCACAAACCGACATATTCAGCTTAGTGCCATTGGCTGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCt                                                      <  1:820496/93‑1 (MQ=255)
                                     acAATCACACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTg                                                    <  2:5065002/92‑1 (MQ=255)
                                       aaTCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGcc                                                  >  1:238980/1‑93 (MQ=255)
                                       aaTCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGcc                                                  <  2:5862545/93‑1 (MQ=255)
                                       aaTCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGcc                                                  <  1:2356611/93‑1 (MQ=255)
                                        aTCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCg                                                 >  1:2085099/1‑93 (MQ=255)
                                         tctcACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCgg                                                <  2:4638741/93‑1 (MQ=255)
                                         tctcACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCgg                                                <  1:1302490/93‑1 (MQ=255)
                                         tctcACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCgg                                                <  1:2842736/93‑1 (MQ=255)
                                         tctcACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCgg                                                >  1:3342614/1‑93 (MQ=255)
                                           tcACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGc                                              <  1:1090578/93‑1 (MQ=255)
                                           tcACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGc                                              >  2:5594035/1‑93 (MQ=255)
                                            cacaAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGcc                                             <  1:5288282/93‑1 (MQ=255)
                                             acaAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCg                                            <  1:4079103/93‑1 (MQ=255)
                                              caAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGt                                           <  2:4810099/93‑1 (MQ=255)
                                               aaaCCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTc                                          <  2:2877101/93‑1 (MQ=255)
                                               aaaCCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTc                                          <  2:5660975/93‑1 (MQ=255)
                                                aaCCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCg                                         <  2:3120301/93‑1 (MQ=255)
                                                aaCCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCg                                         <  2:1908960/93‑1 (MQ=255)
                                                 agcGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGa                                        >  2:3972973/3‑93 (MQ=255)
                                                 aCCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGa                                        <  1:5667520/93‑1 (MQ=255)
                                                  ccGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGAt                                       <  2:1478554/93‑1 (MQ=255)
                                                   cGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATc                                      >  2:5300594/1‑93 (MQ=255)
                                                   cGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATc                                      >  2:1070196/1‑93 (MQ=255)
                                                    gACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCa                                     >  1:4577148/1‑93 (MQ=255)
                                                    gACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCa                                     <  1:3247746/93‑1 (MQ=255)
                                                    gACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCa                                     >  1:1166046/1‑93 (MQ=255)
                                                    gACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCa                                     <  1:1720922/93‑1 (MQ=255)
                                                     aCATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCAt                                    >  1:2292903/1‑93 (MQ=255)
                                                      cATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCAtt                                   >  2:4311122/1‑93 (MQ=255)
                                                      cATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCAtt                                   >  2:4529065/1‑93 (MQ=255)
                                                       atatTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCAttt                                  >  1:2191630/1‑93 (MQ=255)
                                                        tatTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTTTGGGTCTTGGCAAAACAATTGCCCTGCCGGGCCGTCGGTCATTTTc                                >  1:173867/1‑93 (MQ=255)
                                                        tatTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCAtttt                                 <  1:1279449/93‑1 (MQ=255)
                                                         atTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTc                                <  1:1875349/93‑1 (MQ=255)
                                                         atTCAGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTc                                <  1:5130952/93‑1 (MQ=255)
                                                             aGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGttt                            >  2:3284098/1‑93 (MQ=255)
                                                             aGCTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGttt                            <  2:2182915/93‑1 (MQ=255)
                                                               cTTATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTtat                          <  1:5324582/93‑1 (MQ=255)
                                                                ttATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTtata                         >  2:2488139/1‑93 (MQ=255)
                                                                ttATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTtata                         <  2:736127/93‑1 (MQ=255)
                                                                ttATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTtata                         >  1:2154425/1‑93 (MQ=255)
                                                                 tATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTgtat                        <  1:264228/93‑1 (MQ=255)
                                                                 tATTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTatat                        >  1:899942/1‑93 (MQ=255)
                                                                  aTTGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATg                       <  2:4178686/93‑1 (MQ=255)
                                                                   ttGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGa                      >  1:4697077/1‑93 (MQ=255)
                                                                   ttGCCATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGa                      >  1:2486856/1‑93 (MQ=255)
                                                                     tccATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATc                    <  1:1637195/92‑1 (MQ=255)
                                                                       cATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCAt                  >  2:5202125/1‑93 (MQ=255)
                                                                       cATTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCAt                  <  1:4803439/93‑1 (MQ=255)
                                                                        aTTGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCAtt                 <  2:4892649/93‑1 (MQ=255)
                                                                        aTTGCCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCATCGATCATTTTCGTTTATATGATCAtt                 <  2:5422280/93‑1 (MQ=255)
                                                                         ttGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTg                <  1:3076791/93‑1 (MQ=255)
                                                                          tGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTgg               >  2:908702/1‑93 (MQ=255)
                                                                          tGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTgg               <  1:4277958/93‑1 (MQ=255)
                                                                          tGGCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTgg               <  1:6106156/93‑1 (MQ=255)
                                                                           ggCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGt              <  2:3451797/93‑1 (MQ=255)
                                                                           ggCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGt              <  2:2405717/93‑1 (MQ=255)
                                                                            gCGGTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGtt             >  1:5645060/1‑93 (MQ=255)
                                                                              ggTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGtttt           >  2:6137980/1‑93 (MQ=255)
                                                                               gggCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTa          <  2:1623909/91‑1 (MQ=255)
                                                                               gTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTa          <  1:2610607/93‑1 (MQ=255)
                                                                               gTGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTa          <  1:2684509/93‑1 (MQ=255)
                                                                               gTGCCATGGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTa          <  2:5546397/93‑1 (MQ=255)
                                                                                tGCCATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTAt         <  1:3509184/93‑1 (MQ=255)
                                                                                  ccATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGc       >  2:2809908/1‑93 (MQ=255)
                                                                                   cATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGct      >  1:5658153/1‑93 (MQ=255)
                                                                                   cATTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGct      >  1:3495279/1‑93 (MQ=255)
                                                                                    aTTGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGctt     >  1:1525954/1‑92 (MQ=255)
                                                                                     ttGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGctct    >  1:4369906/1‑93 (MQ=255)
                                                                                     ttGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGctct    >  2:3526654/1‑93 (MQ=255)
                                                                                      tGGTGCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCtt   <  1:2894693/93‑1 (MQ=255)
                                                                                          |                                                                                         
GCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTT  >  NC_000913/4427466‑4427645

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