Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 606,299 A→T V362E (GTG→GAG)  ybdK ← weak gamma‑glutamyl:cysteine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913606,2990AT92.3% 294.3 / 12.2 91V362E (GTG→GAG) ybdKweak gamma‑glutamyl:cysteine ligase
Reads supporting (aligned to +/- strand):  ref base A (4/3);  new base T (45/39);  total (49/42)
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.89e-01

GTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAAT  >  NC_000913/606206‑606393
                                                                                             |                                                                                              
gTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGcc                                                                                          <  2:1005466/100‑1 (MQ=255)
gTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGcc                                                                                          <  2:865414/100‑1 (MQ=255)
gTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGcc                                                                                          <  1:754100/100‑1 (MQ=255)
gTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGcc                                                                                          <  1:367516/100‑1 (MQ=255)
gTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGcc                                                                                          >  2:490128/1‑100 (MQ=255)
 ttAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGccc                                                                                         <  2:724773/100‑1 (MQ=255)
  tAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCa                                                                                        >  2:1378841/1‑100 (MQ=255)
  tAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCa                                                                                        >  1:888920/1‑100 (MQ=255)
  tAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCa                                                                                        >  2:1781753/1‑100 (MQ=255)
     aTGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATc                                                                                     <  2:1851679/100‑1 (MQ=255)
      tGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATc                                                                                     >  2:506659/1‑99 (MQ=255)
      tGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCa                                                                                    <  1:538785/100‑1 (MQ=255)
      tGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCa                                                                                    >  1:1184241/1‑100 (MQ=255)
          aGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGa                                                                                >  2:73936/1‑100 (MQ=255)
          aGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGa                                                                                <  1:22094/100‑1 (MQ=255)
          aGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGa                                                                                >  1:1882158/1‑100 (MQ=255)
          aGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGa                                                                                >  2:1204156/1‑100 (MQ=255)
           gTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAg                                                                               >  1:1302727/1‑100 (MQ=255)
           gTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAg                                                                               >  1:685867/1‑100 (MQ=255)
           gTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAg                                                                               >  1:1153946/1‑100 (MQ=255)
            ttAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAgc                                                                              <  2:1493246/100‑1 (MQ=255)
              aaaaCACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAgccg                                                                            >  1:254484/1‑100 (MQ=255)
              aaaaCACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAgccg                                                                            >  2:1716383/1‑100 (MQ=255)
               aaaCACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAgccgc                                                                           >  1:547525/1‑100 (MQ=255)
               aaaCACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAgccgc                                                                           <  1:1764440/100‑1 (MQ=255)
                aaCACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAgccgcc                                                                          >  2:257948/1‑100 (MQ=255)
                   acCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATc                                                                       <  1:1908609/100‑1 (MQ=255)
                   acCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATc                                                                       <  1:1942468/100‑1 (MQ=255)
                   acCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATc                                                                       >  2:1467992/1‑100 (MQ=255)
                    ccATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCg                                                                      >  1:1286083/1‑100 (MQ=255)
                        tGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGcga                                                                  >  2:1886260/1‑100 (MQ=255)
                         gTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGcgac                                                                 <  1:740105/100‑1 (MQ=255)
                           cAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGcgacga                                                               >  1:1818996/1‑100 (MQ=255)
                                ttccgatcttGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATcg                                                          <  1:387200/91‑1 (MQ=255)
                                atatTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATAGGCGCCGAAATcg                                                          >  1:885872/1‑100 (MQ=255)
                                  atTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATcgcg                                                        <  2:281738/100‑1 (MQ=255)
                                  atTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATcgcg                                                        <  2:996260/100‑1 (MQ=255)
                                   ttcttcTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCgcgc                                                       <  1:1204162/100‑1 (MQ=255)
                                    tcttctGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCa                                                      <  1:1022888/100‑1 (MQ=255)
                                    tcttctGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCGTACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCa                                                      >  1:173398/1‑100 (MQ=255)
                                       tctGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTa                                                   >  1:292194/1‑100 (MQ=255)
                                         tGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGc                                                 <  1:780042/100‑1 (MQ=255)
                                         tGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATcgagatcggaa                                                 >  2:387199/1‑91 (MQ=255)
                                          gTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCt                                                >  1:596896/1‑100 (MQ=255)
                                          gTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCt                                                >  2:1962649/1‑100 (MQ=255)
                                           tAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTg                                               >  1:1682684/1‑100 (MQ=255)
                                            aaGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTgc                                              >  1:1321110/1‑100 (MQ=255)
                                            aaGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTgc                                              <  1:490129/100‑1 (MQ=255)
                                              ggCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCCGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTgcgt                                            >  1:61410/1‑99 (MQ=255)
                                              ggCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTgcgc                                            >  1:855427/1‑100 (MQ=255)
                                               gCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCt                                           <  1:65954/100‑1 (MQ=255)
                                                cAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCtt                                          <  1:1276399/100‑1 (MQ=255)
                                                  aTTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCttca                                        <  1:296897/100‑1 (MQ=255)
                                                    tCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCttcatt                                      >  1:1802625/1‑100 (MQ=255)
                                                    tCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCttcatt                                      <  2:1952909/100‑1 (MQ=255)
                                                     cccGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATTAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCttcattc                                     <  2:173397/100‑1 (MQ=255)
                                                      ccGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCttcattca                                    <  2:579032/100‑1 (MQ=255)
                                                         gCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGAc                                 <  1:1781768/100‑1 (MQ=255)
                                                         gCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGAc                                 <  2:1818980/100‑1 (MQ=255)
                                                            tAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGc                              <  2:1286077/100‑1 (MQ=255)
                                                              gTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTg                            >  1:1395109/1‑100 (MQ=255)
                                                                cACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTgac                          >  1:100686/1‑100 (MQ=255)
                                                                cACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTgac                          >  2:796573/1‑100 (MQ=255)
                                                                cACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTgac                          <  2:1882140/100‑1 (MQ=255)
                                                                 aCCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTgacg                         <  2:1330376/100‑1 (MQ=255)
                                                                  ccGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTgacga                        >  1:1946526/1‑100 (MQ=255)
                                                                     gCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCt                     <  1:73936/100‑1 (MQ=255)
                                                                      cccAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTACATTCAGACCGCTGACGACCTg                    <  2:1184235/100‑1 (MQ=255)
                                                                      cccAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTgacg                         <  2:1735101/95‑1 (MQ=255)
                                                                      cccAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTg                    >  1:13341/1‑100 (MQ=255)
                                                                       ccAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTgg                   >  1:687509/1‑100 (MQ=255)
                                                                         aGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCg                 >  1:1921599/1‑100 (MQ=255)
                                                                          naTCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGa                >  2:695834/2‑100 (MQ=255)
                                                                          gATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGa                <  1:1468001/100‑1 (MQ=255)
                                                                           aTCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGAt               >  2:1702123/1‑100 (MQ=255)
                                                                            tctcACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATg              >  1:394677/1‑100 (MQ=255)
                                                                            tctcACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATg              >  1:734432/1‑100 (MQ=255)
                                                                             ctcACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCTGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGc             >  2:603572/1‑100 (MQ=255)
                                                                             ctcACAATACTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGc             >  2:640652/1‑100 (MQ=255)
                                                                              tcACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCa            >  1:1089542/1‑100 (MQ=255)
                                                                               cacaATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGTGCTGCa            <  2:888917/99‑1 (MQ=255)
                                                                                acacTGCTTTTTCTCCAGCCCACTCAGCGCGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAgg          <  2:885869/100‑1 (MQ=255)
                                                                                acaATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAgg          >  1:693817/1‑100 (MQ=255)
                                                                                acaATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAgg          <  1:1699317/100‑1 (MQ=255)
                                                                                  aaTGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGc        <  1:1681778/100‑1 (MQ=255)
                                                                                   aTGCGTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGcc       <  2:1321101/100‑1 (MQ=255)
                                                                                     gCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTc     <  1:1959996/100‑1 (MQ=255)
                                                                                      cTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCa    >  1:1624722/1‑100 (MQ=255)
                                                                                       tttttCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCaa   >  1:1959523/1‑100 (MQ=255)
                                                                                        ttttCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAAt  >  1:1391887/1‑100 (MQ=255)
                                                                                        ttttCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAAt  <  2:1786883/100‑1 (MQ=255)
                                                                                             |                                                                                              
GTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAAT  >  NC_000913/606206‑606393

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