Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,164,998 G→A *241* (TAG→TAA)  baeR → response regulator in two‑component regulatory system with BaeS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,164,9980GA97.6% 265.3 / ‑1.4 82coding (723/723 nt)baeRresponse regulator in two‑component regulatory system with BaeS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (40/40);  minor base . (1/1);  total (41/41)
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.55e-01

TAACCGACCGCACCATCGACAGCCACATTAAAAACCTGCGCCGCAAGCTGGAATCTCTCGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAGTTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAG‑CAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGG  >  NC_000913/2164862‑2165133
                                                                                                                                         |                                                                                                                                        
taACCGACCGCACCATCGACAGCCACATTAAAAACCTGCGCCGCAAGCTGGAATCTCTCGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAAtt                                                                                                                                        >  2:1649372/1‑139 (MQ=255)
        cgCACCATCGACAGCCACATTAAAAACCTGCGCCGCAAGCTGGAATCTCTCGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGAc                                                                                                                                >  1:1748142/1‑139 (MQ=255)
        cgCACCATCGACAGCCACATTAAAAACCTGCGCCGCAAGCTGGAATCTCTCGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGAc                                                                                                                                >  1:1726712/1‑139 (MQ=255)
         gCACCATCGACAGCCACATTAAAAACCTGCGCCGCAAGCTGGAATCTCTCGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACa                                                                                                                               <  2:1726712/139‑1 (MQ=255)
          cACCATCGACAGCCACATTAAAAACCTGCGCCGCAAGCTGGAGTCTCTCGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAAGTTTAGCGACat                                                                                                                              >  2:1228049/1‑139 (MQ=255)
               tCGACAGCCACATTAAAAACCTGCGCCGCAAGCTGGAATCTCTCGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTAttt                                                                                                                         >  1:1291822/1‑139 (MQ=255)
                        acaTTAAAAACCTGCGCCGCAAGCTGGAATCTCTCGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCg                                                                                                                >  2:735918/1‑139 (MQ=255)
                         caTTAAAAACCTGCGCCGCAAGCTGGAATCTCTCGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAg                                                                                                                                    >  1:252651/1‑118 (MQ=255)
                         caTTAAAAACCTGCGCCGCAAGCTGGAATCTCTCGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAg                                                                                                                                    <  2:252651/118‑1 (MQ=255)
                         caTTAAAAACCTGCGCCGCAAGCTGGAATCTCTCGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCgg                                                                                                               <  2:1291822/139‑1 (MQ=255)
                             aaaaaCCTGCGCCGCAAGCTGGAATCTCTCGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGAt                                                                                                           <  2:1748142/139‑1 (MQ=255)
                              aaaaCCTGCGCCGCAAGCTGGAATCTCTCGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATg                                                                                                          <  1:735918/139‑1 (MQ=255)
                                   cTGCGCCGCAAGCTGGAATCTCTCGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGtt                                                                                                     <  2:1372470/139‑1 (MQ=255)
                                      cgcggcAAGCTGGAATCTCTTGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGGCGGAGATGGCGTTCAg                                                                                                  >  2:1277548/1‑139 (MQ=255)
                                      cgccgcAAGCTGGAATCTCTCGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTTGGTTACCGCTGGGAAGCTGACGCCTGCCGCATCGTTTAATTTTAGCGa                                                                                                                                 >  2:409670/1‑108 (MQ=255)
                                      cgccgcAAGCTGGAATCTCTCGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAg                                                                                                                                    >  1:231961/1‑105 (MQ=255)
                                      cgccgcAAGCTGGAATCTCTCGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAg                                                                                                                                    <  2:231961/105‑1 (MQ=255)
                                      cgccgcAAGCTGGAATCTCTCGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATg                                                                                                          <  1:1529208/131‑1 (MQ=255)
                                      cgccgcAAGCTGGAATCTCTCGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATg                                                                                                          <  2:409260/131‑1 (MQ=255)
                                      cgccgcAAGCTGGAATCTCTCGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATg                                                                                                          <  2:1062150/131‑1 (MQ=255)
                                      cgccgcAAGCTGGAATCTCTCGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATg                                                                                                          >  2:1529208/1‑131 (MQ=255)
                                      cgccgcAAGCTGGAATCTCTCGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATg                                                                                                          >  1:1062150/1‑131 (MQ=255)
                                      cgccgcAAGCTGGAATCTCTCGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATg                                                                                                          <  1:641782/131‑1 (MQ=255)
                                      cgccgcAAGCTGGAATCTCTCGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATg                                                                                                          >  2:641782/1‑131 (MQ=255)
                                      cgccgcAAGCTGGAATCTCTCGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATg                                                                                                          >  1:409260/1‑131 (MQ=255)
                                      cgccgcAAGCTGGAATCTCTCGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAg                                                                                                  >  2:163706/1‑139 (MQ=255)
                                      cgccgcAAGCTGGAATCTCTCGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAg                                                                                                  >  2:1093991/1‑139 (MQ=255)
                                      cgccgcAAGCTGGAATCTCTCGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAg                                                                                                  >  2:216180/1‑139 (MQ=255)
                                      cgccgcAAGCTGGAATCTCTCGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCACATCGTTTAATTTTAGCGa                                                                                                                                 <  1:409670/108‑1 (MQ=255)
                                            aaGCTGGAATCTCTCGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATg                                                                                                          <  1:1239066/125‑1 (MQ=255)
                                            aaGCTGGAATCTCTCGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATg                                                                                                          >  2:1239066/1‑125 (MQ=255)
                                                    aTCTCTCGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCTTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCaa                                                                                    <  1:1719351/139‑1 (MQ=255)
                                                        ctcGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCgg                                                                                       >  2:1752885/1‑132 (MQ=255)
                                                        ctcGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCgg                                                                                       <  1:1752885/132‑1 (MQ=255)
                                                            aCGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAg                                                                                                                                    >  1:1284819/1‑83 (MQ=255)
                                                            aCGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAg                                                                                                                                    <  2:1284819/83‑1 (MQ=255)
                                                                   ccAGTCCTTTTTCCCGCGCCGTTTAT‑GCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCAt                                                                     <  2:380563/138‑1 (MQ=255)
                                                                          tttAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATg                                                                                                          <  2:712399/95‑1 (MQ=255)
                                                                          tttAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATg                                                                                                          >  1:712399/1‑95 (MQ=255)
                                                                            tAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAG‑Caa                                                           >  1:1138685/1‑139 (MQ=255)
                                                                                  gcgcCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAAtttt                                                                                                                                      >  1:590516/1‑60 (MQ=255)
                                                                                  gcgcCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAAtttt                                                                                                                                      <  2:590516/60‑1 (MQ=255)
                                                                                  gcgcCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAG‑CAAAc                                                         >  1:152894/1‑136 (MQ=255)
                                                                                  gcgcCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAG‑CAAAc                                                         <  2:152894/136‑1 (MQ=255)
                                                                                  gcgcCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAG‑CAAACTGa                                                      >  2:448962/1‑139 (MQ=255)
                                                                                  gcgcCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAG‑CAAACTGa                                                      >  2:128347/1‑139 (MQ=255)
                                                                                  gcgcCGTTTATGCCG‑CCGTTACCGCCGGGAAGCCCGCCCCTGGCGGCTCGGTTTATTTTAGCCGCATTATTTTGTTTGGCGGGGGTTGCGGTCAGACCAGATCCGGCAAAATTTATCCACGCAAGGTCAGCCAAACTgt                                                      >  1:380563/1‑138 (MQ=255)
                                                                                    gcCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTA‑ATTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGa                                                                                                 <  1:83446/94‑1 (MQ=255)
                                                                                    gcCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTA‑ATTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGa                                                                                                 >  2:83446/1‑94 (MQ=255)
                                                                                      cGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAG‑CAAACTGACAtt                                                  <  2:322564/139‑1 (MQ=255)
                                                                                           aTGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAG‑CAAACTGACATTAtctc                                             >  2:1079949/1‑139 (MQ=255)
                                                                                              gCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACa                                                                                  >  1:913383/1‑100 (MQ=255)
                                                                                              gCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACa                                                                                  <  2:913383/100‑1 (MQ=255)
                                                                                                gTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCgg                                                                                       <  1:1710230/93‑1 (MQ=255)
                                                                                                gTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCgg                                                                                       >  2:1710230/1‑93 (MQ=255)
                                                                                                gTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAG‑CAAACTGACATTATCTCTCCGc                                        >  2:1263394/1‑139 (MQ=255)
                                                                                                gTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAG‑CAAACTGACATTATCTCTCCGc                                        >  1:431532/1‑139 (MQ=255)
                                                                                                  cGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAG‑CAAACTGACATTATCTCTCCGccc                                      >  1:387073/1‑139 (MQ=255)
                                                                                                      tACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAG‑CAAACTGACATTATCTCTCCGCCCTggc                                  <  1:1649372/139‑1 (MQ=255)
                                                                                                             gggAAGCCGACGCCTGCCGCATCGTTTAATTTTAg                                                                                                                                    >  1:530000/1‑35 (MQ=38)
                                                                                                             gggAAGCCGACGCCTGCCGCATCGTTTAATTTTAg                                                                                                                                    <  2:530000/35‑1 (MQ=255)
                                                                                                              ggAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATcc                                                                            <  1:792247/90‑1 (MQ=255)
                                                                                                              ggAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATcc                                                                            >  2:792247/1‑90 (MQ=255)
                                                                                                              ggAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATc                                                                             >  1:719348/1‑89 (MQ=255)
                                                                                                              ggAAGCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATc                                                                             <  2:719348/89‑1 (MQ=255)
                                                                                                                  gCCGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAG‑CAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGaaa                      <  1:448962/139‑1 (MQ=255)
                                                                                                                   ccGACGCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAG‑CAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGaaaa                     <  1:1317856/139‑1 (MQ=255)
                                                                                                                    cGACGCCTGCCGCATCGTTCAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAATATTATCCCACGCATGGTCAG‑CAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAAt                    <  1:1228049/139‑1 (MQ=255)
                                                                                                                        gCCTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAG‑CAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCt                <  1:163706/139‑1 (MQ=255)
                                                                                                                         ccTGCCGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAG‑CAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTc               <  1:1297424/139‑1 (MQ=255)
                                                                                                                             ccGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGtt                                                                                                     >  1:1191036/1‑50 (MQ=255)
                                                                                                                             ccGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGtt                                                                                                     <  2:1191036/50‑1 (MQ=255)
                                                                                                                             ccGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGtt                                                                                                     >  1:1629185/1‑50 (MQ=255)
                                                                                                                             ccGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGtt                                                                                                     <  2:1629185/50‑1 (MQ=255)
                                                                                                                             ccGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACa                                                                                  <  2:418578/69‑1 (MQ=255)
                                                                                                                             ccGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACa                                                                                  >  1:418578/1‑69 (MQ=255)
                                                                                                                             ccGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACAt                                                                                 <  2:104033/70‑1 (MQ=255)
                                                                                                                             ccGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACAt                                                                                 >  1:104033/1‑70 (MQ=255)
                                                                                                                              cGCATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAG‑CAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACt          <  1:216180/139‑1 (MQ=255)
                                                                                                                                cATCGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAG‑CAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACttt        >  1:198679/1‑139 (MQ=255)
                                                                                                                                   cGTTTAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAG‑CAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAAt     <  2:387073/139‑1 (MQ=255)
                                                                                                                                      ttAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAG‑CAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCgg  <  1:1277548/139‑1 (MQ=255)
                                                                                                                                      ttAATTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAG‑CAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCgg  <  1:1019034/139‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                        
TAACCGACCGCACCATCGACAGCCACATTAAAAACCTGCGCCGCAAGCTGGAATCTCTCGACGCCGAACAGTCATTTAT‑CCGCGCCGTTTATGGCGTCGGTTACCGCTGGGAAGCCGACGCCTGCCGCATCGTTTAGTTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAG‑CAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGG  >  NC_000913/2164862‑2165133

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