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% 299.6 / NA 73coding (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 . (41/32);  total (41/32)

CAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTG  >  NC_000913/4447729‑4447998
                                                                                                                                          |                                                                                                                                   
cAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGAt                                                                                                                                     >  1:950758/1‑139 (MQ=255)
  gCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGAttg                                                                                                                                   <  2:263863/139‑2 (MQ=255)
   cGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGAt                                                                                                                                     >  1:1663977/1‑136 (MQ=255)
   cGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGAt                                                                                                                                     <  2:1663977/136‑1 (MQ=255)
    ggCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGAttgct                                                                                                                                 <  2:727383/139‑4 (MQ=255)
             ttGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGtt                                                                                                                       <  2:552088/139‑1 (MQ=255)
                    gTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGttt                                                                                                                >  1:1251481/1‑139 (MQ=255)
                     tCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGtttt                                                                                                               >  2:575341/1‑139 (MQ=255)
                     tCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGtttt                                                                                                               >  1:440248/1‑139 (MQ=255)
                          aaGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACaa                                                                                                          >  2:1211581/1‑139 (MQ=255)
                                       cGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATTAGTTCGCGTTTTAACAATGCGAATATTTGt                                                                                             >  2:877934/1‑139 (MQ=255)
                                       cGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGt                                                                                             >  1:1612854/1‑139 (MQ=255)
                                          gggTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTa                                                                                          <  2:1251481/139‑1 (MQ=255)
                                               tatTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCa                                                                                     >  1:1645904/1‑139 (MQ=255)
                                                      ctctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACg                                                                              >  1:1238049/1‑139 (MQ=255)
                                                                   cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGAttgc                                                                                                                                  >  2:550928/1‑73 (MQ=255)
                                                                   cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGAttgc                                                                                                                                  <  1:550928/75‑3 (MQ=255)
                                                                   cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGc                                                                 >  1:569093/1‑139 (MQ=255)
                                                                   cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGc                                                                 >  2:1731057/1‑139 (MQ=255)
                                                                   cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTCGCCACAAACAAGGCa                                                                >  2:13764/1‑139 (MQ=255)
                                                                   cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGAGTAGACCAGGCACTGGATTGCTCTATCAGATCGAGTTTTAACAATGCGAATATTTGTCGACGGCAGTTTACGCCACAAACAAGGc                                                                 >  1:1508980/1‑139 (MQ=255)
                                                                     cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                               <  1:965671/139‑1 (MQ=255)
                                                                     cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                               <  2:950758/139‑1 (MQ=255)
                                                                     cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                               >  2:1605953/1‑139 (MQ=255)
                                                                     cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                               <  2:1645904/139‑1 (MQ=255)
                                                                     cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                               <  2:891449/139‑1 (MQ=255)
                                                                     cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTAGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                               >  1:1497926/1‑139 (MQ=255)
                                                                          tacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTc                                                          <  2:1612854/139‑1 (MQ=255)
                                                                          tacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTc                                                          <  2:425707/139‑1 (MQ=255)
                                                                              ttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAatat                                                                                                 >  2:1465942/1‑96 (MQ=255)
                                                                              ttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAatat                                                                                                 <  1:1465942/96‑1 (MQ=255)
                                                                                aTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCATTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGa                                                    >  1:431145/1‑139 (MQ=255)
                                                                                aTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACgg                                                                                       <  1:1039252/104‑1 (MQ=255)
                                                                                aTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACgg                                                                                       >  2:1039252/1‑104 (MQ=255)
                                                                                 tCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGAt                                                   <  1:1211581/139‑1 (MQ=255)
                                                                                    ccTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGAc                                                <  1:1731057/139‑1 (MQ=255)
                                                                                      tGATGCCTAAGCTATATCTGGAAGCCGTCTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTCTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCa                                              >  1:222913/1‑139 (MQ=255)
                                                                                       gATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCaa                                             >  2:1360843/1‑139 (MQ=255)
                                                                                         tGCCTAAGCTATATCTGGAAGCCGTGTCTAGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATg                                           <  2:1497926/139‑1 (MQ=255)
                                                                                           ccTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCTAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGcc                                         <  1:1284278/139‑1 (MQ=255)
                                                                                           actagacTATATTTGGAACCTGTGTCAAGTATAGACCAGGCACTGGATTGCTCTAGCAGGTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGcc                                         <  1:877934/133‑1 (MQ=255)
                                                                                             tAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCa                                       >  2:1334731/1‑139 (MQ=255)
                                                                                               aGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGt                                     >  1:1577200/1‑139 (MQ=255)
                                                                                                    tatCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAatat                                                                                                 >  1:1706302/1‑74 (MQ=255)
                                                                                                    tatCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAatat                                                                                                 <  2:1706302/74‑1 (MQ=255)
                                                                                                    tatCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGt                                                           <  1:574664/112‑1 (MQ=255)
                                                                                                    tatCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGt                                                           >  2:574664/1‑112 (MQ=255)
                                                                                                    tatCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTg                                >  1:303423/1‑139 (MQ=255)
                                                                                                     atCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGtt                                                                                                                       >  2:724806/1‑51 (MQ=255)
                                                                                                     atCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGtt                                                                                                                       <  1:724806/51‑1 (MQ=255)
                                                                                                     atCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCca                                                                           >  2:1164989/1‑95 (MQ=255)
                                                                                                     atCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCca                                                                           <  1:1164989/95‑1 (MQ=255)
                                                                                                          gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTaac                                                                                                            >  2:895004/1‑55 (MQ=255)
                                                                                                          gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTaac                                                                                                            <  1:895004/57‑3 (MQ=255)
                                                                                                          gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTaa                                                                                                             <  2:1564901/56‑2 (MQ=255)
                                                                                                          gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTaa                                                                                                             >  1:1564901/1‑55 (MQ=255)
                                                                                                          gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCa                                                         >  1:1008130/1‑108 (MQ=255)
                                                                                                          gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCa                                                         <  2:1008130/108‑1 (MQ=255)
                                                                                                          gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGAt                          >  1:173427/1‑139 (MQ=255)
                                                                                                          gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGAt                          >  2:47067/1‑139 (MQ=255)
                                                                                                           aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGtt                                                                                                                       >  1:776691/1‑45 (MQ=255)
                                                                                                           aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGtt                                                                                                                       <  2:776691/45‑1 (MQ=255)
                                                                                                           aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGAtt                         >  1:1592769/1‑139 (MQ=255)
                                                                                                            aGCCGTGTCTGGCGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTc                                                          <  1:313772/105‑1 (MQ=255)
                                                                                                            aGCCGTGTCTGGCGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTc                                                          >  2:313772/1‑105 (MQ=255)
                                                                                                              ccGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCa                      >  2:1447844/1‑139 (MQ=255)
                                                                                                                 tgtCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTa                   >  1:134465/1‑139 (MQ=255)
                                                                                                                  gtCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTAt                  <  1:575341/139‑1 (MQ=255)
                                                                                                                  gtCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTAt                  >  1:24216/1‑139 (MQ=255)
                                                                                                                     tGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGa               <  2:569093/139‑1 (MQ=255)
                                                                                                                          tAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACt          <  2:134465/139‑1 (MQ=255)
                                                                                                                            gACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCAATTTCAGTATCGATAACTAc        <  2:1341471/139‑1 (MQ=255)
                                                                                                                             aCCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTCGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCa      <  1:13764/139‑1 (MQ=255)
                                                                                                                               cAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAg                                                                   <  2:897519/77‑1 (MQ=255)
                                                                                                                               cAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAg     <  2:1577200/139‑1 (MQ=255)
                                                                                                                               cAGGCACTGGATTGCTCTATCAGTTCGAGCTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAg                                                                   >  1:897519/1‑77 (MQ=255)
                                                                                                                                 ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGtt   >  1:968815/1‑139 (MQ=255)
                                                                                                                                  gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCAc                                                        >  1:362809/1‑85 (MQ=255)
                                                                                                                                  gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCAc                                                        <  2:362809/85‑1 (MQ=255)
                                                                                                                                  gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTg  >  2:1424657/1‑139 (MQ=255)
                                                                                                                                          |                                                                                                                                   
CAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTG  >  NC_000913/4447729‑4447998

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: