Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,447,869 Δ1 bp coding (3758/3780 nt) tamB → translocation and assembly module for autotransporter export, inner membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,447,8670T.99.1% 433.5 / ‑1.5 108coding (3756/3780 nt)tamBtranslocation and assembly module for autotransporter export, inner membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (56/50);  minor base A (1/0);  total (58/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 109 total aligned reads displayed.

AGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTT  >  NC_000913/4447730‑4447997
                                                                                                                                         |                                                                                                                                  
aGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGAtt                                                                                                                                   >  2:285239/1‑139 (MQ=255)
            ttGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGtt                                                                                                                      <  1:1277525/139‑1 (MQ=255)
              gCCAGGTCTGCAAGTGAAATACGGAGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCg                                                                                                                    >  1:1395559/1‑139 (MQ=255)
                  ggTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATGGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGtt                                                                                                                >  1:392312/1‑139 (MQ=255)
                  ggTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTGTCAGTTCGAGtt                                                                                                                >  1:544459/1‑139 (MQ=255)
                  ggTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGtt                                                                                                                >  2:1060481/1‑139 (MQ=255)
                    tCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCCCTTGTTCCTTAGGCTTTTTCTGAAGCCCTTTTTTGGTGTAACCAAGGCCTGGGATGGTTTTATAAGTCCGAGtttt                                                                                                              >  2:134521/1‑139 (MQ=255)
                        cAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACa                                                                                                          <  1:426924/139‑1 (MQ=255)
                          aGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGtt                                                                                                                      >  1:1636675/1‑125 (MQ=255)
                          aGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGtt                                                                                                                      <  2:1636675/125‑1 (MQ=255)
                                  aCGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAatat                                                                                                <  1:285239/139‑1 (MQ=255)
                                     gCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGtt                                                                                                                      >  1:657779/1‑114 (MQ=255)
                                     gCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGtt                                                                                                                      <  2:657779/114‑1 (MQ=255)
                                          ggTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTAc                                                                                        >  1:459704/1‑139 (MQ=255)
                                            tataTTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTGTCAGTTCGAGTTTTAACAATGCGAATGTTTGTCTACgg                                                                                      <  2:544459/139‑1 (MQ=255)
                                                     ctctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGt                                                                                  >  2:1129169/1‑134 (MQ=255)
                                                     ctctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGt                                                                                  <  1:1129169/134‑1 (MQ=255)
                                                     ctctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGCTTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACg                                                                             >  1:1014372/1‑139 (MQ=255)
                                                      tctATAGCAACACTCGCGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCGCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGc                                                                            >  1:1374586/1‑139 (MQ=255)
                                                        tataGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCca                                                                          >  1:163869/1‑139 (MQ=255)
                                                        tataGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCca                                                                          >  1:250692/1‑139 (MQ=255)
                                                               acacTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaacaa                                                                   >  1:500278/1‑139 (MQ=255)
                                                                  cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGtt                                                                                                                      >  1:690937/1‑85 (MQ=255)
                                                                  cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGtt                                                                                                                      <  2:690937/85‑1 (MQ=255)
                                                                  cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGc                                                                >  1:199644/1‑139 (MQ=255)
                                                                  cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGc                                                                >  1:1441723/1‑139 (MQ=255)
                                                                  cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGc                                                                >  2:763674/1‑139 (MQ=255)
                                                                    cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGtt                                                                                                                      >  1:730617/1‑83 (MQ=255)
                                                                    cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGtt                                                                                                                      <  2:730617/83‑1 (MQ=255)
                                                                    cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                              <  2:1014372/139‑1 (MQ=255)
                                                                    cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                              <  2:392312/139‑1 (MQ=255)
                                                                    cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                              <  1:763674/139‑1 (MQ=255)
                                                                    cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                              <  1:1692189/139‑1 (MQ=255)
                                                                    cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                              <  1:1060481/139‑1 (MQ=255)
                                                                       gttacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAg                                                           >  1:1042857/1‑139 (MQ=255)
                                                                         tacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTc                                                         <  2:500278/139‑1 (MQ=255)
                                                                               aTCGCCTGATGCCTAAGCTATATCTGGTAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGtt                                                                                                                      <  1:1452909/72‑1 (MQ=255)
                                                                               aTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGtt                                                                                                                      >  2:1452909/1‑72 (MQ=255)
                                                                                tCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTAt                                                                                                                           >  2:1512820/1‑66 (MQ=255)
                                                                                tCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                              >  2:1559761/1‑127 (MQ=255)
                                                                                tCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                              <  1:1559761/127‑1 (MQ=255)
                                                                                 cGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGtt                                                                                                                      <  1:1076329/70‑1 (MQ=255)
                                                                                   ccTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGAc                                               <  2:1441723/139‑1 (MQ=255)
                                                                                   ccTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGAc                                               >  2:1002318/1‑139 (MQ=255)
                                                                                     tGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCa                                             >  2:586901/1‑139 (MQ=255)
                                                                                     tGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGCCACTGGATGACCa                                             >  1:127488/1‑139 (MQ=255)
                                                                                      gATGCCTAAAGTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCaa                                            <  2:127488/139‑1 (MQ=255)
                                                                                          ccTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGcc                                        <  2:1222439/139‑1 (MQ=255)
                                                                                                   tatCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaa                                                                      >  2:1208537/1‑100 (MQ=255)
                                                                                                   tatCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaa                                                                      <  1:1208537/100‑1 (MQ=255)
                                                                                                   tatCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTg                               >  2:1198509/1‑139 (MQ=255)
                                                                                                    atCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCaa                                            <  1:1175035/125‑1 (MQ=255)
                                                                                                    atCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCaa                                            >  2:1175035/1‑125 (MQ=255)
                                                                                                         gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGAc                                               >  1:1668009/1‑117 (MQ=255)
                                                                                                         gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGAc                                               <  2:1668009/117‑1 (MQ=255)
                                                                                                          aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGaa                                                                                                   <  1:2006320/64‑1 (MQ=255)
                                                                                                          aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGaa                                                                                                   >  2:2006320/1‑64 (MQ=255)
                                                                                                          aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTa                                                                               <  1:648611/84‑1 (MQ=255)
                                                                                                          aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTa                                                                               >  2:648611/1‑84 (MQ=255)
                                                                                                             ccGTGTCTGGTTTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                              <  1:1704490/98‑1 (MQ=255)
                                                                                                             ccGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                              >  2:1704490/1‑98 (MQ=255)
                                                                                                             ccGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGAt                                                  <  2:1200512/110‑1 (MQ=255)
                                                                                                             ccGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGAt                                                  >  1:1200512/1‑110 (MQ=255)
                                                                                                             ccGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCa                     <  2:199644/139‑1 (MQ=255)
                                                                                                             ccGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCa                     >  1:2005001/1‑139 (MQ=255)
                                                                                                                tgtCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAGTGCCCAGTTACTGGGCGATTTCAGTa                  >  1:27645/1‑139 (MQ=255)
                                                                                                                tgtCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGAt                         >  1:969730/1‑132 (MQ=255)
                                                                                                                 gtCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTAt                 <  1:134521/139‑1 (MQ=255)
                                                                                                                 gtCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTAt                 >  1:419949/1‑139 (MQ=255)
                                                                                                                   cTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAGTGCCCAGTTACTGGGCGATTTCAGTATCg               <  2:27645/139‑1 (MQ=255)
                                                                                                                      gtgtAGACCAGGCACTGGATTGCCCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTACAGTATCGATa            >  1:1022240/1‑139 (MQ=255)
                                                                                                                         tAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACt         <  2:250692/139‑1 (MQ=255)
                                                                                                                            aCCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTgg                                                    <  1:405853/93‑1 (MQ=255)
                                                                                                                            aCCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTgg                                                    >  2:405853/1‑93 (MQ=255)
                                                                                                                            aCCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACt         <  1:49581/136‑1 (MQ=255)
                                                                                                                            aCCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACt         >  2:49581/1‑136 (MQ=255)
                                                                                                                            aCCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTAcc      >  2:209161/1‑139 (MQ=255)
                                                                                                                              cAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAg    <  2:163869/139‑1 (MQ=255)
                                                                                                                               aGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAatat                                                                                                >  2:1428769/1‑46 (MQ=255)
                                                                                                                               aGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAatat                                                                                                <  1:1428769/46‑1 (MQ=255)
                                                                                                                                ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTa                                                                                         >  1:1139643/1‑52 (MQ=255)
                                                                                                                                ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTa                                                                                         <  2:1139643/52‑1 (MQ=255)
                                                                                                                                ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGAt                                                  <  1:1666978/91‑1 (MQ=255)
                                                                                                                                ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGAt                                                  >  2:1666978/1‑91 (MQ=255)
                                                                                                                                ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGAt             >  2:1494924/1‑128 (MQ=255)
                                                                                                                                ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGAt             <  1:1494924/128‑1 (MQ=255)
                                                                                                                                ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACt         >  2:1034099/1‑132 (MQ=255)
                                                                                                                                ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACt         <  1:1034099/132‑1 (MQ=255)
                                                                                                                                ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGtt  >  1:1378096/1‑139 (MQ=255)
                                                                                                                                ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGtt  >  2:423281/1‑139 (MQ=255)
                                                                                                                                 gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCAc                                                       >  2:2041454/1‑85 (MQ=255)
                                                                                                                                 gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGAt                                                  <  2:1609617/90‑1 (MQ=255)
                                                                                                                                 gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGAt                                                  >  1:1609617/1‑90 (MQ=255)
                                                                                                                                 gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCaa                                            <  1:1112274/96‑1 (MQ=255)
                                                                                                                                 gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCaa                                            >  2:1112274/1‑96 (MQ=255)
                                                                                                                                 gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATg                                          >  2:634068/1‑98 (MQ=255)
                                                                                                                                 gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGc                            >  2:1350500/1‑112 (MQ=255)
                                                                                                                                 gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGc                            <  1:1350500/112‑1 (MQ=255)
                                                                                                                                 gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGt                   <  1:693651/121‑1 (MQ=255)
                                                                                                                                 gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGt                   >  2:693651/1‑121 (MQ=255)
                                                                                                                                         |                                                                                                                                  
AGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTT  >  NC_000913/4447730‑4447997

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: