Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,400,051 T→G E346A (GAA→GCA)  mreB ← cell wall structural complex MreBCD, actin‑like component MreB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,400,0510TG98.5% 694.8 / 0.2 197E346A (GAA→GCA) mreBcell wall structural complex MreBCD, actin‑like component MreB
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base G (107/86);  total (109/88)
Fisher's exact test for biased strand distribution p-value = 5.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

Only 100 of 199 total aligned reads displayed.

GCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACT  >  NC_000913/3399955‑3400146
                                                                                                |                                                                                               
gCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCg                                                                                               <  2:2878425/99‑1 (MQ=255)
  cACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCt                                                                                             >  1:3982780/1‑99 (MQ=255)
   aCGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTg                                                                                            >  2:2127983/1‑99 (MQ=255)
    cGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGa                                                                                           >  1:3931055/1‑99 (MQ=255)
      gCTAAGAAGTGACGTCATAAGTTATGCGTATTCGCGTATCAGACCAGGCAGGGTAAACAGGCACTGCCCCTGCCTGCATCCGATTACTCTGCGCTGAAc                                                                                         <  2:4227661/99‑1 (MQ=255)
      gCTAAAAATTGGCTTCATAAGTTATGCGTCTTCTCGTATCAGACCAGGCAGGGTAAACTGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAAc                                                                                         <  2:443963/99‑1 (MQ=255)
      gCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAAc                                                                                         <  1:1950156/99‑1 (MQ=255)
        tAAAAATTGGCTTCATAAGTTATGTGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAg                                                                                       >  1:4808506/1‑99 (MQ=255)
           aaaTTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTc                                                                                    >  1:999167/1‑99 (MQ=255)
              ttGTCTTCATAAGTTATGCGTATTCGCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTcgcc                                                                                 <  1:4543710/99‑1 (MQ=255)
               tGGCTTCATAAGTTGTGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTCCCCCTGCCTGCATCCGATTACCCTGCGCTGCACAGGTcgccg                                                                                <  2:1343703/99‑1 (MQ=255)
               tGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTcgccg                                                                                <  1:2203917/99‑1 (MQ=255)
                 gCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCgccgcc                                                                              <  2:907495/99‑1 (MQ=255)
                   ttCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGt                                                                            <  2:5206708/99‑1 (MQ=255)
                   ttCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGt                                                                            <  1:3780338/99‑1 (MQ=255)
                      aTAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCa                                                                         <  2:2453527/99‑1 (MQ=255)
                      aTAAGTTATGCGGATTCTCGTATCAGACCAGGCAGCGTAAACAGACACATCCCCTGCCTGCATCCGAATACTATGCGCTGGACAGGTCGCCGCCGTCCa                                                                         >  2:4490530/1‑99 (MQ=255)
                        aaGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATg                                                                       >  1:3518965/1‑99 (MQ=255)
                         aGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGt                                                                      <  2:2380136/99‑1 (MQ=255)
                            tATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGa                                                                   >  1:3212518/1‑99 (MQ=255)
                             aTGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACGTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGAt                                                                  <  2:3249195/99‑1 (MQ=255)
                               gCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCa                                                                <  2:5390217/99‑1 (MQ=255)
                                cGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCAt                                                               >  1:2803718/1‑99 (MQ=255)
                                cGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCAt                                                               >  1:3380964/1‑99 (MQ=255)
                                 gTATTCTCGTATCAGACGAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCTGATTACTCTGCGATGATCAGGTCGACGACGTGCATGTCGATCAtt                                                              >  2:1023684/1‑99 (MQ=255)
                                 gTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCAtt                                                              <  1:7474/99‑1 (MQ=255)
                                 gTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCAtt                                                              >  1:3672749/1‑99 (MQ=255)
                                 gTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCAtt                                                              >  2:3218752/1‑99 (MQ=255)
                                   aTTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTc                                                            <  1:4666190/99‑1 (MQ=255)
                                    ttCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTcc                                                           <  2:3794606/99‑1 (MQ=255)
                                     tctcGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGAGTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCa                                                          <  1:1611810/99‑1 (MQ=255)
                                      ctcGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTAGCCGCCGTGCATGTCGATCATTTCCAg                                                         >  1:1336901/1‑99 (MQ=255)
                                       tcGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAgc                                                        <  1:836182/99‑1 (MQ=255)
                                       tcGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAgc                                                        >  1:940474/1‑99 (MQ=255)
                                       tcGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAgc                                                        >  1:4729637/1‑99 (MQ=255)
                                          tATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCt                                                     >  1:5556229/1‑99 (MQ=255)
                                          tATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCt                                                     <  1:2866727/99‑1 (MQ=255)
                                          tATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCt                                                     >  1:3868437/1‑99 (MQ=255)
                                           aTCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCtt                                                    >  2:1935657/1‑99 (MQ=255)
                                            tCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCttt                                                   <  2:631612/99‑1 (MQ=255)
                                            tCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCttt                                                   <  2:1260196/99‑1 (MQ=255)
                                               gACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTgcc                                                >  2:1642484/1‑99 (MQ=255)
                                               gACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTgcc                                                <  1:308366/99‑1 (MQ=255)
                                               gACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTgcc                                                >  1:5250514/1‑99 (MQ=255)
                                               gACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCAATCATTTCCAGCGCTTTgcc                                                <  2:5471870/99‑1 (MQ=255)
                                               gACCAGGCAGGGTAAACAGACACGTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTgcc                                                <  1:1113148/99‑1 (MQ=255)
                                                aCCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTgccg                                               >  2:4825848/1‑99 (MQ=255)
                                                aCCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTgccg                                               >  1:4163059/1‑99 (MQ=255)
                                                 tctggcaggGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTgccgc                                              <  1:2389625/96‑1 (MQ=255)
                                                 ccAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTgccgc                                              <  2:1496510/99‑1 (MQ=255)
                                                   aggcaggGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCa                                            >  2:799093/1‑99 (MQ=255)
                                                   aggcacgGTAAACAGACACTCCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCCTGCATGTCGATCATTTCCAGCGCTTTGCCGCCa                                            <  1:5366778/99‑1 (MQ=255)
                                                    ggcaggGTAAACAGACACTTCCCCTGCCTGCATCCGCTTACTCTGCGCTGACCAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCAc                                           <  2:464486/99‑1 (MQ=255)
                                                     gcatgGTAAACAGGCACTTCCCCTGCCTTCTTCCGATTTCTCTGCGTTGAAGAGGTCGCCGCCGTGCATGTCGTTTTTTTCCAGCGCTCGGCGGCCAcc                                          <  2:2199484/99‑1 (MQ=255)
                                                       aggGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACcgc                                        <  2:4712219/99‑1 (MQ=255)
                                                         ggTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACcgcgc                                      <  2:5130131/99‑1 (MQ=255)
                                                         ggTAAACAGACACTTACCCTGCCAGCATCCGATTCATCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACcgcgc                                      <  2:3456458/99‑1 (MQ=255)
                                                          gTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACcgcgcg                                     >  1:315804/1‑99 (MQ=255)
                                                           tAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCgcgcgc                                    >  2:2865827/1‑99 (MQ=255)
                                                            aaaCAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGcc                                   <  1:484590/99‑1 (MQ=255)
                                                            aaaCAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGcc                                   >  2:3213979/1‑99 (MQ=255)
                                                             aaCAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCca                                  >  2:1711511/1‑99 (MQ=255)
                                                             aaCAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCca                                  <  2:1556343/99‑1 (MQ=255)
                                                               cAGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCcaca                                >  2:760878/1‑99 (MQ=255)
                                                                aGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCcacac                               >  1:619775/1‑99 (MQ=255)
                                                                aGACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCcacac                               >  1:1020374/1‑99 (MQ=255)
                                                                 gACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCcacaca                              <  1:5267310/99‑1 (MQ=255)
                                                                 gACACTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCCCCCCCGCGCGCcacaca                              >  2:3159410/1‑99 (MQ=255)
                                                                  acacTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAg                             <  1:1047241/99‑1 (MQ=255)
                                                                  acacTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAg                             <  1:2057887/99‑1 (MQ=255)
                                                                  acacTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAg                             <  1:2519155/99‑1 (MQ=255)
                                                                    acctCCCCTGCCTGCATCCGATTCCTCTGCGCTGCACAGGTCGCCCCCGCGCATGCCGCTCATTTCCCGCGCTTTCCCCCCACCGCCCGCCACACAGGt                           >  2:5237811/4‑99 (MQ=255)
                                                                    acTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGt                           <  1:2909577/99‑1 (MQ=255)
                                                                     cTTCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTc                          >  1:5031876/1‑99 (MQ=255)
                                                                      ttCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCa                         >  2:181730/1‑99 (MQ=255)
                                                                      ttCCCCTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCa                         >  2:3557149/1‑99 (MQ=255)
                                                                        ccccTGCCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATGTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGc                       <  2:438178/99‑1 (MQ=255)
                                                                         cccGGCCTGCGTCCGGTTGCTGTGCGCTGAACAGGTCGCCGCCGGGCGTGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCg                      <  1:4690122/99‑1 (MQ=255)
                                                                          cctgcctgCATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCgg                     <  1:2622874/99‑1 (MQ=255)
                                                                           ctgcctgcATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCggg                    >  1:5402584/1‑99 (MQ=255)
                                                                            tgcctgcATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGt                   >  1:905830/1‑99 (MQ=255)
                                                                            tgcctgcATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGAGCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGt                   <  2:1168245/99‑1 (MQ=255)
                                                                             gcctgcATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGTCGCCACCGCGCGCCACACAGGTCAGCGGGTc                  <  2:4968956/99‑1 (MQ=255)
                                                                             gcctgcATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTc                  <  1:1601699/99‑1 (MQ=255)
                                                                              cctgcATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACCCAGGTCAGCGGGTCt                 >  1:5293102/1‑99 (MQ=255)
                                                                               ctgcATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGGGGGTGtt                >  1:4224768/1‑99 (MQ=255)
                                                                                tgcATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTc               <  1:412088/99‑1 (MQ=255)
                                                                                  cATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAg             >  1:3403889/1‑99 (MQ=255)
                                                                                  cATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAg             >  1:3005399/1‑99 (MQ=255)
                                                                                  cATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAg             >  1:440742/1‑99 (MQ=255)
                                                                                  cATCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAg             >  1:4905481/1‑99 (MQ=255)
                                                                                   aTCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGc            <  1:1756630/99‑1 (MQ=255)
                                                                                   aTCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGc            <  1:4703219/99‑1 (MQ=255)
                                                                                    tCCGATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGca           >  1:4788269/1‑99 (MQ=255)
                                                                                       gATTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGcaaca        >  1:430781/1‑99 (MQ=255)
                                                                                        aTTACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGcaacaa       <  1:3385709/99‑1 (MQ=255)
                                                                                         ttACTCTGCGCTGAACAGGTCGCCGCCCTGCAGGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCaacaac      <  2:3327794/99‑1 (MQ=255)
                                                                                          tACTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACg     <  2:2680886/99‑1 (MQ=255)
                                                                                           aCTCTGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGa    >  2:1596941/1‑99 (MQ=255)
                                                                                             tctGCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGCTCTTCAGCAACAACGACt  >  2:2969236/1‑99 (MQ=255)
                                                                                                |                                                                                               
GCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACT  >  NC_000913/3399955‑3400146

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