Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 485,069 (G)5→6 coding (85/648 nt) acrR → acrAB operon repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913485,0641.G97.1% 533.1 / 7.9 173coding (80/648 nt)acrRacrAB operon repressor
Reads supporting (aligned to +/- strand):  ref base . (3/2);  new base G (87/81);  total (90/83)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 203 total aligned reads displayed.

AATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCA‑GGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAGTCGG  >  NC_000913/484972‑485153
                                                                                             |                                                                                          
 aTTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCAcagca‑g                                                                                           >  1:2936014/1‑93 (MQ=255)
  ttCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCTATGTGGCTCTACGTCTTTTCTCACAGCA‑gg                                                                                          <  2:2388377/93‑1 (MQ=255)
  ttCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCA‑gg                                                                                          >  1:4151684/1‑93 (MQ=255)
   tCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCA‑ggg                                                                                         >  1:900640/1‑93 (MQ=255)
     aCGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCA‑ggggg                                                                                       >  2:1636071/1‑93 (MQ=255)
     aCGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCA‑ggggg                                                                                       >  2:3706774/1‑93 (MQ=255)
       gAACATATGGAACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCCCAGCA‑ggggggc                                                                                     >  2:4065792/1‑91 (MQ=255)
         aCATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCA‑ggggggtat                                                                                   >  1:3523013/1‑89 (MQ=255)
         aCATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCA‑ggggggtat                                                                                   >  1:2791039/1‑89 (MQ=255)
           atatGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑Tatca                                                                                 <  2:4162625/93‑1 (MQ=255)
            tatGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTATTTTCTCACAGCAGGGGGG‑Tatcat                                                                                >  1:4528497/1‑93 (MQ=255)
             atGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑Tatcatc                                                                               >  1:4074602/1‑93 (MQ=255)
             atGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑Tatcatc                                                                               <  1:2401534/93‑1 (MQ=255)
              tGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATcc                                                                              >  2:4421037/1‑93 (MQ=255)
               ggCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCa                                                                             <  2:3004997/93‑1 (MQ=255)
                 cACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCAcc                                                                           <  2:4098332/93‑1 (MQ=255)
                  aCGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCt                                                                          >  1:923610/1‑93 (MQ=255)
                      aaaaCCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCt                                                                      >  1:654604/1‑93 (MQ=255)
                      aaaaCCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCt                                                                      <  2:2237371/93‑1 (MQ=255)
                        aaCCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATGATCCACCTCGCTgg                                                                    >  2:1854246/1‑93 (MQ=255)
                        aaCCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTgg                                                                    <  2:1071127/93‑1 (MQ=255)
                        aaCCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTgg                                                                    >  2:3156649/1‑93 (MQ=255)
                          ccAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGc                                                                  <  2:3418496/93‑1 (MQ=255)
                          ccAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGc                                                                  >  1:4540042/1‑93 (MQ=255)
                            aaaCAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCga                                                                <  1:3283522/93‑1 (MQ=255)
                            aaaCAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCga                                                                <  1:3754349/93‑1 (MQ=255)
                             aaCAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCgag                                                               >  2:1776542/1‑93 (MQ=255)
                              aCAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCgaga                                                              <  1:1807437/93‑1 (MQ=255)
                              aCAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCgaga                                                              <  1:1469575/93‑1 (MQ=255)
                                 agaagCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTg                                                           >  1:1756524/1‑93 (MQ=255)
                                 agaagCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTg                                                           >  1:1492320/1‑93 (MQ=255)
                                  gaagCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGc                                                          <  2:619100/93‑1 (MQ=255)
                                       gcAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAag                                                     >  2:1700653/1‑93 (MQ=255)
                                        cAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAagc                                                    <  2:3278308/93‑1 (MQ=255)
                                        cAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAagc                                                    >  2:3070492/1‑93 (MQ=255)
                                        cAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAagc                                                    >  2:1343549/1‑93 (MQ=255)
                                        cAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAagc                                                    >  2:4515313/1‑93 (MQ=255)
                                        cAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAagc                                                    <  1:2120867/93‑1 (MQ=255)
                                        cAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAagc                                                    <  2:1990477/93‑1 (MQ=255)
                                         aaGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAagca                                                   <  1:30594/93‑1 (MQ=255)
                                         aaGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAagca                                                   <  2:4643837/93‑1 (MQ=255)
                                          aGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAagcag                                                  >  2:2979551/1‑93 (MQ=255)
                                          aGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAagcag                                                  <  2:4569811/93‑1 (MQ=255)
                                          aGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAagcag                                                  >  1:3345912/1‑93 (MQ=255)
                                           gAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAagcaac                                                 >  1:1713367/1‑91 (MQ=255)
                                            aaaCGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCt                                                <  2:3996352/93‑1 (MQ=255)
                                             aaCGCGCCAACACATCCTCGATGTGTCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTg                                               <  2:3982527/93‑1 (MQ=255)
                                             aaCGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTg                                               >  2:3119425/1‑93 (MQ=255)
                                             aaCGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTg                                               >  1:3885719/1‑93 (MQ=255)
                                              aCGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTgg                                              <  1:3136372/93‑1 (MQ=255)
                                               cgcgCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGc                                             <  1:3594602/93‑1 (MQ=255)
                                                 cgcCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGt                                           >  1:1070753/1‑93 (MQ=255)
                                                  gcCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGtt                                          <  1:806171/93‑1 (MQ=255)
                                                   ccAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTa                                         <  2:249714/93‑1 (MQ=255)
                                                   ccAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTa                                         >  2:4257850/1‑93 (MQ=255)
                                                   ccAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCA‑GGGGGTTATCATCCCCCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTa                                         <  2:3699863/93‑1 (MQ=255)
                                                    cAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTAc                                        <  1:1435883/93‑1 (MQ=255)
                                                    cAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTAc                                        >  2:335088/1‑93 (MQ=255)
                                                     aaCACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCATCTGGCGTTAcg                                       >  2:3200140/1‑93 (MQ=255)
                                                      acacATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTAcgc                                      <  2:2144967/93‑1 (MQ=255)
                                                         caTCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTAcgcgcg                                   >  2:1880982/1‑93 (MQ=255)
                                                          aTCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCgg                                  <  1:58874/93‑1 (MQ=255)
                                                          aTCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCgg                                  <  2:2614931/93‑1 (MQ=255)
                                                          aTCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCgg                                  >  1:3643667/1‑93 (MQ=255)
                                                           tCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCggg                                 >  2:2033484/1‑92 (MQ=255)
                                                           tCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGt                                 >  2:3593803/1‑93 (MQ=255)
                                                           tCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGt                                 >  2:3632982/1‑93 (MQ=255)
                                                           tCCTCGATGTGGCTCTACGTCTTTTCTCACAGCA‑GGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTg                                >  2:157824/1‑93 (MQ=255)
                                                              tCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCACGGTGCa                              >  2:3627666/1‑93 (MQ=255)
                                                               cGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCaa                             <  1:4291126/93‑1 (MQ=255)
                                                                gATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAAt                            <  2:4581915/93‑1 (MQ=255)
                                                                 aTGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATc                           >  2:1163993/1‑93 (MQ=255)
                                                                 aTGTGGCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATc                           <  2:454261/93‑1 (MQ=255)
                                                                     ggCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACt                       <  1:2480640/93‑1 (MQ=255)
                                                                     ggCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACt                       <  1:1672747/93‑1 (MQ=255)
                                                                     ggCTCTACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACt                       >  1:2535775/1‑93 (MQ=255)
                                                                     ggCTCTACGTCTTTTCTCACAGCA‑GGGGGTTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACt                       <  1:618458/93‑1 (MQ=255)
                                                                      gCTCTCCGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTg                      <  1:1155161/93‑1 (MQ=255)
                                                                        tctACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGCTGCAATCTACTGGc                    <  2:2863319/93‑1 (MQ=255)
                                                                         ctACGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCa                   <  1:1951795/93‑1 (MQ=255)
                                                                           aCGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCAtt                 >  1:995925/1‑93 (MQ=255)
                                                                            cGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCAttt                >  2:3497540/1‑93 (MQ=255)
                                                                            cGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCAttt                <  2:20620/93‑1 (MQ=255)
                                                                            cGTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCAttt                >  1:3650535/1‑93 (MQ=255)
                                                                             gTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCAtttt               >  2:2146429/1‑93 (MQ=255)
                                                                             gTCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCAtttt               >  1:2806333/1‑93 (MQ=255)
                                                                              tCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTa              >  1:3211731/1‑93 (MQ=255)
                                                                              tCTTTTCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTa              >  2:4203152/1‑93 (MQ=255)
                                                                                ttttCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTaaa            <  1:2736800/93‑1 (MQ=255)
                                                                                ttttCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTaaa            <  1:2253640/93‑1 (MQ=255)
                                                                                  ttCTCACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGa          <  2:1465695/93‑1 (MQ=255)
                                                                                   tctcACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGAc         >  1:1817442/1‑93 (MQ=255)
                                                                                     tcACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACaa       >  2:2349377/1‑93 (MQ=255)
                                                                                     tcACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACaa       <  1:4405332/93‑1 (MQ=255)
                                                                                     tcACAGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACaa       <  2:835384/93‑1 (MQ=255)
                                                                                      cacaGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAg      <  1:472561/93‑1 (MQ=255)
                                                                                      cacaGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAg      >  1:4375586/1‑93 (MQ=255)
                                                                                      cacaGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAg      >  1:2220767/1‑93 (MQ=255)
                                                                                       acaGCAGGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAGt     <  2:581977/93‑1 (MQ=255)
                                                                                        cagcagGGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAGTc    <  2:413159/93‑1 (MQ=255)
                                                                                             |                                                                                          
AATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCA‑GGGGG‑TATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAGTCGG  >  NC_000913/484972‑485153

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: ‑