Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,091,351 C→T G168D (GGC→GAC)  rhaD ← rhamnulose‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,091,3510CT99.7% 1288.4 / NA 354G168D (GGC→GAC) rhaDrhamnulose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (142/210);  total (143/211)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 367 total aligned reads displayed.

GCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGCCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATC  >  NC_000913/4091261‑4091440
                                                                                          |                                                                                         
  ggTCGCCTGGCCTATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTc                                                                                       <  1:3256849/93‑1 (MQ=255)
   gTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTcc                                                                                      >  1:5753579/1‑93 (MQ=255)
      cccTGGCCCATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCcac                                                                                   <  1:2700147/92‑1 (MQ=255)
       gcTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCcaca                                                                                  <  1:4498780/92‑1 (MQ=255)
         tGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGt                                                                                >  2:9809085/1‑93 (MQ=255)
         tGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGt                                                                                <  2:7324393/93‑1 (MQ=255)
         tGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGt                                                                                >  1:6487593/1‑93 (MQ=255)
         tGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGt                                                                                >  1:5590660/1‑93 (MQ=255)
         tGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGCCTTCCCACAGt                                                                                >  1:1854512/1‑93 (MQ=255)
          ggCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGtt                                                                               <  1:9096760/93‑1 (MQ=255)
            ccGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTgg                                                                             >  2:2604060/1‑93 (MQ=255)
            ccGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTgg                                                                             <  2:8513671/93‑1 (MQ=255)
              gATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCg                                                                           >  1:8854502/1‑93 (MQ=255)
              gATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCg                                                                           >  1:759201/1‑93 (MQ=255)
              gATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCg                                                                           <  1:5284715/93‑1 (MQ=255)
              gATTTAGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCg                                                                           >  2:5011292/1‑93 (MQ=255)
              gATGTCGTCCGTGCCGGGCTCGATCCACGGCAAATTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCg                                                                           <  1:1691593/93‑1 (MQ=255)
               aTTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGa                                                                          <  1:1060796/93‑1 (MQ=255)
               aTTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGa                                                                          >  1:4139800/1‑93 (MQ=255)
               aTTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGa                                                                          <  2:8845914/93‑1 (MQ=255)
               aTTCCGTCCGTGCAGCGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGa                                                                          <  2:9017532/93‑1 (MQ=255)
                   cGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGa                                                                      >  2:1758786/1‑93 (MQ=255)
                   cGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGa                                                                      <  2:4515306/93‑1 (MQ=255)
                     tCCGTCCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAg                                                                    <  2:3998647/93‑1 (MQ=255)
                        gTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGTAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGAcc                                                                 <  1:7145639/93‑1 (MQ=255)
                          gCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGTGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACcgc                                                               <  1:533943/93‑1 (MQ=255)
                          gCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACcgc                                                               >  2:2976141/1‑93 (MQ=255)
                           ccGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACcgcg                                                              <  1:3735257/93‑1 (MQ=255)
                            cGGGCCCCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCgg                                                             <  2:6322075/93‑1 (MQ=255)
                            cGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATCCCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCgg                                                             <  1:6208614/93‑1 (MQ=255)
                            cGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCgg                                                             <  2:4762533/93‑1 (MQ=255)
                              ggCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGgtg                                                           >  2:1451905/1‑93 (MQ=255)
                               gCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGgtgt                                                          <  2:1730615/93‑1 (MQ=255)
                               gCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGgtgt                                                          <  2:4240781/93‑1 (MQ=255)
                                cACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGCCCGCGgtgtg                                                         >  1:8340409/1‑92 (MQ=255)
                                cACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTc                                                         >  2:10031880/1‑93 (MQ=255)
                                   cATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGAc                                                                  >  2:559068/1‑81 (MQ=255)
                                    aTCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGttt                                                     <  2:4096643/93‑1 (MQ=255)
                                     tCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGtttt                                                    >  2:6951010/1‑93 (MQ=255)
                                     tCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGCTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCgagtt                                                    >  2:5332723/1‑89 (MQ=255)
                                      ccACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTc                                                   <  1:3291852/93‑1 (MQ=255)
                                      ccACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTc                                                   <  1:2352584/93‑1 (MQ=255)
                                       cACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCa                                                  >  1:7863197/1‑93 (MQ=255)
                                        ccGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCaa                                                 <  1:8471042/92‑1 (MQ=255)
                                         cGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAg                                                >  2:9380575/1‑93 (MQ=255)
                                          ggCAAAATGCCACCGCCATCCGGGAATCCCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGt                                               <  2:9193939/93‑1 (MQ=255)
                                           gCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTa                                              <  1:134996/93‑1 (MQ=255)
                                            cAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTAc                                             <  1:6660586/93‑1 (MQ=255)
                                             aaaaTGCCAACGCGATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACa                                            <  1:5098566/93‑1 (MQ=255)
                                             aaaaTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACa                                            <  1:2444727/93‑1 (MQ=255)
                                              aaaTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACAt                                           >  2:7064381/1‑93 (MQ=255)
                                               aaTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATa                                          <  1:5227994/93‑1 (MQ=255)
                                               aaTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATa                                          <  2:7399751/93‑1 (MQ=255)
                                                aTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAg                                         <  1:491120/93‑1 (MQ=255)
                                                   ccAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTg                                      <  2:7795175/93‑1 (MQ=255)
                                                    cAACGCCATCGGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGa                                     <  1:8282535/93‑1 (MQ=255)
                                                    cAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGa                                     <  2:2789468/93‑1 (MQ=255)
                                                     aaCGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAg                                    >  1:6043785/1‑93 (MQ=255)
                                                      aCGCCATCCGGGAATACCCCCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAgg                                   <  1:2973523/93‑1 (MQ=255)
                                                       cGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAggg                                  <  1:3599737/93‑1 (MQ=255)
                                                         ccATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCg                                <  2:9723262/93‑1 (MQ=255)
                                                         ccATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCg                                <  2:1145253/93‑1 (MQ=255)
                                                          cATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGa                               >  2:5452958/1‑93 (MQ=255)
                                                            tCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATc                             >  2:8330173/1‑93 (MQ=255)
                                                            tCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATc                             <  2:8008433/93‑1 (MQ=255)
                                                            tCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATc                             >  1:9746621/1‑93 (MQ=255)
                                                            tCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATc                             <  1:8901991/93‑1 (MQ=255)
                                                             ccGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCa                            >  2:2418784/1‑93 (MQ=255)
                                                              cGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAg                           <  1:9332801/93‑1 (MQ=255)
                                                                ggAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGt                         <  2:235746/93‑1 (MQ=255)
                                                                ggAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGt                         <  2:8686545/93‑1 (MQ=255)
                                                                ggAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGt                         <  1:4703907/93‑1 (MQ=255)
                                                                ggAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGt                         >  2:5883422/1‑93 (MQ=255)
                                                                ggAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGt                         <  2:5773318/93‑1 (MQ=255)
                                                                gaaaTACCCCCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGt                         <  2:6761555/91‑1 (MQ=255)
                                                                  aaTACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtg                       >  2:3215716/1‑93 (MQ=255)
                                                                   aTACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtgg                      <  1:4372082/93‑1 (MQ=255)
                                                                    tACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtggt                     >  1:5797240/1‑93 (MQ=255)
                                                                     accaccAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtggtg                    <  1:8556809/93‑1 (MQ=255)
                                                                     accaccAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtggtg                    <  1:2674132/93‑1 (MQ=255)
                                                                        accaGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGCTCAGGTTGggggcg                 >  1:1588508/1‑93 (MQ=255)
                                                                         ccaGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcgt                <  1:7274210/93‑1 (MQ=255)
                                                                         ccaGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcgt                <  1:7941730/93‑1 (MQ=255)
                                                                         ccaGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcgt                <  2:6274972/93‑1 (MQ=255)
                                                                          caGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcgtg               <  1:6922787/93‑1 (MQ=255)
                                                                          caGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcgtg               >  2:7649156/1‑93 (MQ=255)
                                                                          caGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGGTGgtggcgtg               >  1:3861421/1‑93 (MQ=255)
                                                                            gtcacTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcgtggc             <  2:3589664/91‑1 (MQ=255)
                                                                              cacTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAg           <  1:9211893/93‑1 (MQ=255)
                                                                               ccTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGt          <  2:2020776/92‑1 (MQ=255)
                                                                               acTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGt          <  2:8569269/93‑1 (MQ=255)
                                                                                cTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTg         >  1:7010690/1‑93 (MQ=255)
                                                                                 tCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGc        >  2:3751182/1‑93 (MQ=255)
                                                                                  cGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGAGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCa       >  1:206882/1‑93 (MQ=255)
                                                                                    gTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATg     >  1:1115500/1‑93 (MQ=255)
                                                                                     tgctgtCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCTGTGCATGa    <  1:8876858/93‑1 (MQ=255)
                                                                                      gctgtCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGCGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGAt   >  1:8081797/1‑93 (MQ=255)
                                                                                      gctgtCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGAt   >  1:8653685/1‑93 (MQ=255)
                                                                                      gctgtCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGAt   <  2:4713304/93‑1 (MQ=255)
                                                                                       ctgtCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATc  >  1:4948715/1‑93 (MQ=255)
                                                                                          |                                                                                         
GCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGCCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATC  >  NC_000913/4091261‑4091440

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