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% 359.7 / NA 87coding (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 . (45/42);  total (45/42)

TTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTA  >  NC_000913/4447742‑4448000
                                                                                                                             |                                                                                                                                     
ttGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGtt                                                                                                                         <  1:726518/139‑1 (MQ=255)
      ggTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGAt                                                                                                                                       <  2:1272126/120‑1 (MQ=255)
      ggTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGAt                                                                                                                                       >  1:1272126/1‑120 (MQ=255)
      ggTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGtt                                                                                                                   >  2:194028/1‑139 (MQ=255)
      ggTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGtt                                                                                                                   >  1:632685/1‑139 (MQ=255)
        tCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGAt                                                                                                                                       >  1:578058/1‑118 (MQ=255)
        tCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGAt                                                                                                                                       <  2:578058/118‑1 (MQ=255)
              aGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAAt                                                                                                           <  2:647174/139‑1 (MQ=255)
                      aCGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAatat                                                                                                   <  2:278610/139‑1 (MQ=255)
                         gCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATTCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTaac                                                                                                              >  2:1202490/1‑123 (MQ=255)
                         gCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTaac                                                                                                              <  1:1202490/125‑3 (MQ=255)
                            tGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCt                                                                                             >  2:47455/1‑139 (MQ=255)
                                      tGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGttt                                                                                   <  1:194028/139‑1 (MQ=255)
                                        aCTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGtt                                                                                                                         >  2:113667/1‑99 (MQ=255)
                                        aCTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGtt                                                                                                                         <  1:113667/99‑1 (MQ=255)
                                         ctctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGtt                                                                                                                         <  2:1137412/98‑1 (MQ=255)
                                         ctctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGtt                                                                                                                         >  1:1137412/1‑98 (MQ=255)
                                         ctctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACg                                                                                >  2:1284331/1‑139 (MQ=255)
                                            tataGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCca                                                                             >  1:1244222/1‑139 (MQ=255)
                                              taGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCcaca                                                                           >  2:1603781/1‑139 (MQ=255)
                                              taGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCcaca                                                                           <  2:794458/139‑1 (MQ=255)
                                              taGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCcaca                                                                           >  1:63007/1‑139 (MQ=255)
                                                      cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTaac                                                                                                              <  1:1496408/96‑3 (MQ=255)
                                                      cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTaac                                                                                                              >  2:1496408/1‑94 (MQ=255)
                                                      cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGcc                                                                              <  2:1175241/128‑1 (MQ=255)
                                                      cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGcc                                                                              >  1:1175241/1‑128 (MQ=255)
                                                      cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGc                                                                   >  1:1020600/1‑139 (MQ=255)
                                                      cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGc                                                                   >  1:952749/1‑139 (MQ=255)
                                                      cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGc                                                                   >  1:296054/1‑139 (MQ=255)
                                                         acgttacgttATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAg                                                                                                                     <  2:1214313/86‑1 (MQ=255)
                                                         acgttacgttATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAg                                                                                                                     >  1:1214313/1‑86 (MQ=255)
                                                         acgttacgttATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAAc                                                                >  2:1547521/1‑139 (MQ=255)
                                                          cgttacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACa                                                               <  2:296054/139‑1 (MQ=255)
                                                             tacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTc                                                            <  2:1471635/139‑1 (MQ=255)
                                                               cgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCAc                                                          >  2:1340051/1‑139 (MQ=255)
                                                                   aTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGa                                                      <  2:952749/139‑1 (MQ=255)
                                                                    tCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                 >  1:427767/1‑127 (MQ=255)
                                                                    tCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                 <  2:427767/127‑1 (MQ=255)
                                                                    tCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGAt                                                     <  2:1132502/139‑1 (MQ=255)
                                                                       ccTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGAc                                                  >  1:1581731/1‑139 (MQ=255)
                                                                       ccTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGAc                                                  >  1:1577142/1‑139 (MQ=255)
                                                                          gATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTTTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCaa                                               <  1:47455/139‑1 (MQ=255)
                                                                               cTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGccc                                          <  2:1244222/139‑1 (MQ=255)
                                                                                  aGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGttt                                                                                   >  2:72354/1‑95 (MQ=255)
                                                                                  aGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGttt                                                                                   <  1:72354/95‑1 (MQ=255)
                                                                                        atCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGAt                                                                                                                                       >  2:1024136/1‑38 (MQ=255)
                                                                                        atCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGAt                                                                                                                                       <  1:1024136/38‑1 (MQ=255)
                                                                                        atCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGcc                                                                              <  2:316501/94‑1 (MQ=255)
                                                                                        atCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGcc                                                                              >  1:316501/1‑94 (MQ=255)
                                                                                        atCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTgg                                 >  1:142180/1‑139 (MQ=255)
                                                                                          cTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGc                               <  2:142180/139‑1 (MQ=255)
                                                                                          cTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGc                               <  2:63007/139‑1 (MQ=255)
                                                                                            ggAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGa                             >  1:1228895/1‑139 (MQ=255)
                                                                                             gAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGcc                                           >  2:1211661/1‑124 (MQ=255)
                                                                                             gAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGcc                                           <  1:1211661/124‑1 (MQ=255)
                                                                                             gAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGAt                            >  2:1434756/1‑139 (MQ=255)
                                                                                             gAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGAt                            <  2:1415162/139‑1 (MQ=255)
                                                                                             gAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGAt                            >  2:1284979/1‑139 (MQ=255)
                                                                                              aaGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaa                                                                         >  1:1494114/1‑93 (MQ=255)
                                                                                              aaGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaa                                                                         <  2:1494114/93‑1 (MQ=255)
                                                                                              aaGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGAt                                                     <  1:854234/113‑1 (MQ=255)
                                                                                              aaGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGAt                                                     >  2:854234/1‑113 (MQ=255)
                                                                                              aaGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGAtt                           >  1:654360/1‑139 (MQ=255)
                                                                                                 ccGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGAcc                                                 >  2:1219026/1‑114 (MQ=255)
                                                                                                 ccGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGAcc                                                 <  1:1219026/114‑1 (MQ=255)
                                                                                                   gtgtCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGt                      >  1:562022/1‑139 (MQ=255)
                                                                                                    tgtCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTa                     >  1:1562396/1‑139 (MQ=255)
                                                                                                     gtCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTAt                    <  2:1228895/139‑1 (MQ=255)
                                                                                                        tGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGa                 <  1:1434756/139‑1 (MQ=255)
                                                                                                         ggTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAACGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGAt                <  1:851685/139‑1 (MQ=255)
                                                                                                          gtgtAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATa               <  2:632685/139‑1 (MQ=255)
                                                                                                             tAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACt            <  1:1603781/139‑1 (MQ=255)
                                                                                                             tAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACt            <  1:1408387/139‑1 (MQ=255)
                                                                                                             tAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACt            <  1:1547521/139‑1 (MQ=255)
                                                                                                              aGACTAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGc                               <  2:898885/119‑1 (MQ=255)
                                                                                                              aGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGc                               >  1:898885/1‑119 (MQ=255)
                                                                                                                aCCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCa                        >  2:39966/1‑124 (MQ=255)
                                                                                                                aCCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCa                        <  1:39966/124‑1 (MQ=255)
                                                                                                                aCCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTAcc         >  2:603835/1‑139 (MQ=255)
                                                                                                                 ccAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCa        <  2:1020600/139‑1 (MQ=255)
                                                                                                                  cAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAg       >  2:121318/1‑139 (MQ=255)
                                                                                                                    ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                 <  2:1176167/79‑1 (MQ=255)
                                                                                                                    ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                 >  1:1176167/1‑79 (MQ=255)
                                                                                                                    ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTGACTGGATGACCAATGcc                                           >  2:858405/1‑101 (MQ=255)
                                                                                                                    ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGcc                                           <  1:858405/101‑1 (MQ=255)
                                                                                                                     gCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaaca                                                                       >  1:727884/1‑72 (MQ=255)
                                                                                                                     gCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaaca                                                                       <  2:727884/72‑1 (MQ=255)
                                                                                                                     gCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGc                               <  2:579532/112‑1 (MQ=255)
                                                                                                                     gCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGc                               >  1:579532/1‑112 (MQ=255)
                                                                                                                     gCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCa                        >  1:217208/1‑119 (MQ=255)
                                                                                                                     gCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCa                        <  2:217208/119‑1 (MQ=255)
                                                                                                                     gCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTg    >  1:802440/1‑139 (MQ=255)
                                                                                                                       aCTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGta  >  2:349331/1‑139 (MQ=255)
                                                                                                                             |                                                                                                                                     
TTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTA  >  NC_000913/4447742‑4448000

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