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.100.0% 469.6 / NA 111coding (3756/3780 nt)tamBtranslocation and assembly module for autotransporter export, inner membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (59/52);  total (59/52)

Only 100 of 115 total aligned reads displayed.

ATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTA  >  NC_000913/4447737‑4448000
                                                                                                                                  |                                                                                                                                     
aTGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTAt                                                                                                                              >  1:662416/1‑139 (MQ=255)
        ccAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGtt                                                                                                                         <  1:1258148/136‑1 (MQ=255)
        ccAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGtt                                                                                                                         >  2:1258148/1‑136 (MQ=255)
           ggTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGtt                                                                                                                   >  1:126937/1‑139 (MQ=255)
            gTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCATGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGttt                                                                                                                  >  1:567237/1‑139 (MQ=255)
            gTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGttt                                                                                                                  >  2:898918/1‑139 (MQ=255)
             tCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGtttt                                                                                                                 >  1:1748655/1‑139 (MQ=255)
               tGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTaa                                                                                                               <  1:898918/139‑2 (MQ=255)
                 cAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCATGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACa                                                                                                             <  2:567237/139‑1 (MQ=255)
                      gAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCg                                                                                                        <  2:1748655/139‑1 (MQ=255)
                              gCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTAt                                                                                                                              >  2:1385133/1‑109 (MQ=255)
                                 tGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCt                                                                                             >  2:843601/1‑139 (MQ=255)
                                    gTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACg                                                                                          <  1:1712278/139‑1 (MQ=255)
                                             aCTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGtt                                                                                                                         >  1:50033/1‑99 (MQ=255)
                                             aCTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGtt                                                                                                                         <  2:50033/99‑1 (MQ=255)
                                                 tataGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACgg                                                                                         >  1:411463/1‑127 (MQ=255)
                                                 tataGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACgg                                                                                         <  2:411463/127‑1 (MQ=255)
                                                 tataGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCca                                                                             >  2:1191189/1‑139 (MQ=255)
                                                   taGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCcaca                                                                           >  2:1701055/1‑139 (MQ=255)
                                                   taGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCcaca                                                                           >  1:218335/1‑139 (MQ=255)
                                                        acacTCATGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGtt                                                                                                                         >  1:966308/1‑88 (MQ=255)
                                                           cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGtt                                                                                                                         <  1:748929/85‑1 (MQ=255)
                                                           cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGtt                                                                                                                         >  2:748929/1‑85 (MQ=255)
                                                           cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACgg                                                                                         <  1:1275749/117‑1 (MQ=255)
                                                           cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACgg                                                                                         >  2:1275749/1‑117 (MQ=255)
                                                           cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGc                                                                   >  2:1426619/1‑139 (MQ=255)
                                                           cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGc                                                                   >  1:307224/1‑139 (MQ=255)
                                                           cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGc                                                                   >  1:1331569/1‑139 (MQ=255)
                                                           cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGc                                                                   >  1:1004089/1‑139 (MQ=255)
                                                             cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                 <  2:1678931/139‑1 (MQ=255)
                                                             cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                 <  2:641995/139‑1 (MQ=255)
                                                             cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                 <  2:662416/139‑1 (MQ=255)
                                                                gttacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAg                                                                                                                     >  1:720790/1‑84 (MQ=255)
                                                                gttacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAg                                                              >  2:525183/1‑139 (MQ=255)
                                                                    cgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAAt                                                                                                           <  1:1751302/90‑1 (MQ=255)
                                                                    cgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAAt                                                                                                           >  2:1751302/1‑90 (MQ=255)
                                                                    cgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCAc                                                          >  1:1437972/1‑139 (MQ=255)
                                                                    cgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCAc                                                          >  1:666792/1‑139 (MQ=255)
                                                                    cgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTG‑‑GTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACgg                                                                                         >  1:1075119/1‑106 (MQ=255)
                                                                    cgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTG‑‑GTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACgg                                                                                         <  2:1075119/106‑1 (MQ=255)
                                                                     gttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGc                                                                                                         <  1:109617/91‑1 (MQ=255)
                                                                     gttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGc                                                                                                         >  2:109617/1‑91 (MQ=255)
                                                                        aTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGa                                                      <  2:1331569/139‑1 (MQ=255)
                                                                         tCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGAt                                                     <  1:525183/139‑1 (MQ=255)
                                                                         tCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGAt                                                     <  2:126937/139‑1 (MQ=255)
                                                                          cGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGtt                                                                                                                         <  1:1119631/70‑1 (MQ=255)
                                                                          cGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGtt                                                                                                                         >  2:1119631/1‑70 (MQ=255)
                                                                          cGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAGTGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATg                                                    >  1:1812620/1‑139 (MQ=255)
                                                                               gATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCaa                                               <  1:1426619/139‑1 (MQ=255)
                                                                               gATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCaa                                               <  2:666792/139‑1 (MQ=255)
                                                                                aTGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGAttg                                                                                                                                     >  1:1211436/1‑52 (MQ=255)
                                                                                aTGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGAttg                                                                                                                                     <  2:1211436/53‑2 (MQ=255)
                                                                                aTGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGAttg                                                                                                                                     <  2:470134/53‑2 (MQ=255)
                                                                                aTGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGAttg                                                                                                                                     >  1:470134/1‑52 (MQ=255)
                                                                                aTGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGtt                                                                                                                         >  1:1716645/1‑64 (MQ=255)
                                                                                   ccTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGcc                                           <  1:1701055/139‑1 (MQ=255)
                                                                                    cTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGccc                                          <  2:665369/139‑1 (MQ=255)
                                                                                             atCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGcc                                                                              >  1:1264126/1‑94 (MQ=255)
                                                                                             atCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCaa                                               <  2:411680/125‑1 (MQ=255)
                                                                                             atCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTgg                                 >  1:36813/1‑139 (MQ=255)
                                                                                             atCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTgg                                 <  2:307224/139‑1 (MQ=255)
                                                                                                 ggAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                 <  2:409832/103‑1 (MQ=255)
                                                                                                 ggAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                 <  1:1139023/103‑1 (MQ=255)
                                                                                                 ggAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                 >  1:409832/1‑103 (MQ=255)
                                                                                                 ggAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGa                             >  1:1551212/1‑139 (MQ=255)
                                                                                                  gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGc                                                                                                         <  2:1250641/62‑1 (MQ=255)
                                                                                                  gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGc                                                                                                         >  1:1250641/1‑62 (MQ=255)
                                                                                                  gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGc                                                                   <  2:268787/100‑1 (MQ=255)
                                                                                                   aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                 >  2:676885/1‑101 (MQ=255)
                                                                                                   aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                 <  1:676885/101‑1 (MQ=255)
                                                                                                   aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTc                                                            <  2:1531610/106‑1 (MQ=255)
                                                                                                   aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTc                                                            >  1:1531610/1‑106 (MQ=255)
                                                                                                   aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCaa                                               <  1:961916/119‑1 (MQ=255)
                                                                                                   aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCaa                                               >  2:961916/1‑119 (MQ=255)
                                                                                                   aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGAtt                           >  2:639557/1‑139 (MQ=255)
                                                                                                      ccGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCa                        >  2:1195910/1‑139 (MQ=255)
                                                                                                        gtgtCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGt                      >  2:1329261/1‑139 (MQ=255)
                                                                                                         tgtCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGt                                                             >  1:313962/1‑99 (MQ=255)
                                                                                                         tgtCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGt                                                             <  2:313962/99‑1 (MQ=255)
                                                                                                         tgtCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTa                     >  2:1656347/1‑139 (MQ=255)
                                                                                                                  tAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACt            <  2:1677453/139‑1 (MQ=255)
                                                                                                                    gACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGa                                                   <  1:749843/98‑1 (MQ=255)
                                                                                                                    gACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGa                                                   >  2:749843/1‑98 (MQ=255)
                                                                                                                       cAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAg       <  1:1191189/139‑1 (MQ=255)
                                                                                                                         ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACg                                                                                >  1:1604240/1‑64 (MQ=255)
                                                                                                                         ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACg                                                                                <  2:1604240/64‑1 (MQ=255)
                                                                                                                         ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGAcc                                                 >  1:1562322/1‑95 (MQ=255)
                                                                                                                         ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGAcc                                                 <  2:1562322/95‑1 (MQ=255)
                                                                                                                         ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCaa                                               <  1:711977/97‑1 (MQ=255)
                                                                                                                         ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCaa                                               >  2:711977/1‑97 (MQ=255)
                                                                                                                         ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGcc                                           <  1:364504/101‑1 (MQ=255)
                                                                                                                         ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGtt     >  2:1553387/1‑139 (MQ=255)
                                                                                                                          gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAAt                                                                                                           <  1:33050/36‑1 (MQ=255)
                                                                                                                          gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAAt                                                                                                           >  2:33050/1‑36 (MQ=255)
                                                                                                                          gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaaca                                                                       >  2:1176778/1‑72 (MQ=255)
                                                                                                                          gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaaca                                                                       <  1:1176778/72‑1 (MQ=255)
                                                                                                                          gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                 >  1:954118/1‑78 (MQ=255)
                                                                                                                          gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                 <  2:954118/78‑1 (MQ=255)
                                                                                                                            aCTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGta  <  1:639557/139‑1 (MQ=255)
                                                                                                                            aCTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGta  <  1:628792/139‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                     
ATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTA  >  NC_000913/4447737‑4448000

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