Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,131,524 A→G 15.7% S295S (AGT→AGC yghQ ← putative inner membrane polysaccharide flippase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,131,5240AG15.7% 392.1 / 65.5 183S295S (AGT→AGCyghQputative inner membrane polysaccharide flippase
Reads supporting (aligned to +/- strand):  ref base A (64/90);  new base G (21/8);  total (85/98)
Fisher's exact test for biased strand distribution p-value = 3.90e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.24e-03

Only 100 of 198 total aligned reads displayed.

CACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACC  >  NC_000913/3131425‑3131623
                                                                                                   |                                                                                                   
 accaGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAAc                                                                                                    <  2:1474369/100‑1 (MQ=255)
    aGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACttt                                                                                                 <  2:3001021/100‑1 (MQ=255)
     gAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCCCGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACttt                                                                                                 <  1:601479/99‑1 (MQ=255)
       aTACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTAc                                                                                              <  2:2864926/100‑1 (MQ=255)
       aTACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTAc                                                                                              <  2:786622/100‑1 (MQ=255)
       aTACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTAc                                                                                              <  1:30544/100‑1 (MQ=255)
         aCCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAGTTTTTAccc                                                                                            <  2:3229340/100‑1 (MQ=255)
           cGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAGTTTTTACCCAg                                                                                          <  1:1993795/100‑1 (MQ=255)
           cGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAGTTTTTACCCAg                                                                                          <  2:2916085/100‑1 (MQ=255)
           cGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAg                                                                                          <  2:2344551/100‑1 (MQ=255)
           cGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAg                                                                                          <  1:3682696/100‑1 (MQ=255)
             aTTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCa                                                                                        <  2:3842686/100‑1 (MQ=255)
               tCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCaaa                                                                                      <  2:2444376/100‑1 (MQ=255)
               tCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGACCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCaaa                                                                                      <  1:228140/100‑1 (MQ=255)
                ccgcccgcCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAGTTTTTACCCAgaggt                                                                                      <  1:117313/99‑5 (MQ=39)
                ccgcccgcCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAAc                                                                                     >  2:2090546/1‑100 (MQ=255)
                 cgcgcgccAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAAcc                                                                                    <  1:2549216/100‑1 (MQ=255)
                 cgcccgccAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAAcc                                                                                    >  1:447209/1‑100 (MQ=255)
                 cgcccagcAGTAAACCGCACTTCACACCTAACACGCCCGGTCTGGTGGCGCGCGGATCTAAACGCATCTCCTCTAGGTAAAAGTTTTTACCCAgagg                                                                                       <  1:2908155/97‑4 (MQ=255)
                  gcccgccAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAAccc                                                                                   >  1:3282434/1‑100 (MQ=255)
                   cccgccAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCg                                                                                  >  1:1863096/1‑100 (MQ=255)
                    ccgccAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCCCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGc                                                                                 <  1:2153056/100‑1 (MQ=255)
                       ccAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCTCCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTgg                                                                              <  1:2339630/100‑1 (MQ=255)
                       ccAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTgg                                                                              <  2:527996/100‑1 (MQ=255)
                        cAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGc                                                                             >  1:3564430/1‑100 (MQ=255)
                         aGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCg                                                                            >  2:3388679/1‑100 (MQ=255)
                              aCCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTcc                                                                        >  1:2873704/1‑99 (MQ=255)
                                cGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCcggc                                                                     <  1:2480990/100‑1 (MQ=255)
                                cGGACTTCACACCTAACAACCACAATCTGGTGGTGCGCGGATCTAAACGCATCACCTCTAGGTAAAAAATTTTACCCAGCAAACCCGCTGGCg                                                                            >  1:1609280/1‑93 (MQ=255)
                                    cTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGgcggcg                                                                 <  1:3034018/100‑1 (MQ=255)
                                     ttCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGt                                                                >  2:3136389/1‑100 (MQ=255)
                                       cacaCCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCg                                                              <  2:4151644/100‑1 (MQ=255)
                                        acacCTAACAACCCCGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCga                                                             <  2:854756/100‑1 (MQ=255)
                                            cTAACAACCACGGTCTGGTGGTGCGGGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCgaagaa                                                         <  1:622728/100‑1 (MQ=255)
                                            cTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCgaagaa                                                         <  1:2688164/100‑1 (MQ=255)
                                             tAACAACCACGGTCTGGTGGTGCGCGGATATAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGCTCCGGCGGCGTCgaaga                                                          >  1:1020431/1‑98 (MQ=255)
                                                caacCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTc                                                     <  1:208864/100‑1 (MQ=255)
                                                  acCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCAt                                                   >  2:4260398/1‑100 (MQ=255)
                                                  acCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCAt                                                   >  1:2740260/1‑100 (MQ=255)
                                                   ccACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATg                                                  >  2:3544398/1‑100 (MQ=255)
                                                   ccACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATg                                                  >  1:427213/1‑100 (MQ=255)
                                                    cACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATgg                                                 >  2:4320569/1‑100 (MQ=255)
                                                     aCGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGc                                                <  2:1897629/100‑1 (MQ=255)
                                                       ggTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGa                                              >  2:56731/1‑100 (MQ=255)
                                                        gTCTGGTGGTGCGCGGATCTAAACGCGTCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGc                                                <  1:2636786/97‑1 (MQ=255)
                                                              tggtgCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTaaa                                       <  1:76564/100‑1 (MQ=255)
                                                               ggtgCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAAt                                      <  2:2195555/100‑1 (MQ=255)
                                                               ggtgCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAAt                                      >  2:72924/1‑100 (MQ=255)
                                                                gtgCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATa                                     <  2:888205/100‑1 (MQ=255)
                                                                 tgCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATaa                                    <  1:4168980/100‑1 (MQ=255)
                                                                  gcgcGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAAc                                   >  1:1331941/1‑100 (MQ=255)
                                                                  gcgaGGATCTAAACGCACCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATa                                     <  2:3987674/98‑1 (MQ=255)
                                                                     cgGATCTAAACGCATCACCTCTGGGTCAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCcg                                >  2:2676339/1‑100 (MQ=255)
                                                                     cgGATCTAAACGCATCACCTCTGGGGAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAAcc                                  >  1:3896709/1‑98 (MQ=255)
                                                                         tCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCcggcgg                            <  1:2418379/100‑1 (MQ=255)
                                                                         tCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCcggcg                             <  1:374404/99‑1 (MQ=255)
                                                                          cTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCggcggc                           >  2:3812803/1‑100 (MQ=255)
                                                                           tAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCggcggc                           >  2:1838694/1‑99 (MQ=255)
                                                                           tAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCa                          >  2:2027805/1‑100 (MQ=255)
                                                                            aaaCGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAg                         >  1:1390115/1‑100 (MQ=255)
                                                                            aaaCGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAg                         >  2:1492994/1‑100 (MQ=255)
                                                                            aaaCGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAg                         >  1:4221387/1‑100 (MQ=255)
                                                                              ncGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCg                       >  1:3591377/2‑100 (MQ=255)
                                                                               cGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGc                        <  2:1632848/98‑1 (MQ=255)
                                                                               cGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCgg                      <  1:1727789/100‑1 (MQ=255)
                                                                                gCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCggg                     <  2:4280238/100‑1 (MQ=255)
                                                                                 cATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGt                    <  2:2613201/100‑1 (MQ=255)
                                                                                 cATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGt                    >  2:3345325/1‑100 (MQ=255)
                                                                                 cATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCACCGGGt                    >  1:4065615/1‑100 (MQ=255)
                                                                                  aTCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTc                   <  2:1548949/100‑1 (MQ=255)
                                                                                   tCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTcc                  <  2:1606657/100‑1 (MQ=255)
                                                                                    cACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATACCCCGGCGGCAGCGGGTcc                  <  1:3872464/99‑1 (MQ=255)
                                                                                     aCCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTa                <  1:501098/100‑1 (MQ=255)
                                                                                      ccTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCg                       <  1:805359/92‑1 (MQ=255)
                                                                                       ctctGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAac              <  2:4024086/100‑1 (MQ=255)
                                                                                       ctctGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAac              >  1:366014/1‑100 (MQ=255)
                                                                                       ctctGGGTAAAAACATTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAac              <  1:2994839/100‑1 (MQ=255)
                                                                                        tctGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAaca             >  1:752530/1‑100 (MQ=255)
                                                                                         ctGGGTAAAAGTTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAacgc            >  2:2264840/1‑98 (MQ=255)
                                                                                         ctGGGTAAAAGTTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAacac            >  1:92574/1‑100 (MQ=255)
                                                                                         ctGGGTAAAAGTTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAacac            >  2:3869597/1‑100 (MQ=255)
                                                                                         ctGGGTAAAAGTTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAacac            >  2:4081396/1‑100 (MQ=255)
                                                                                         ctGGGTAAAAGTTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAacac            >  2:4164351/1‑100 (MQ=255)
                                                                                         ctGGGTAAAAGTTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAacac            >  2:1297294/1‑100 (MQ=255)
                                                                                         ctGGGTAAAAGTTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAacac            >  2:1255989/1‑100 (MQ=255)
                                                                                         ctGGGTAAAAGTTTTTACCCAGCAAACCCGCTAGCGTTCCGGCGGCGTCGAAGAATGGCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAaca             >  1:2066450/1‑99 (MQ=255)
                                                                                          tGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACg           <  1:4257502/100‑1 (MQ=255)
                                                                                          tGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCCTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGTGTCGTAACACg           <  1:2313017/100‑1 (MQ=255)
                                                                                           gggTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGa          >  1:3204688/1‑100 (MQ=255)
                                                                                           gggTAAAAACTTTTACCCAGCAAACACGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATATTAAATAACACGGCGGCAGCGGGTCCTAACACGa          >  1:4002480/1‑100 (MQ=255)
                                                                                           aggTAAAAACTTTTACCCAGCAAACCCGCTCGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGACCTAACACGa          >  1:2885848/2‑100 (MQ=255)
                                                                                            ggTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGAt         >  1:3348555/1‑100 (MQ=255)
                                                                                              tAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATcc       <  1:2262586/100‑1 (MQ=255)
                                                                                                naaaCTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATcccc     >  2:2767002/2‑100 (MQ=255)
                                                                                                aaaaCTTTTACCCATCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATcccc     >  1:2582494/1‑100 (MQ=255)
                                                                                                aaaaCTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATcccc     >  2:44635/1‑100 (MQ=255)
                                                                                                aaaaCTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATcccc     <  1:2172048/100‑1 (MQ=255)
                                                                                                aaaaCTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATcccc     >  2:3632043/1‑100 (MQ=255)
                                                                                                   aCTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCccacc  <  2:1952221/100‑1 (MQ=255)
                                                                                                   aCTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCccacc  <  1:3530663/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACC  >  NC_000913/3131425‑3131623

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: