Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,533,519 T→C T122A (ACG→GCG)  yddE ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,533,5190TC99.6% 839.4 / NA 259T122A (ACG→GCG) yddEhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (111/147);  minor base A (1/0);  total (112/147)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 272 total aligned reads displayed.

AGCCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGTACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGTGC  >  NC_000913/1533429‑1533608
                                                                                          |                                                                                         
aGCCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCAc                                                                                         >  2:2469887/1‑93 (MQ=255)
aGCCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCAc                                                                                         >  2:2778308/1‑93 (MQ=255)
 gCCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCAcc                                                                                        <  1:5332251/93‑1 (MQ=255)
 gCCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCAcc                                                                                        <  2:283518/93‑1 (MQ=255)
    tGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGGCCGGCGCACCTTg                                                                                     >  2:3431264/1‑93 (MQ=255)
    tGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTg                                                                                     <  1:4542628/93‑1 (MQ=255)
    tGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTg                                                                                     >  1:2924100/1‑93 (MQ=255)
     ggCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGt                                                                                    <  1:661873/93‑1 (MQ=255)
     ggCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGt                                                                                    >  2:1700466/1‑93 (MQ=255)
      gCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGtt                                                                                   <  1:1328037/93‑1 (MQ=255)
        aGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTcc                                                                                 <  1:5209761/93‑1 (MQ=255)
        aGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTcc                                                                                 <  2:4432803/93‑1 (MQ=255)
        aGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTcc                                                                                 >  2:5002426/1‑93 (MQ=255)
         gAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCa                                                                                >  2:1277167/1‑93 (MQ=255)
            tatCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCg                                                                             >  1:1469066/1‑93 (MQ=255)
            tatCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCg                                                                             >  1:162826/1‑93 (MQ=255)
              tCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGaa                                                                           >  1:3822143/1‑93 (MQ=255)
               cGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGaaa                                                                          <  1:2597357/93‑1 (MQ=255)
               cGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGaaa                                                                          <  2:4796348/93‑1 (MQ=255)
                gTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAAt                                                                         <  2:5109098/93‑1 (MQ=255)
                 tCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATa                                                                        >  1:6166772/1‑93 (MQ=255)
                  ccTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGACCGGCGCACCTTGTTCCCGCGAAATAc                                                                       >  2:5403035/1‑93 (MQ=255)
                  ccTCAGTGAGAAGAAGCGCGTTTATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATAc                                                                       <  2:5068947/93‑1 (MQ=255)
                   cTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATaag                                                                      >  1:6051641/1‑91 (MQ=255)
                   cTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACg                                                                      >  2:2257687/1‑93 (MQ=255)
                    tCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGa                                                                     >  2:6247218/1‑93 (MQ=255)
                    tCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGa                                                                     >  1:3596112/1‑93 (MQ=255)
                    tCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGa                                                                     >  1:4091327/1‑93 (MQ=255)
                    tCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGa                                                                     >  2:3197446/1‑93 (MQ=255)
                     cAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGAt                                                                    <  1:130799/93‑1 (MQ=255)
                     cAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGAt                                                                    <  1:1090437/93‑1 (MQ=255)
                      aGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATa                                                                   >  1:3042708/1‑93 (MQ=255)
                      aGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATa                                                                   >  1:5284222/1‑93 (MQ=255)
                       gttagaTGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATaa                                                                  <  1:1218200/90‑1 (MQ=255)
                       gTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATaa                                                                  >  2:3554147/1‑93 (MQ=255)
                       gTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATaa                                                                  <  1:2239751/93‑1 (MQ=255)
                        tGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAat                                                                 <  2:594697/93‑1 (MQ=255)
                        tGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAat                                                                 <  2:663372/93‑1 (MQ=255)
                          agaTGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAatca                                                               <  1:1143691/93‑1 (MQ=255)
                          agaTGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAatca                                                               <  1:3780059/93‑1 (MQ=255)
                           gaTGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAatcat                                                              <  1:565704/93‑1 (MQ=255)
                           gaTGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAatcat                                                              >  1:318024/1‑93 (MQ=255)
                            aTGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCTGCGAAATACGATAatcatc                                                             <  1:6533303/93‑1 (MQ=255)
                            aTGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAatcatc                                                             <  1:2071313/93‑1 (MQ=255)
                             tGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCg                                                            >  1:56920/1‑93 (MQ=255)
                             tGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCg                                                            >  1:5818434/1‑93 (MQ=255)
                             tGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCg                                                            >  1:3978399/1‑93 (MQ=255)
                             tGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCg                                                            >  2:6490556/1‑93 (MQ=255)
                              gAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGt                                                           >  2:2023951/1‑93 (MQ=255)
                               aaGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGtt                                                          <  2:867118/93‑1 (MQ=255)
                                aGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTtg                                                         <  1:4802327/93‑1 (MQ=255)
                                aGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTtg                                                         >  1:5618027/1‑93 (MQ=255)
                                aGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTtg                                                         >  1:6343547/1‑93 (MQ=255)
                                 gcgcGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTtgt                                                        <  2:2732395/93‑1 (MQ=255)
                                 gcgcGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTtgt                                                        <  2:1320175/93‑1 (MQ=255)
                                  cgcgTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTtgtg                                                       <  2:1543503/93‑1 (MQ=255)
                                   gcgTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCTGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGc                                                      <  2:4697806/93‑1 (MQ=255)
                                   gcgTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTtgtgt                                                      >  2:5573188/1‑92 (MQ=255)
                                    cgTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCt                                                     <  2:1053800/93‑1 (MQ=255)
                                     gTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCtt                                                    <  1:3542461/93‑1 (MQ=255)
                                      ttGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCttt                                                   >  1:6447979/1‑93 (MQ=255)
                                       tGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCtttt                                                  <  1:1890252/93‑1 (MQ=255)
                                        gATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTc                                                 >  1:2813706/1‑93 (MQ=255)
                                           aaTCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGAt                                              <  1:165226/93‑1 (MQ=255)
                                            aTCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATa                                             <  1:5951936/93‑1 (MQ=255)
                                            aTCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATa                                             >  2:4723720/1‑93 (MQ=255)
                                            aTCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATa                                             >  2:2847162/1‑93 (MQ=255)
                                            aTCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATa                                             <  2:5237283/93‑1 (MQ=255)
                                            aTCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATa                                             <  1:805377/93‑1 (MQ=255)
                                            aTCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATa                                             <  2:5340370/93‑1 (MQ=255)
                                               gccgcACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTc                                          >  2:3879478/1‑93 (MQ=255)
                                                     cGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGa                                    <  2:2854218/93‑1 (MQ=255)
                                                     cGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGa                                    <  2:1854389/93‑1 (MQ=255)
                                                          ttCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGttt                               <  1:5108068/93‑1 (MQ=255)
                                                           tCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGtttt                              >  1:5385438/1‑93 (MQ=255)
                                                             aCCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTcc                            <  2:5690407/93‑1 (MQ=255)
                                                              ccTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCt                           >  1:2473643/1‑93 (MQ=255)
                                                                ttCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGc                         <  2:4038312/93‑1 (MQ=255)
                                                                 tCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGcc                        <  1:90712/93‑1 (MQ=255)
                                                                  ccAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCa                       <  2:2586689/93‑1 (MQ=255)
                                                                     gCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCg                    >  2:4269879/1‑93 (MQ=255)
                                                                     gCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCg                    >  1:3388/1‑93 (MQ=255)
                                                                     gCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCg                    >  1:3135150/1‑93 (MQ=255)
                                                                         tGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGt                >  1:4005811/1‑93 (MQ=255)
                                                                          ggCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGtt               >  1:4461761/1‑93 (MQ=255)
                                                                           gCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGttt              <  2:6136525/93‑1 (MQ=255)
                                                                              cAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGcc           >  2:2903714/1‑93 (MQ=255)
                                                                               aaaGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCa          >  2:2131488/1‑93 (MQ=255)
                                                                               aaaGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCa          >  2:5890793/1‑93 (MQ=255)
                                                                                aaGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAg         >  2:6555630/1‑93 (MQ=255)
                                                                                 aGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGTCAGa        >  1:4480486/1‑93 (MQ=255)
                                                                                 aGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGa        <  1:1432751/93‑1 (MQ=255)
                                                                                 aGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGa        <  2:2260931/93‑1 (MQ=255)
                                                                                  gCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGAt       <  1:4951270/93‑1 (MQ=255)
                                                                                     cGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGt    <  2:5841613/93‑1 (MQ=255)
                                                                                     cGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGt    <  2:3146055/93‑1 (MQ=255)
                                                                                     cGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGt    >  2:3746183/1‑93 (MQ=255)
                                                                                      ggCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGTg   >  1:6727830/1‑93 (MQ=255)
                                                                                       gcgcACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGTGc  >  1:3128953/1‑93 (MQ=255)
                                                                                       gcgcACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGTGc  <  2:3915194/93‑1 (MQ=255)
                                                                                          |                                                                                         
AGCCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGTACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGTGC  >  NC_000913/1533429‑1533608

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